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The  base-motion  isolation  in  navigation,  guidance,  and  gunfire 
control  systems  is  a relatively  new  practice  in  the  history  of 
technology. 

In  nature,  however,  all  vertebrates,  which  include  Homo  sapiens, 
possess  built-in  base-motion  isolation  systems  in  the  form  of 
vestibuloocular  reflexes  to  stabilize  tlie  visual  axis  despite  head 
rotations.  This  reflex  was  developed  early  in  evolution,  literally 
millions  of  years  ago,  as  a means  of  survival  to  obtain  food  and  to 
protect  against  danger  by  allowing  enough  time  for  the  visual  obser- 
vation of  a stationary  environment  during  rotatory  head  motion. 

This  report  presents  the  mathematical  model  of  the  vestibuloocular 
reflex  wliich  is  formulated  from  the  viewpoint  of  cybernetics.  It  also 
presents  the  model's  simulation  results  and  discusses  the  comparison 
with  pliysiological  data.  The  model  is  constructed  in  such  a way  that 
the  signals  from  it  are  compatible  with  the  physiological  data  at  the 
system  output  levels  Moreover,  the  model's  internal  signals  are 
similar  to  the  neural  discharge  patterns  observed  at  the  corresponding 
locations  of  the  central  nervous  system. 

The  mathematical  tools  and  the  approaches  used  in  the  formulation  of 
the  model  are  similar  to  those  used  in  the  modeling  of  the  weapon  systems. 
Particularly,  the  way  the  noise  problem  is  treated  in  the  model  may  be 
useful  in  simulations  of  certain  weapon  systems. 
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CHAPTER  I 


INTRODUCTION 

1 . 1 An  elementary  description  of  the  vestibuloocular  reflex 
The  function  of  the  vestibuloocular  reflex  'VOR)  is  to 
stabilize  the  visual  axis  in  the  environment  during  a head- 
rotation  by  compensating  for  it.  It  does  this  by  making  an 
equal  and  opposite  eye-rotation  in  the  orbit  with  the  same 
time  course  as  the  head  movement. 

One  way  to  demonstrate  the  speed  and  accuracy  of  the 
VOR  is  to  try  to  read  a book  held  stationary  in  one's  hands 
while  one's  head  is  oscillating  at  about  one  or  two  Hz. 

Next,  reverse  the  procedure  and  try  to  read  a book  being 
oscillated  by  one's  hand  at  similar  angular  amplitudes  and 
frequencies  while  one's  head  is  held  stationary.  In  the 
latter  case,  because  the  eye  cannot  keep  up  with  the  book 
which  causes  its  image  to  slip  on  the  retina,  visual  acuity 
is  so  low  that  reading  is  virtually  impossible.  In  the 
former  case,  visual  acuity  is  good  (normal)  because  the 
visual  axes  are  being  automatically  stabilized  with  respect 
to  the  stationary  book  by  the  VOR. 

The  VOR  isolates  the  eyeball  from  the  disturbances  of 
the  head's  motions  by  an  open-loop  reflex  action,  unlike 
some  mechanical  control  systems  which  isolate  a so-called 
stable  platform  from  the  disturbances  of  the  base  motions 
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of  ships  and  rockets  by  closed-loop,  or  negative  feedback, 
control  systems. 

The  reflex  is  phylogenetically  very  old.  It  developed 
in  the  early  stages  of  the  evolution  of  vertebrates  as  a 
means  of  survival  by  stabilizing  vision  of  the  environment 
in  spite  of  head  and  body  motions.  The  reflex  is  anatomically 
and  physiologically  similar  for  widely  different  species. 

This  allows  a fairly  direct  inference  from  animal  experiments 
to  man. 

For  most  normal  head  movements,  the  required  eye  move- 
ment is  small  (e.g.,  20  30°)  and  well  within  the  oculo- 

motor range.  However,  because  of  mechanical  limitations  of 
eye  rotation,  if  head  rotation  is  continued  beyond  a certain 
limit,  the  eye  must  be  reset.  This  is  done  by  a rapid  move- 
ment to  a new  position  in  the  direction  of  the  head  motion. 

If  the  head  is  continuously  rotated  in  one  direction  (say 
to  the  left)  with  a step  of  constant  angular  velocity  H 
(Figure  1.1a),  the  eye  position  0 as  a function  of  time  t 
(Figure  1.1b)  will  describe  a rhythmic  saw-tooth  motion  with 
alternate  periods  of  low  velocity  in  the  direction  opposite 
to  the  head  motion,  thus  compensating  the  head  motion,  and 
fast  velocity  in  the  direction  of  the  head  motion,  thus 
catching  up  with  the  head  by  rapid  resetting  movements. 

This  typo  of  eye  motion  is  called  rotatory  or  vestibular 
nystagmus.  The  compensatory  slow  eye  motion  is  called  the 
slow-phase  (SP)  and  the  quick  returning  motion  is  called 
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Figure  1.1 


For  step  input  head  velocity  H (a) , 
the  eye  position  0(t)  describes  a 
rhythmic  saw-tooth  motion  (b) , with 
alternate  slow  and  quick  eye  veloc- 
ities 9 in  opposite  directions  (c) . 


the  quick -phase  (QP)  of  the  nystagmus.  The  profile  of  the 
eye  velocity  0 during  rotatory  nystagmus  is  shown  in  Fig- 
ure 1.1c  to  indicate  that  during  each  slow-phase  the  eye 
velocity  is  equal  and  opposite  to  head  velocity. 

The  fact  that  the  mean  eye  position  (Figure  1.1b)  shifts 
in  the  direction  of  the  QP's  is  one  of  the  major  points  of 
investigation  in  this  study  and  will  be  discussed  in  detail 
subsequently . 

This  reflex  is  designed  to  assist  vision  but  is  suffi- 
ciently automatic  that  it  works  almost  as  well  without  it. 

To  eliminate  the  modifications  introduced  by  vision,  the  VOR 
is  usually  studied  in  isolation  by  rotating  subjects  in  the 
dark.  This  thesis  concerns  only  the  latter  situation. 

1.2  A brief  sketch  of  vestibuloocular  physiology 

There  are  two  groups  of  semicircular  canals  (SCC)  on 
each  side  of  the  skull  located  within  the  bony  labyrinth  of 
the  inner  ear.  Each  group  consists  of  three  approximately 
mutually  orthogonal  canals.  The  three  canals  on  each  side 
extend  from  a common  compartment  called  the  utricle  which  is 
contained  in  the  egg-shaped  "vestibule"  of  the  inner  ear. 

The  expression  vestibuloocular  reflex  is  derived  Irom  the 
fact  that  the  reflex  is  created  by  the  action  of  the  vesti- 
bular organ  upon  the  oculomotor  system. 
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For  horizontal  head  rotation  (yaw)  in  the  normal  upright 
position  the  lateral  semicircular  canal  pair  (one  on  each 
side  of  the  head)  is  primarily  stimulated.  The  experimental 
data  on  which  this  research  is  based  are  obtained  by  rotatory 
tests  in  the  horizontal  plane  of  subjects  in  the  normal 
upright  position.  For  this  reason,  we  will  consider  for 
simplicity  that  only  the  lateral  semicircular  canals  are 
stimulated  in  this  study.  There  is  no  reason  to  suppose 
that  the  data  cannot  be  simply  extrapolated  to  the  other 
canals  for  rotations  in  the  roll  and  pitch  directions. 

The  semicircular  canal  (Figure  1.2)  is  a circular, 
closed,  membranous  tube  encased  inside  rigid  Ixane.  The 
tube's  internal  radius  is  about  0.15  mm  for  man  (Melvill 
Jones  and  Milsum,  1965) . The  membranous  tube  is  filled 
with  a fluid  called  endolymph.  The  membranous  labyrinth  is 
submerged  in  a clear  fluid  called  perilymph  contained  in 
the  bony  labyrinth  in  the  petrous  portion  of  each  temporal 
bone. 

Each  semicircular  canal  completes  a hydrodynamic 
circuit  through  the  utricle  to  which  the  three  canals  on 
each  side  of  the  head  are  joined.  Each  canal  completes 
about  two-thirds  of  a circle  and  is  interrupted  by  an 
enlarged  specialized  compartment  called  the  ampulla  located 
adjacent  to  the  utricle.  The  ampulla  contains  a transverse 
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Figure  1.2 

The  figure  shows  a schematic  repre- 
sentation of  the  semicircular  canal 
SSC,  which  is  a circular,  closed, 
membranous  tube  encased  inside  rigid 
bone.  Head  angular  acceleration 
causes  the  endolymph  to  flow.  This 
causes  a deflection  of  the  cupula, 
which  in  turn  creates  a signal  which 
passes  into  the  central  nervous 
system. 
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crest  called  the  crista  ampullaris  which  contains  a sensory 
epithelium  consisting  of  hair  cells.  The  cupula,  a mass  of 
gelatinous  material  which  receives  the  processes  of  the  hair 
cells,  covers  the  crista  ampullaris  and  forms  a fluid  tight, 
elastic  valve  across  the  canal  in  the  ampulla.  Head  angular 
acceleration  causes  the  endolymph  to  flow  (in  the  direction 
opposite  to  the  head  motion)  because  of  its  inertial  reaction 
which  causes  a deflection  of  the  cupula.  This,  in  turn, 
bends  the  hairs  and  depolarizes  the  sensory  hair  cells.  This 
modulates  the  discharge  rate  of  action  potentials  of  first- 
order  vestibular  neurons.  This  signal  then  passes  into  the 
central  nervous  system. 

When  the  two  lateral  semicircular  canals  are  stimulated 
(in  push-pull)  by  head  angular  acceleration,  the  signals  from 
them  are  t:ransmitted  to  the  vestibular  nuclei  (VN)  through 
the  non-auditory  (vestibular)  portion  of  the  Vlllth  cranial 
nerve  (Figure  1.3).  From  the  vestibular  nuclei  the  signal 
reaches  the  oci’lomotor  nucleus  III  and  abducens  nucleus  VI 
(for  simplicity  these  will  often  be  lumped  together  and 
referred  to  simply  as  the  oculomotor  nuclei  (OMN) ) through 
two  pathways,  one,  direct  (DIR)  and  another,  indirect,  through 
the  paramedian  pontine  reticular  formation  (PPRF) . The 
signals  from  the  oculomotor  nuclei  innervate  the  extraocular 
muscles  (EOM) ; the  medial  rectus  via  the  Illrd  nerve  and 


Figure  1.3 


This  figure  shows  schematically  the 
neural  signal  flow  of  the  vesti- 
buloocular  reflex  (restricted  to 
che  slow-phase,  excluding  the  quick- 
phase)  in  the  bilateral  push-pull 
arrangement.  If  we  view  the  VOR 
.as  a system,  the  input  to  the  system 
is  the  head  acceleration,  while  the 
output  from  the  system  is  the  eye 


position. 


Figure  1.3  Neural  Signal  Flow  Diagram  (bilateral  push-pull  arrangement) 


lateral  rectus  via  the  Vlth  nerve.  The  former  pulls  the  eye 
medially,  the  latter  laterally;  together,  in  push-pull,  they 
effect  horizontal  eye  movements.  There  is  no  stretch  reflex 
in  the  extraocular  muscles  and  therefore  no  local  feedback 
around  the  eye  plant  (Keller  and  Robinson,  1971).  A simpli- 
fied diagram  with  two  bilateral  push-pull  elements  (canals, 
nuclei,  fiber  tracts  and  muscles)  is  shown  in  Figure  1.3 
(Robinson,  1975). 

There  are  crossed  inhibitory  pathways  at,  at  least, 
two  levels.  Each  vestibular  nerve  inhibits  the  contralateral 
VN  (Shimazu  and  Precht,  1966)  and  each  output  of  the  brain- 
stem centers  inhibits  antagonist  motoneurons  as  well  as 
exciting  agonist  motoneurons.  The  inhibitory  signal  reinforces 
the  excitation  in  a push-pull  manner  at  each  level  because 
it  is  180  degrees  out  of  phase  with  the  excitory  signal  and 
is  subtracted  from  it.  This  is  called  disinhibition.  Thus 
each  VN  cell,  for  example,  is  excited  from  the  ipsilateral 
canal  and  disinhibited  by  the  contralateral  canal.  These 
connections  do  not  alter  the  basic  signal  processing  and  have 
been  left  out  for  simplicity. 

As  shown  in  Figure  1.3,  the  left  channel  and  the  right 
channel  of  the  VOR  are  exactly  symmetrical  to  each  other. 

Since  both  channels  have  the  same  input  (head  acce] eration) 
and  the  same  output  (eye  position) , we  have  constructed  an 
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equivalent  single  channel  diagram  (Figure  1.4)  to  simplify 
the  analysis  without  loss  of  generality.  In  Figure  1.4, 
the  equivalent  canal  output  can  go  either  positive  or  neg- 
ative and  the  equivalent  muscle  can  push  as  well  as  pull. 

1 . 3 The  statement  of  the  problem 

The  slow-phase  portion  of  the  VOR  has  been  studied  by 
memy  people,  and  there  is  a general  consensus  of  how  it  works 
and  how  it  can  be  represented  by  a mathematical  model.  This 
will  be  presented  in  more  detail  later,  but  essentially,  the 
canals  when  stimulated  by  head  angular  acceleration,  produce 
a neural  signal  proportional  to  head  velocity.  The  velocity 
signal  is  integrated  in  the  pontine  reticular  formation  to 
produce  a signal  proportional  to  head  angular  position.  This 
actuates  the  eyeball  to  move  in  the  opposite  direction  so 
that  eye  angular  position  is  equal  but  opposite  to  head 
angular  position.  The  velocity  signal  is  also  fed  directly 
to  the  eye  plcint  to  offset  the  lag  generated  by  the  visco- 
elastic nature  of  the  eye  plant.  Although  a few  details 
remain  to  be  worked  out,  and  there  are  still  problems 
associated  with  specific  neuroanatomical  connections,  the 
general  signal  processing  events  are  fairly  well  understood 
in  the  modeling  of  the  slow-phase  system. 

However,  the  modeling  of  the  quick-phase  part  is 
a different  story.  As  of  this  writing,  very  little  qu^mtitive 
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Figure  1.4 


This  figure  shows  schematically  a 
single  channel  neural  signal  flow 
of  the  vestibuloocular  reflex 
(restricted  to  the  slow-phase 
only)  equivalent  to  the  bilateral 
push-pull  arrangement  of  Figure  1.3 
This  is  done  to  simplify  the 
analysis  without  loss  of  generality 


Neural  Signal  Flow  Diagram  (equivalent 


descriptions  exist  as  to  what  determines  the  time  of  the 
quick-phase  initiation  and  the  magnitude  of  a quick-phase, 
and  no  model  has  been  presented  which  tries  to  tie  its 
elements  as  closely  as  possible  to  establish  neurophysiological 
results  and  tries  to  evolve  its  organization  on  a teleological 
basis. 

Sugie  and  Melvill  Jones  (1971)  and  Schmid  and  Lardini  (1976) 
have  proposed  VOR  models  with  quick-phases.  Their  models  are 
discussed  in  detail  in  section  8.4  and  compared  to f the  present 
model  but,  in  general,  in  one  way  or  another  they  have  proposed 
circuits  and  elements  which  are  contrary  to  current  neurophysi- 
ological findings.  More  important,  they  do  not  deal'  with  the 
stochastic  behavior  of  the  vestibuloocular  reflex  which  is  a 
major  interest  in  this  development.  The  only  stochastic 
analysis  of  which  we  are  aware  is  that  of  Cheng  (1972)  who 
used  the  technique  of  point  process  analysis  to  analyze  the 
time  intervals  between  successive  quick-phases.  He  has 
described  a simple  stochastic  model  which  could  simulate 
various  observed  inter-saccadic  interval  characteristics  of 
vestibular  and  optokinetic  nystagmus.  However,  his  analysis 
is  purely  statistical  in  nature  and  made  no  attempt  to  model 
the  VOR. 

The  goal  of  this  thesis  is  to  construct  a mathematical 
model  for  the  VOR  which  includes  quick-phases  and  which  is 
as  consistent  as  possible  with  recent  advances  in  neuro- 
physiology. First,  it  must  reasoncibly  simulate  the 
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deterministic  characteristics  of  nystagmus  and  then  include  the 
stochastic  behavior  of  the  VOR  nystagmus  by  introducing 
appropriate  random  signals  at  neurophysiologically  appropriate 
points  in  the  deterministic  model.  The  final  model  determines 
not  only  when  a quick-phase  should  be  triggered,  but  also 
determines  where  it  should  carry  the  eye. 

Experimental  data  were  obtained  from  a cat  rotated  by 
head  velocity  steps  and  sine  waves  of  various  amplitudes  and 
frequencies.  These  results  formed  the  data  base  for  the 
model.  The  topology  of  the  model  was  then  set  up  by  inference 
from  a variety  of  physiological  observations  and  the  param- 
eters adjusted  to  give  a best  match  on  a digital  computer 
simulation  for  the  deterministic  and  stochastic  behavior 
of  the  reflex  that  was  experimentally  observed. 
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CflAPTER  II 


METHODS 


2 . 1 Experimental  data 

All  experimental  data  used  in  this  research  were  obtained 
from  a single  cat.  Data  from  cat  were  readily  available 
from  an  ongoing  experiment  in  the  laboratory.  A fairly 
direct  inference  from  cat's  VOR  to  man's  can  be  made  because 
the  reflex  is  anatomically  and  physiologically  similar  for 
widely  different  species. 

The  fact  that  the  data  came  from  one  cat  instead  of 
many  cats  presents  no  disadvantage  in  this  study  because 
the  main  interest  is  in  the  topology  and  phenomenology  of 
the  VOR  in  general  and  not  in  the  quantitative  details  of 
the  VOR  of  any  particular  species.  Nevertheless,  it  was 
important  to  show  that  the  model  could  be  adjusted  to 
quantitatively  fit  actual  biological  data.  The  data  from 
one  cat  formed  a self-consistent  data  base.  The  data  from 
a large  population  of  cats  were  simply  not  of  interest. 

Two  types  of  inputs  (head  motions)  were  used.  For  step 
velocity  inputs,  nominal  head  velocities  of  approximately 
10,  20,  40,  and  60  dog/sec  were  used.  Although  some  ou*-put 
data  from  input  head  velocity  of  5 deg/sec  were  available, 
they  were  too  limited  to  be  used  for  any  firm  conclusions 
but  did  indicate  some  trends.  For  sinusoidal  inputs. 
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frequencies  of  0.05,  0.25,  and  1.2  Hz  with  peak  head  veloci- 
ties ranging  from  10  to  30  deg/ sec  were  used. 

2 . 2 Experimental  setup 

Cat's  eye  movements  are  measured  by  techniques  originally 
developed  by  Robinson  (1963)  and  further  refined  by  Fuchs 
and  Robinson  (1966).  Briefly,  in  this  method,  a coil  of  fine 
wire,  surgically  implanted  upon  the  eyeball  beneath  the  bulbar 
conjunctiva  and  recti  insertions,  is  subject  to  horizontal 
and  vertically  alternating  magnetic  fields  90°  out  of  phase 
with  each  other.  A voltage  is  induced  in  the  coil  which  is 
a function  of  the  angle  between  the  eye  position  and  magnetic 
field.  This  voltage  is  phase-detected  against  both  horizontal 
and  vertical  reference  voltages,  resulting  in  two  d-c  voltages 
proportional  to  the  horizontal  and  vertical  angles  of  eye 
position. 

The  technique  makes  it  possible  to  record  eye  movements 
conveniently  and  accurately  in  the  experimental  animal  in 
long-term  experiments  on  oculomotor  function.  The  method 
has  a resolution  of  1.5  min  of  arc,  is  linear  within  6%  of 
a full  scale  of  + 20  deg,  and  possesses  a system  bandwidth 
of  1,000  Hz, 

A platform  was  attached  to  the  cat's  skull  by  three 
stainless  steel  bolts.  The  cat  was  restrained  in  a close- 
fitting  box  and  its  head  immobilized  by  fastening  the 
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platform  to  a shelf  extending  from  the  box.  The  box  was 


placed  in  the  magnetic  field  coil  frame  which  could  rotate 
in  the  horizontal  plane.  Head  (frame)  angular  velocity  was 
measured  by  a tachometer  generator.  Eye  velocity  was  obtained 
by  differentiating  eye  position  over  the  bandwidth  0-16  Hz. 

All  experiments  were  done  in  total  darkness.  Care  was  taken 
to  keep  the  cat  alert. 

2 . 3 Theoretical  calculation 

The  mathematical  model  formulated  in  this  research  was 
programmed  in  FORTRAN  IV  language.  The  simulations  of  the 
model  were  performed  by  a CDC  6700  digital  computer.  The 
output  graphs  were  plotted  by  a Calcomp  plotter.  The 
computation  interval  was  fixed  at  0.01  sec  which  is  much 
smaller  than  quick -phase  durations  (which  are  generally  more 
than  0.1  sec  for  the  cat),  the  periods  of  input  sinusoids 
(1.2  Hz,  maximum)  and  any  system  time  constant  (which  were 
larger  than  about  0.1  sec). 

For  the  numerical  integration  of  the  differential 
equations,  all  with  the  form  of  x(t)  = a x(t)  + b u(t)  where 
a and  b are  constants,  the  time  domain  solutions  were  directly 
used.  At  one  point  a simple  version  of  the  predictor-corrector 
method  for  the  numerical  integration  was  tried.  The  two 
methods  were  virtually  indistinguishable  at  the  accuracy 
releveuit  to  our  investigations  which  is  about  0.1  degrees 
of  eye  position. 
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Gaussian  rcuidom  numbers  were  generated  sequentially 
by  the  computer  at  each  computation  interval  and  fed  into 
the  model  at  specified  points  after  passing  through  appro- 
priate filters  to  meet  the  characteristics  of  randomness 
shown  by  the  data  (for  detail,  see  section  6.3). 


CHAPTER  III 


THE  COMPENSATORY  SLOW-PHASE 

3 . 1 The  general  description 

From  the  viewpoint  of  control  systems  analysis,  the 
input  to  the  VOR  system  is  the  head  angular  acceleration 
H with  respect  to  the  environment  (specifically,  an  inertial 
reference  frame)  and  the  output  from  the  system  is  the  eye 
angular  position  6 with  respect  to  the  head.  The  receptor 
organs  of  the  VOR,  which  are  stimulated  by  the  head  acceler- 
ation, are  the  semicircular  canals.  The  effector  organ  of 
the  VOR  system,  which  produces  the  response  in  the  form  of 
eye  position,  is  the  oculomotor  plant  (OMP)  which  consists 
of  the  eyeball,  the  extraocular  muscles  and  associated 
suspensory  tissues. 

The  pathway  between  the  receptor  and  effector  organs 
consists  of  two  branches,  one  direct  from  the  vestibular 
nuclei  to  the  eye  muscle  motoneurons  and  the  other  through 
the  paramedian  pontine  reticular  forma*-ion  (PPRF)  . 

3 . 2 The  semicircular  canal* 

The  semicircular  canal  (SCO  is  considered  as  a rigid 
circular  tube  with  a very  small  uniform  cross-section.  Each 
canal  completes  a hydrodynamic  circuit  through  the  utricles 
(which  it  shares  in  common  with  the  other  two  canals  on  its 

♦References:  Melvill  Jones,  and  Milsum,  1965;  Outerbridge,  1969 
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side  of  the  head)  with  negligible  coupling  with  the  other 
canals.  The  fluid  flow  within  the  canal  is  laminar  because 
of  the  small,  smooth  bore  of  the  canal.  The  canal  has  high 
viscous  damping  as  indicated  by  the  small  Reynolds  number 
of  the  flow.  It  is  common  to  use  lumped  analysis  and  ignore 
the  flow  distribution  within  the  canal.  Thus,  the  fluid  is 
considered  to  rotate  as  a ring  relative  to  head. 

The  following  analysis  is  based  on  angular  rotation  of 
the  canal  in  a single  plane.  Consider  a counter-clockwise 
rotatioi  (defined  as  the  positive  direction  of  rotation)  of 
the  canal  and  hence  the  head  as  described  in  Figure  1.2.  The 
damping  torque  (positive  when  acting  counter-clockwise) 
on  the  fluid  ring  is  assumed  to  be 

= -b0(F/H)  (3.2.1) 

d 

where  0(F/H)  is  the  angular  velocity  of  the  fluid  relative 
to  the  head  and  b is  a positive  proportionality  constant. 

Note  that  0(F/H)  is  clockwise  relative  to  the  head  in  this 
case  because  of  inertial  reaction.  The  elastic  restoring 
torque  on  the  fluid  ring  is  assumed  to  be 

M = -k0(F/H)  (3.2.2) 

e 

where  9 (F/H)  is  the  angular  displacement  of  the  fluid  relative 
to  the  head  and  k is  a positive  proportionality  constant. 

By  Newton's  second  law, 

j'e(F/I)  = + M = -b&(F/H)  - k0(F/H)  (3.2.3) 

a e 
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in  which  J is  the  fluid  ring  moment  of  inertia  and  0(F/1) 
is  the  angular  acceleration  of  the  fluid  relative  to  inertial 
space.  By  denoting  the  angular  acceleration  of  the  head 
relative  to  space  by  0(H/I),  it  follows  that 

6(F/i)  = 0(H/1)  + 0(F/H)  (3.2.4) 

Assuming  that  the  cupular  deflection  relative  to  the  head 
0(C/H)  is  proportional  to  the  angular  displacement  of  the 
head  relative  to  fluid  0(H/F), 

0(C/H)  = a0(H/F)  = -a6(F/H)  (3.2.5) 

where  a i'^  a positive  constant,  which  reflects  the  fact  that 
the  area  of  the  canal  is  not  equal  to  the  area  of  the  ampulla. 
Substituting  (3.2.4)  and  (3.2.5)  into  (3.2.3) 

J0(C/H)  + b0(C/H)  + )c6(C/H)  = aJ0(H/I)  (3.2.6) 

Replacing  cupular  deflection  0(C/H)  by  4)  and  0(H/I)  by  H in 
(3.2.6)  , 

J<t>  h<t>  + k4>  = aJH  (3.2.7) 

which  is  identical  with  the  equation  given  by  Melvill  Jones 
and  Milsum  (1965) . It  is  the  same  form  as  the  differential 
equation  for  a torsion  pendulum  which  describes  the  torsional 
vibrations  of  an  elastic  shaft  with  a circular  rotor  rigidly 
attached  to  it.  The  canal  system  is  probably  ten  times  more 
than  critically  damped. 

Since  the  system  represented  by  (3.2.7)  is  very  over- 
damped (b^  >>  J)c) , the  transfer  function  (in  Laplace  trans- 
form notation)  of  the  SCC  may  be  approximated  by 
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(3.2.8) 


H(s)  (T  s + 1)  (T  s + 1) 
J 2 


where 


1 k 


T 

2 b 


For  man  T is  thought  to  be  about  10  ~ 16  sec  and  T about 
1 2 

0.003  sec.  For  cat  T is  found  to  be  about  4 sec  and  T 

1 2 

about  0.0015  sec  (Melvill  Jones  and  Milsum,  1971).  For  the 

frequency  range  below  5 Hz  corresponding  to  normal  daily 

activities,  where  T s < 1,  the  transfer  function  may  be 
2 

further  approximated  as 


aT  T 
1 2 


H(s)  T s + 1 
1 

Further,  in  the  range  above  about  0.05  Hz,  where  sT^  > 1 
which  includes  most  normal  head  movements. 


iii  = (aT  ) i 
IKS)  ^ 


which  is  a pure  integrator  with  gain  (aT  ) . 

2 


It  follows  that  in  this  frequency  range. 


')i(t)  “ / H(t)dt  « H(t) 
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Thus,  the  output  of  the  SCC  is  proportional  to  head  velocity 
over  the  range  of  natural  head  movements,  th.e  role  of  SCC 
being  that  of  an  integratng  accelerometer  or  velocity  trans- 
ducer. Indeed,  the  experimentally  measured  firing  rate 
modulation  of  primary  vestibular  afferents  is  proportional 
to  head  velocity  over  this  frequency  range  (Fernandez  and 
Goldberg,  1971). 

The  gain  constants  of  the  internal  neural  signal  proc- 
essing elements  are  all,  for  practical  purposes,  indeter- 
minable because  the  signals  all  consist  of  the  firing  rates 
of  large  populations  of  neurons,  only  a few  of  which  can  be 
observed  at  any  one  time.  Only  the  final  overall  gain  is 
important  and  interior  gains  may  be  adjusted  arbitrarily  so 
long  as  the  total  gain  is  kept  correct.  Thus,  the  gain  of 
the  transfer  function  (3.2.10)  of  the  SCC  may  be  arbitrarily 
adjusted  so  that  (3.2.10)  will  behave  like  a pure  integrator 
with  a unity  scale  factor  at  midband  frequencies.  That  is, 

T is  relabeled  the  cupula  long-time  constant  T and  (aT  T ) 
is  replaced  by  T^.  It  follows  from  (3.2.10)  that 

T 

4^(3)  c 

~ sT  +1  (3.2.13) 

H (s)  C 

which  reduces  to  1/s  for  sT  >1.  The  break  frequency  for 

c 

l/(sT  + 1)  is  about  0.016  Hz  with  T = 10  sec  (for  man), 
c c 

Thus,  for  all  frequencies  of  normal  head  rotation  above 
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0.016  Hz,  the  canals  integrate  head  acceleration  and  produce 

a signal  proportional  to  head  velocity.  For  cat  for  which 

T is  about  4 sec,  the  break  frequency  is  about  0.04  Hz. 
c 


3.3  The  oculomotor  plant 

The  oculomotor  plant  consists  of  the  eyeball,  the  extra- 
ocular muscles  and  associated  suspensory  tissues.  The 
relationship  between  eye  position  6 in  the  head  eind  the 
discharge  rate  of  the  oculomotor  neurons  has  been  found 
by  recording  these  variables  simultaneously  in  alert  monkeys 
(Robinson,  1'370;  Robinson  and  Keller,  1972) . To  a first 
order  approximation, 

R^  k{H-0^)  + re  (3.3.1) 

where  0^  is  the  threshold  eye  angle  at  which  the  neuron  is 

recruited  into  activity  and  k and  r are  coefficients  which 

relate  R to  0 and  0.  One  can  rearrange  the  cibove  equation 
e 

by  redefining  R as  the  change  in  rate  from  the  baseline 
value  (~k0^)  thus  dropping  this  term  out. 

Factoring  out  k, 

R = k(0  + f 6)  (3.3.2) 

e k 

This  gives  the  transfer  function  (arbitrarily  iticiking  the 
d.c.  gain  factor  unity) 


0 (s)  _ 1 

R (s)  sT  + 1 
e e 


(3.3.3) 
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where  the  time  constant  T , which  is  r/k,  has  a neuron 

e 

population  mean  of  about  0.2  sec. 

There  is  no  stretch  reflex  in  the  extraocular  muscles 
(Keller  and  Robinson,  1971) . This  means  that  there  is  no 
local  feedback  around  the  plant  which  could  conceivably 
change  its  transfer  function  for  different  eye  movement 
tasks . 


3 . 4 The  central  pathways 

From  the  viewpoint  of  control  systems  analysis,  the 
whole  purpose  of  the  VOR  is  to  recover  the  head  angular 
position  signal  by  integrating  the  head  angular  acceleration 
signal  with  respect  to  time  twice. 

The  first  integration  is  performed  by  the  SCC.  Another 
integration  is  needed  somewhere.  The  oculomotor  plant 
cannot  be  an  integrator  in  the  bandwidth  within  which  the 
VOR  works  properly  (down  to  0.02  Hz)  because  the  plant  will 
not  introduce  much  phase  lag  below  cibout  0.8  Hz  (for  T^  = 0.2). 
This  has  been  further  verified  by  rotating  monkeys  sinu- 
soidally and  measuring  the  phase  shift  between  eye  position 
and  the  sinusoidally  modulated  firing  rates  of  ocular  moto- 
neurons (Skavenski  and  Robinson,  1973) . 

The  equation  R^  = k(6  + 0.20)  in  the  last  section 
suggests  two  central  pathways  from  the  vestibular  nuclei  VN 
to  the  oculomotor  nuclei  OMN.  One,  a direct  pathway  with 
a pure  gain  of  T^  = 0.2,  furnishes  the  0.20  term  (raw 
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velocity  information)  in  the  above  equation  to  the  OMN.  The 
other  pathway,  possibly  throuqh  the  reticular  formation, 
furnishes  the  6 term  in  the  above  equation  to  the  OMN  after 
integrating  the  velocity  signal  with  respect  to  time.  Thus, 
this  path  fulfills  the  role  of  the  needed  second  integrator 
to  make  the  VOR  work.  The  indirect  pathway  will  thus  be 
referred  to  as  the  neural  integrator  NI.  The  transfer  func- 
tion of  the  direct  pathway  is  simply  its  gain  T . The 

e 

transfer  function  of  the  indirect  pathway  NI  is  assumed  to 
be  a leaky  integrator  because  it  is  unlikely  to  be  ideal 
and  there  is  more  direct  evidence  that  it  behaves  this  way. 

It  is  given  by 

T 

TF(NI)  = (3.4.1) 

n 

which  approaches  a pure  integrator  1/s  for  frequencies 
greater  than  about  0.006  Hz  because  of  its  rather  large  time 
constant  of  about  25  sec  (Becker  and  Klein,  1973) . 

Thus  for  normal  head  movements  whose  spectral  components 
lie  well  above  0.01  Hz,  the  transfer  function  for  the  central 
pathways  between  VN  and  OMN  is 

sT  + 1 

T + - = — (3.4.2) 

e s s 
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This  conveniently  reduces  the  transfer  function 
between  the  VN  and  the  eye  position  to  that  of  a pure  inte- 
grator as  shown  below  (see  Figure  3.1  in  the  next  section): 


(3.4.3) 


One  may  interpret  the  role  of  the  direct  pathway  as  creating 
the  phase  lead  needed  to  offset  the  phase  lag  that  would 
otherwise  be  created  by  the  plant. 


3 . 6 The  compensatory  s low-phase  model 

Putting  all  tliese  elements  together,  the  total  transfer 
function  of  the  compensatory  slow-phase  model  for  frequencies 
above  0.006  Hz  (corresponding  to  T^  = 25  sec)  is 


sT 

c 

sT  + 1 
c 

(3.5.1) 


For  frequ'^ncies  above  0.01  Hz  (corresponding  to  = 16  sec), 
it  further  simplifies  to 


6(s) 

H(s) 


= -1 


(3.5.2) 


Equation  (3.5.2)  meets  the  requirement  of  the  slow-phase 
VOR  which  is  to  recover,  in  the  opposite  direction,  the 
head  angular  position  signal. 
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For  cat,  the  cupula  time  constant  of  the  SCC  is 
experimentally  found  to  be  about  4 sec  (Melvill  Jones  and 
Milsum,  1971) . This  was  determined  by  directly  recording 
from  vestibular  neurons.  This  gives  a breaJc- frequency  of 
cibout  0.04  Hz  for  the  cat,  which  is  about  four  times  as 
high  as  that  of  0.01  Hz,  corresponding  to  T^  = 16  sec  for 
mem.  But  it  turns  out  that  the  VOR  works  properly  at  a 
frequency  lower  than  0.04  Hz  for  cat.  According  to  a recent 
investigation  (Robinson,  1976)  of  the  Bode  diagreun  of  the 
VOR  of  about  ten  cats  tested  in  the  dark,  the  actual  trcinsfer 
function  between  eye  position  0 (output)  and  head  position 
H (input)  is  given  by 


0(s) 

H(s) 


sT 


-(0.9) 


sT'  + 1 

» c 


(3.5.3) 


in  which  t'  - 12  sec  and  T ==50  sec.  The  factor  sT  /(sT  + 1) 
c a a/  a 

was  noticed  by  others  (e.g.,  '<"oung,  J969)  and  appears  to 
reflect  adaptation.  This  factor  cannot  be  the  neural  integrator 
(3.4.1)  for  the  following  reason.  When  quick-phase  signals 
are  added  to  the  slow-phase  signal  (see  Figure  5.1)  they 
chop  up  the  latter  and  remove  all  of  its  low  frequency 
components.  This  means  that  a poor  low  frequency  response 
of  the  NI  will  not  be  detected.  The  chopping  action  of  the 
quick-phase  pulses  makes  even  a badly  leaky  integrator  look 
perfect.  Since  the  effect  of  sT ^(sT^+1)  is  observable,  it 
must  occur  before  the  site  of  quick-phase  injection  and 
adaptation  is  the  most  likely  cause.  The  break  frequency 
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of  sT  /(sT  + 1)  is  about  0.003  Hz  which  is  well  below  the  lower 
a/  a 

frequency  limit  for  proper  functioning  of  the  reflex  and  is  also 

below  the  frequency  range  of  interest  in  this  study.  Thus,  we 

may  replace  sT  /(sT  +1)  by  unity  for  our  purpose.  Equa- 
<1/  a 

tion  (3.5.3)  now  may  be  approximated  as 


0 (s) 
H(s) 


-(0.9)s^ 


sT'  + 1 


1^ 

s 


(3.5.4) 


\ c / 

If  we  identify  s with  the  conversion  of  head  position  to 
head  acceleration  and  1/s  in  (3.5.4)  with  the  neural  inte- 
grator system  described  in  (3.4.3),  it  appears  that 
T^^(sT'  + 1)  must  correspond  to  some  sort  of  equivalent 
see  with  an  effective  time  constant  of  t'  - 12  sec. 


c 

It  follows  that  there  must  exist  in  the  brain  a compen- 
sation process  for  the  See  represented  by  a transfer  function 
G (s)  such  that 


c 


T T 

c G (s)  = g - 

sT  + 1 c sT'  + 1 

c c 

from  which  we  obtain 


(3.5.5) 


T (sT  + 1) 

c'  ' T (sT-  + 1) 
c c 

where  T is  the  experimentally  measured  cupula 
c 

while  t'  is  the  effective  time  constant  for  an 
c 

see  determined  from  the  reflex  response. 


(3.5.6) 

time  constant 
equivalent 


It  is  unknown  where  and  how  the  activity  corresponding  to 

G (s)  as  given  l)y  hapi  ens  neurophysiologically . Our 

c 

use  of  = 12  instead  of  = 4 is  purely  pracjmatic. 

The  complete  compensatory  slow-phase  model  is  shown  in 
Figure  3.1,  in  which  the  reflex  gain  of  0.9  in  the  darkness 
is  inserted  between  the  output  of  the  effective  SCC  and  the 
input  of  the  neural  integrator.  All  animals  have  a reflex 
gain  less  than  one  in  the  dark.  The  reason  for  this  is  not 
known.  The  reflex  gain  is  brought  up  to  one  in  the  light 
presumably  because  of  the  contribution  from  the  smooth 
pursuit  or  optokinetic  systems. 
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Figure  3.1 

This  figure  shows  the  compensatory 
slow-phase  model.  The  SCC  is 
replaced  by  the  "effective"  SCC 
which  incorporates  a compensation 
process  in  the  brain  for  the  SCC 
to  yield  a proper  reflex  time 
constant.  The  reflex  gain  of 
0.9  in  the  darkness  is  inserted 
between  the  effective  SCC  and  the 
neural  integrator. 


Figure  3.1  Compensatory  Slow-Phase  Model 


CHAPTER  IV 


THE  MODELS  OF  SACCADES  AND  QUICK-PHASES 

4 . 1 General  description  of  saccades  and  quick-phases 

Light  enters  the  eye  through  the  lens  projecting  an 
image  on  a photosensitive  sheet  called  the  retina.  However, 
the  resolving  power  of  the  retina  is  not  uniform.  To  see 
an  object  with  best  visual  acuity  one  must  move  his  eye  to 
place  its  image  on  a special  area  in  the  retina  called  the 
fovea  which  has  a high  density  of  photoreceptors.  Target 
acquisition  in  this  manner  is  done  with  a fast  discrete  eye 
movement  called  a saccade.  The  extraocular  muscles  are 
among  the  fastest  in  the  body  and  saccades  are  their  fastest 
product.  During  a saccade  the  agonist  muscle  contracts 
almost  maximally,  and  the  antagonist  muscle  is  inhibited 
completely.  As  an  example  of  the  speeds  involved,  a 10 
degree  human  saccade  lasts  about  45  milliseconds  and  has  a 
peak  velocity  of  about  400  deg/sec. 

It  was  qualitatively  and  quantitatively  substantiated 
that  saccades  and  vestibular  quick-phases  are  essentially 
the  same  movement,  and  therefore,  probably  the  product  of 
the  same  neural  mechanism.  This  was  done  by  comparing  the 
shape  of  their  time  courses  and  the  relationship  between 
the  amplitude  and  duration  of  saccades  and  the  quick-phases 
of  vestibular  nystagmus  in  the  alert  monkey  and  showing  that 
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they  are  almost  the  same.  The  differences  were  so  small 
that  it  seems  unnecessary  to  hypothesize  separate  circuits 
for  the  two  types  of  rapid  eye  movements  (Ron,  Robinson, 
and  Skavensky,  1972) . It  will  bo  shown  subsequently  that 
rapid  eye  movements  are  made  by  cells  in  the  brainstem 
called  bursters.  Such  cells  burst  briefly  at  high  discharge 
rates  during  saccades  and  quick-phases.  No  one  has  ever 
reported  a cell  that  burst  for  one  type  of  movement  but  not 
the  other.  There  is  no  evidence  that  these  movements  are 
different  and  this  idea  is  generally  accepted. 

Historically,  the  term  saccade  denotes  a rapid  eye 
movement  to  a visual  target  while  the  term  quick-phase 
denotes  a rapid  eye  movement  made  by  the  VOR  even  in  the 
dark.  There  is  the  tacit  assumption  that  quite  different 
central  mechanisms  are  involved  although  the  final  movement 
machinery  is  the  same.  That  is,  no  distinction  exists  between 
saccadic  and  quick-phase  signals  at  the  level  of  the  neural 
pulse  generator  (the  bursters)  to  effect  rapid  eye  m<  vement, 
but  separate  signals,  one  for  the  saccade  and  the  other  for 
the  quick-phase  have  been  assumed  to  enter  this  pre-nuclear 
pulse  generating  system.  However,  vision  is  not  always 
necessary  for  saccades.  They  are  stimulated  by  auditory  and 
tactile  stimuli  as  well.  It  is  more  general  to  think  of  a 
saccade  as  an  eye  movement  to  a point  in  space  which  may 
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contain  something  of  interest.  If  that  is  the  case,  why 
not  admit  vestibular  signals  as  being  able  to  determine 
which  point  in  space  is  of  interest?  If  one  is  suddenly 
rotated,  say,  to  the  left,  what  lies  to  the  left  is  suddenly 
of  great  interest  because  that  is  where  one  is  going.  Is 
the  rapid  eye  movement  to  the  left,  which  in  fact  one  always 
makes,  a saccade  or  a quick-phase?  The  idea  will  be  proposed 
that  rapid  eye  movements  are  designed  to  bring  the  eye  to 
a point  in  space  of  interest  and  that  point  can  be  determined 
by  any  sense  modality  including  vestibular.  In  that  sense, 
there  may  not  be  any  distinction  between  a saccade  and  a 
quick-phase  even  in  the  higher  centers  which  control  the 
pre-nuclear  ijulse  generator. 

The  most  noticeable  event  during  a rapid  eye  movement 
is  the  burst  of  high  frequency  discharges  seen  to  occur  in 
the  motoneurons.  The  discharge  rates  of  the  oculomotor 
neurons  during  a QP  are  very  high  and  range  from  400  to  600 
spikes  per  sec.  This  drives  the  eye  at  velocities  up  to 
700  deg/sec  (in  man).  Cat's  QP  eye  velocity  is  much  less 
than  that  of  man.  For  our  cat,  it  was  about  120  deg/sec 
for  saccades  of  amplitudes  of  about  13.5  deg. 

Neurons  that  burst  with  a high  frequency  discharge 
time-locked  with  the  motoneuron  burst,  but  which  are  silent 
otherwise,  are  found  throughout  the  pons  (Sparks  and  Travis, 
1971).  However,  a great  concentration  of  so-called  bursters 
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are  seen  in  the  PPRF  (Cohen  and  Henn,  1972) . They  were 
recently  studied  in  some  detail  by  Keller  (1974).  Unilateral 
lesions  of  the  PPRF  abolish  ipsilateral  saccades  while  lesions 
elsewhere  do  not.  It  is  difficult  to  avoid  the  conclusion 
that  saccades  and  quicl<. -phases  are  generated  by  the  activity 
of  the  bursters  in  the  PPRF. 

The  neural  pulse  of  bursters  is  very  large  during  a 
saccade  and  zero  before  and  after  it,  just  li)ce  eye  velocity. 
Therefore,  it  is  assumed  that  the  eye  velocity  is  proportional 
to  the  pulse  height  (by  this  is  meant  the  firing  rate  in 
spikes/sec) . That  is,  the  pulse  is  an  eye  velocity  command. 

We  must  transform  the  neural  eye  velocity  information  (pulse 
height)  into  eye  position.  Consequently,  the  signal  must  be 
integrated  just  like  the  vestibular  signal.  So  we  use  the 
simplest  path  already  existing  in  the  slow-phase  model.  The 
consequent  arrangement  is  shown  in  Figure  4.1.  The  pathway 
between  the  VN  and  the  output  of  the  OMP,  that  is,  the  eye 
position  (as  shown  in  Figure  4.1)  is  defined  for  our  purpose 
as  the  extended  final  common  path  because  it  is  a common  path 
for  both  the  SP  and  QP  of  nystagmus.  The  neural  network 
that  produces  the  saccadic  pulse  is  termed  the  pulse  generator 
(PG) . A neural  pulse  generated  by  the  PG  enters  at  the  VN 
and  becomes  a pulse-step  (weighted  sum  of  pulse  and  step) 
at  the  OMN;  the  pulse  being  required  to  m.ove  the  eye  rapidly 
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Figure  4. I 

This  figure  shows  a conceptual  arrangement 
for  saccade  generation.  A neural  pulse 
generated  by  the  pulse  generator  is 
injected  at  the  VN  level  and  processed 
(integrated)  through  the  existing 
pathway  resulting  in  a saccade. 


against  orbital  viscous  forces  and  the  step  to  hold  the  new 
eye  position  against  elastic  restoring  forces.  The  relative 
weight  between  the  pulse  and  step  at  the  OMN  is  0.2  and  1.0 
respectively,  so  that  their  sum  gives  the  step  (after  it  is 
passed  through  the  OMP)  which  we  call  a saccade.  This  is 
assured  because  the  equivalent  transfer  function  between 
the  VN  and  eye  position  is  1/s  or  a pure  integrator  as 
explained  in  the  last  chapter  on  the  slow-phase  model. 

This  scheme  is  generally  accepted  among  oculomotor 
neurophysiologists.  The  question  now  is:  How  is  tne  pulse 

generated?  How  is  the  pulse  size  made  correct  for  the 
desired  saccade  size?  This  will  be  discussed  in  the 
following  sections. 

4 . 2 Processing  of  information 

The  problem  confronting  us  is  how  the  burst  size  is 
determined  for  a given  required  saccade  size.  There  exist 
two  fundamentally  different  views  as  to  what  information 
the  PG  uses  to  generate  a saccade  (Robinson,  1973;  Zee 
et  al.,  1976).  The  first  view  is  that  the  saccadic  commands 
coming  to  the  PG  are  coded  in  the  distance  by  which  the  eye 
must  be  moved.  The  second  view  is  that  the  saccadic  commands 
are  coded  in  the  position  in  the  orbit  to  which  the  eye  must 
be  brought. 
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The  former  view  is  illustrated  in  Figure  4.2.  The  PG 
does  not  have  to  know  the  initial  and  the  final  coordinates 
of  the  eye  but  rather  it  has  to  recognize  only  the  difference 
between  two  points  on  the  retina  and  then  drives  the  eye 
with  a high  velocity  by  a pulse  the  area  of  which  (pulse 
height  times  pulse  duration)  is  equal  to  the  retinal  differ- 
ence. This  view  is  called  a retinotopic  organization 
because  calculations  are  presumed  to  be  made  only  in  a 
retinal  coordinate  system. 

The  latter  view  is  shown  in  Figure  4.3.  Based  on  the 
coordinates  of  the  eye  position  at  the  beginning  of  a saccade 
and  the  coordinate  of  eye  position  that  is  desired,  the  pulse 
generator  (a  negative  feedback  system)  drives  the  eye  with 
a high  velocity  until  it  reaches  its  goal.  This  is  called 
a spatially  organized  system  because  the  main  signals  are 
specified  in  a head-coordinate  reference  frame.  The  former 
hypothesis  superficially  appears  to  be  simplest  but  only 
because  it  glosses  over  many  of  the  details.  It  is  inherent 
so  far  in  most  control  systems  models  of  saccadic  eye  move- 
ments (Young  and  Stark,  1963;  Robinson,  1973) . 

However,  we  know  that  the  brain  makes  saccades  when 
there  is  no  retinal  error.  For  example,  we  make  saccades 
in  the  dark  to  tactile  and  auditory  stimuli  without  any 
retinal  error.  As  another  excimple,  if  a point  source  of 
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^ Figure  4.2 
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This  figure  shows  the  outline  of  the 
control  scheme  for  retinotopic 
organization,  in  which  the  saccadic 
commands  (neural  pulses)  are  coded 
in  distance  (the  difference  between 
two  points  on  the  retina)  by  which 
the  eye  must  be  moved. 


I 
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Figure  4.3 

This  figure  outlines  the  scheme  for 
spatial  organization,  in  which  the 
saccadic  commands  are  coded  in  the 
position  in  the  orbit  to  which  the 
eye  must  be  brought. 
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light  is  briefly  flashed  on  the  fovea  of  an  eye  moving 


during  a saccade,  the  subject,  who  has  been  told  to  track 
the  target,  makes  another  saccade  back  to  the  light  source 
from  the  end  point  of  the  first  saccade  even  though  the 
retinal  error  was  zero  at  the  time  of  the  flash  (Hallet  and 
Lightstone,  1976).  This  is  possible  only  if  the  subject 
could  know  the  position  of  the  target  at  the  time  of  the 
flash  in  space  rather  than  on  the  retina. 

In  retinotopic  organization,  the  saccade  is  assumed  to 
be  ballistic  or  preprogrammed.  That  is,  the  product  of  the 
amplitude  and  duration  of  the  pulse  (which  is  equal  to  the 
magnitude  of  the  saccade)  is  determined  before  it  begins 
and  thus  cannot  be  modified  in  midflight  once  it  has  started. 
However,  it  was  discovered  that  patients  who  made  abnormally 
slow  saccades  could  interrupt  a saccade  in  midflight  to 
acquire  a new  target  (Zee  et  al.,  1976).  This  cast  consid- 
erable doubt  on  the  idea  that  saccades  are  preprogrammed. 

For  these  reasons,  we  have  adopted  in  our  model  the 
spatial  organization  of  saccades  in  which  the  brain  computes 
the  position  of  the  target  (or  the  position  where  the  eye 
wants  to  go)  in  head  coordinates  and  then  simply  drives  the 
eye  rapidly  towards  the  target  (or  its  goal)  until  it  gets 
there. 


The  brain  could  monitor  the  eye  position  in  the  head 
either  by  using  an  "efference  copy"  of  the  eye  position 
derived  from  internally  monitored  commands  to  the  OMN  or  by 
using  direct  afferent  information  from  extraocular  muscle 
spindles  and  other  sensory  receptors.  Most  studies  (e.g., 
Skavenski  et  al.,  1972)  suggest  efference  copy  is  used  to 
compute  the  location  of  seen  objects  in  head  coordinates. 
Thus,  in  our  model,  the  eye  position  in  head  coordinates  is 
obtained  from  the  efference  copy  (0'  in  Figure  4.1)  for 
the  purpose  of  saccadic  pulse  generation. 

4 . 3 Control  scheme  for  the  pulse  generator 

We  want  a mode  of  control  which  drives  the  eye  as 
rapidly  as  possible  to  the  point  where  it  wants  to  go  from 
the  point  where  it  is  now  based  on  the  tnowledge  of  the 
coordinates  of  tx)th  points  in  the  head-coordinate  frame. 

A classic  mode  of  a control  system  which  performs  this 
kind  of  function  is  the  bang-bang  control  system.  We  use 
a modified  version  of  it  to  match  the  behavior  of  our  cat  as 
shown  in  Figure  4.4.  It  consists  of  a simple  negative 
feedback  system  designed  to  drive  error  signal  to  a null. 

Its  forward  path  contains  a high-gain  saturating  amplifier. 
The  reference  input  C(t)  to  the  system  is  the  target  position 
or  the  point  in  space  where  the  eye  wants  to  go.  The  out- 
put of  the  amplifier  is  Integrated  by  the  NI  whose  output 


Figure  4.4 

This  is  a bang-bang  control  system 
with  high-gain  amplifier  modified 
for  our  purpose  from  the  standard 
form.  The  reference  input  C(t)  is 
the  target  position.  The  control 
system  drives  the  efference  copy 
6^(t)  of  the  eye  position  to  the 
target  position  C(t)  by  means  of 
feedback  control. 


» 
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becomes  the  cfference  copy  of  the  eye  position.  It  follows 
that  the  output  of  the  amplifier  must  be  proportional  to 
the  eye  velocity.  For  this  reason,  the  output  of  the  ampli- 
fier is  sometimes  referred  to  as  the  velocity  command.  The 
efference  copy  6^(t)  is  fed  back  to  the  input  to  produce 
the  error  signal  A6(t)  which  is  the  difference  between  the 
target  position  C(t)  and  the  eye  position  copy  0^(t).  When 
the  eye  position  catches  up  with  the  target  position,  the 
error  signal  is  reduced  to  zero  and,  consequently,  the 
velocity  command  vanishes. 

The  high-gain  amplifier  in  Figure  4.4  has  infinite 
gain  at  A0  = 0 so  the  system  as  shown  would  be  unstable. 

At  the  end  of  one  quick-phase  it  could  jump  to  the  opposite 
limb  of  P(A0)  and  so  start  oscillating.  In  order  to  stop 
this,  we  have  to  construct  some  associated  gating  circuitry. 
This  will  be  done  in  the  next  section. 

Also  the  shape  of  the  amplifier  curve  P(A0)  must  be 
chosen  in  sach  a way  that  amplitude  and  duration  of  quick- 
phases  generated  by  the  model  are  compatible  with  those  of 
the  cat.  The  method  of  determining  this  curve  is  also  given 

I 

in  the  next  section. 
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4 . 4 The  description  of  the  quick-phase  generation  mechanism 


The  proposed  mechanism  of  quick -phase  general- ion  is 

shown  in  more  detail  in  Figure  4.5.  The  trigger  signal  which 

initiates  a quick -phase  is  called  the  "WHEN  signal."  The 

trigger  signal  is  assumed  to  be  discrete.  The  continuous 

signal  which  represents  the  target  position  C(t)  where  the 

eye  should  be  driven  is  called  the  "WHERE  signal." 

Referring  to  Figure  4.5,  a discrete  WHEN  signal  actuates 

the  OR-gate,  thus  closing  the  switch  s which  starts  the 

2 

quick-phase.  Once  a quick-phase  is  initiated,  a branch  of 

the  output  P(t)  of  the  amplifier  is  channeled  back  to  the 

OR-gate  so  that,  even  after  the  WHEN  signal  (being  a discrete 

trigger  signal  by  assumption)  disappears,  the  OR-gate  is 

kept  activated  as  long  as  the  output  of  the  amplifier  persists 

at  non-zero  values.  If,  during  a quick-pl.ase , the  efference 

copy  0^(t)  matches  the  WHERE  signal  C(t),  thus  making  A0  = 0, 

the  output  of  the  amplifier  first  becomes  zero  which  in  turn 

opens  the  switch  s by  deactivating  the  OR-gate.  This  makes 
2 

the  QP  timing  autonomous.  Once  started,  it  runs  its  course 
to  conpletion,  independent  of  the  WHEN  signal. 

The  expanded  view  of  the  modified  version  of  the  high- 
gain  amplifier  as  used  in  Figure  4.5  is  shown  in  Figure  4.6. 

It  has  a built-in  arrangement  to  prevent  the  oscillation 
which  would  otherwise  be  induced  because  in  digital  simula- 
tion the  last  increment  of  input  reaching  the  origin  from 


Figure  4.5 

This  figure  shows  the  outline  of  the 

quick-phase  generation  mechanism.  A 

discrete  WHEN  signal  closes  the  switch 

s through  the  OR-gate.  Once  s is 
2 2 

closed,  it  remains  closed  by  the  output 
P(t)  of  the  amplifier  through  the  OR- 
> gate  until  the  efference  copy  6'(t)  of 

the  eye  position  matches  the  WHERE 
signal  C(t)  which  is  the  target  position. 
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Figure  4.6 

This  is  the  expanded  view  of  the  high- 


gain  amplifier  used  in  Figure  4.5.  For 
positive  A0  [where  A0  = C(t)  - 0^(t)], 
the  upper  half  of  the  figure  denoted 
by  P'*'(A0)  is  used.  For  negative  A0, 
the  lower  half  of  the  figure  denoted 
by  P"(A0)  is  used.  This  arrangement 
pertains  only  to  our  computer 
simulation  in  order  to  eliminate  the 
oscillation  at  A0  = 0. 
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THE  ERROR  ANGLE  FROM  THE  EYE  POSITION 
TO  THE  WHERE  CURVE 


Figure  4.6  The  High  Gain  An^lifier  (used  in 
Figure  4.5) 


one  side  of  the  input  axis  generally  overshoots  slightly 
into  the  region  with  the  opposite  sign. 

The  high-gain  amplifier  in  Figure  4.6  is  broken  up  into 
two  curves:  P (A9)  and  P (A6) . If  A0  > 0 when  a quick-phase 

is  triggered,  the  model  picks  the  P^(A0)  curve.  If  A0  < 0 
when  a quick-phase  is  triggered,  the  model  picks  the  P (A0) 
curve.  This  arrangement  pertains  only  to  our  computer 
simulation.  Just  how  the  brain  might  handle  this  problem 
is  unknown. 

The  quick-phase  generation  mechanism  works  in  the 
following  way.  Suppose  the  head  is  turning  to  the  right. 

The  error  angle  from  the  eye  position  to  the  place  where 
it  wants  to  go  by  a quick-phase  would  be  positive,  or  in 
this  case  A0  = C(t)  - 0^(t)  > 0.  So,  the  upper  half  of 
Figure  4.6  applies  for  quick-phase  pulse  generation  caused 
by  the  right  turn  head  motion.  Note  that,  in  this  case, 
once  a quick-phase  is  initiated  by  a positive  A0,  the  quick- 
phase  will  be  terminated  even  if  A0  overshoots  to  the 
negative  region  after  crossing  the  origin  from  the  right 
because  the  output  P(A0)  of  the  amplifier  is  zero  for  A0  < 0 
in  the  upper  half  of  Figure  4.6,  thus  preventing  the  oscil- 
lations. For  the  left  turn  head  motion,  the  error  angle 
is  negative  or  A0  < 0.  So,  for  the  quick-phase  pulse  gen- 
eration, the  lower  left  half  of  Figure  4.6  is  used  in  an 
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exactly  opposite  way  from  the  case  for  that  caused  by  the 
right  turn  head  motion.  The  upper  and  lower  halves  correspond 
to  two  unidirectional  pulse  generators  in  each  half  of 
the  pons. 

The  logic  which  determines  the  WHEN  signals  and  the 
WHERE  signals,  that  is,  the  timings  and  the  anplitudes  of 
quick-phase,  is  the  central  problem  of  this  thesis  cuid  will 
be  discussed  in  detail  in  the  next  chapter. 

The  next  step  is  to  determine  the  input-output  relation- 
ship P(A0)  of  the  amplifier  in  order  to  generate  saccades  with 
appropriate  velocities  for  the  cat.  In  general,  P(A6) 
increases  with  AQ  to  reflect  the  physiological  observation 
that  large  saccades  have  a higher  velocity  than  small  saccades. 
Let  the  input  to  the  amplifier  be  denoted  by  A6(t).  Note  that 
A6(t)  = C(t)  - e"(t)  (4.4.1) 

For  the  steady  state  response  of  the  step  input  head  velocity, 
C(t)  changes  only  slowly  compared  to  the  rapidity  of  quick- 
phases.  Therefore,  it  follows  that 

A0(t)  = -6"(t)  (4.4.2) 

Since  6^(t)  is  the  approximate  input  to  the  neural  integrator, 
it  must  be  approximately  equal  to  the  output  P(t)  of  the 
amplifier  because  during  the  quick-phase  the  only  input  to 
the  neural  integrator  is  the  output  from  the  amplifier.  It 
follows  that 
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p(A0(t))  = e"(t)  = -A0(t) 


(4.4.3) 


The  actual  curve  P(A0)  is  smooth;  that  is,  it  has  a contin- 
uous derivative  (see  Zee  et  al.,  1976:  Keller,  1974)  but 

a straight  line  segment  approximation  is  adequate.  In  fact, 
because  this  model  is  only  concerned  with  QP's  as  events 
and  is  not  concerned  with  great  accuracy  concerning  intra- 
saccadic  dynamics,  we  shall  assume  a simple  linear  relation- 
ship between  the  input  and  the  output  of  the  amplifier  (shown 
in  Figure  4.7  for  A0  > 0)  in  the  form  of 

P+(A0)  = aA0  + b for  A0  > 0 (4.4.4) 

and 


P"(A0)  = aA0  - b for  A0  < 0 (4.4.5) 

Equation  (4.4.4),  using  (4.4.3),  may  be  expressed  as 

A0(t)  = -aA0(t)  - b (4.4.6) 

The  solution  of  (4.4.6)  is 


A0(t)  = 


e~^'^A0(O) 


•/: 


-a(t-T) 


(-b)dx 


(4.4.7) 


or 

A0(t)  = e‘^’^A0(O)  - - (1  - e"^'^)  = e~^^(A0(C)  + - 

a \ a / a 

(4.4.8) 

We  want  to  determine  the  value  of  a )cnowing  A0(O)  at  the 
beginning  t = 0 of  the  quick-phase,  A0(t^)  at  the  end  t = t^ 
of  the  quick-phase  and  the  intercept  point  b of  p'*^(A0)  with 
the  output  axis  in  Figure  4.7.  The  lowest  quick-phase 
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Figure  4.7 

This  is  the  model  used  to  determine 
the  equation  for  the  bang-bang 
cimplifier  (for  A0  > 0)  . The  figure 
describes  the  case  for  right  turn 
head  motion  in  which  A0  > 0.  The 
problem  is  to  determine  the  slope, 
and  the  intercept  point  b for  the 
equation  P^(A0)  = aA0  + b which 
describes  the  input-output 
relationship  of  the  bang-bang 


velocity  6^  from  our  cat  is  40.5  deg/sec  corresponding  to  the 
slow-phase  velocity  of  3.6  deg/sec.  Thus,  we  chose  the  value 
of  40  for  b.  The  largest  mean  quick-phase  cunplitude  in  this 
study  (for  the  highest  mean  input  head  velocity  of  60.6  deg/sec) 
is  13.5  deg  with  the  mean  quick-phase  duration  of  0.115  sec. 

The  quick-phase  duration  of  0.115  sec  for  a 13.5  deg  rapid  eye 
movement  also  coincides  with  the  results  obtained  on  the 
cat's  saccadic  amplitude  - duration  relationship  (an  unpub- 
lished observation  by  D.  A.  Robinson) . 

Substituting  A0(O)  = 13.5,  A6(t)  = A6(t^)  =0,  ^ ~ 

0.115  and  b = 40  in  (4.4.8),  we  have 


0 = e-"^°-'^^Nl3.5t^) 

a a 


(4.4.9) 


Solving  (4.4.10)  for  a (ay  trial  and  error)  we  obtain 
a = 16.2.  It  follows  from  (4.4.4)  and  (4.4.5)  that 

P'^(A0)  = 16.2  A0  + 40  (4.4.10) 

and  similarly  for  P (A6) . 

This  model  will  now  automatically  generate  cat  saccades 
which  have  roughly  the  correct  amplitude  - duration  cuid 


amplitude  - velocity  relationship  for  saccades  of  all  sizes. 


CHAPTER  V 


THE  INTEGRATED  MODEL  OF  THE  SLOW-PHASE  AND  THE  QUICK-PHASE 

5 . 1 The  outline  of  the  integrated  model 

The  first  attempt  to  integrate  the  slow-phase  and  the 
quick-phase  models  is  made  simply  by  summing  the  output  ())(t) 
of  the  see  and  the  output  P(t)  of  the  NPG  as  shown  concep- 
tually in  Figure  5.1.  This  idea  conforms  with  the  observed 
modulation  of  saccade  velocity  during  head  movement  in  a 
study  of  monkey  head-eye  coordination  (Morasso,  et  al.,  1973). 
However,  in  those  studies  it  is  the  case  that  for  monkey 
P{t)  » (()(t)  . 

But  cat  behaves  differently.  For  our  cat,  which  was 
typical,  the  quick-phase  or  saccade  eye  velocity  for  ampli- 
tudes of  13  deg  is  about  120  deg/sec  compared  with  about 
850  deg/sec  for  monkey  (Fuchs  and  Luschei,  1970) . The 
arrangement  of  Figure  5.1  would  work  so  long  as  P(t)  is 
always  much  greater  than  <()(t).  However,  if  P (t)  is  near 
or  less  than  <t)(t)  there  is  trouble.  Suppose,  for  instance, 

P (t)  = 60  deg/sec  and  <t)(t)  = 80  deg/sec.  The  arrangement 
would  not  function  properly  Ijecause  the  eye  would  still  be 
moving  in  the  direction  of  the  slow-phase  instead  of  the 
quick-phase  while  the  neural  pulse  was  being  generated  by 
the  NPG  for  a quick-phase.  Thus,  a proper  quick-phase  could 
not  occur,  contrary  to  fact.  This  raises  a question  as  to 
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Figure  5.1 


This  figure  depicts  the  hypothesis 
that  SP  and  QP  commands  linearly 
summate.  This  hypothesis  is  rejected 


because  it  contradicts  existing 


evidence  and  leads  to  contradictions 


when  the  output  of  the  neural  pulse 


generator  is  near  or  less  than  the 


output  of  the  see. 


I 


whether  the  output  from  the  SCC  should  bo  suppressed  from 
entering  the  extended  final  common  path  (Figure  5.1) 
during  quick-phases. 

For  example  as  supporting  evidence,  many  cells  in  the 
vestibular  nucleus  which  fire  with  a rate  proportional  to 
head  velocity  (the  4>(t)  signal)  pause  during  all  quick-phases 
(Fuchs  and  Kimm,  1975).  This  is  neural  evidence  that  the 
vestibular  signal  is  suppressed  during  '^uick-phases . Further- 
more, a recent  study  (Jurgens  and  Becker  1975)  indicates  that 
the  hypothesis  of  a linear  addition  of  saccades  and  pursuit 
eye  movements  does  not  work.  That  is,  the  smooth  pursuit 
signal  did  not  appear  to  continue  during  a saccade  and 
therefore,  did  not  appear  to  be  simply  added  or  subtracted 
from  the  saccadic  signal.  Using  these  findings  as  supporting 
evidence  for  suppressing  the  slow-phase  signal  during  quick- 
phases,  our  model  implements  this  by  introducing  a switch  s 

1 

at  the  output  of  the  SCC  as  shown  in  Figure  5.2.  s is 

1 

closed  during  s low-phases  and  open  during  quick-phases.  With 
this  arrangement  the  input  to  the  extended  final  common  path 
is  the  signal  from  the  SCC  only  during  slow-phases  and  the 
signal  from  the  NI'G  only  during  quick-phases. 

Incorporating  Figure  5.2  with  Figure  4.5,  we  obtain 
Figure  5.3  which  reriesents  the  outline  of  the  integrated 
model.  The  role  of  the  AND-gate  is  to  assure  a reaction 


Figure  5.2 


This  figure  depicts  the  hypothesis 
that  the  canal  signal  is  suppressed 
by  quick-phases.  The  switch  s^  is 
closed  for  slow-phases  and  open  for 
quick-phases  thus  making  the  output 
from  the  neural  pulse  generator  the 
only  input  to  the  extended  final 
Cf>mmon  path  during  quick-phases. 
This  arrangement  is  adopted  in  this 
study . 
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Sy.  CLOSED  FOR  SLOW-PHASES 
OPEN  FOR  QUICK-PHASES 


Figure  5.2  The  Hypothesis  that  the  Canal 
Suppressed  by  Quick-Phases, 
is  used  in  this  study) 


TO  THE  EXTENDED 
FINAL  COMMON  PATH 


Signal  is 
(This  scheme 
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Figure  5.3 

This  figure  shows  the  outline  of  the 

integrated  model  which  is  achieved  by 

incorporating  Figure  5.2  with  Figure  4.5 

Since  s is  open  and  s is  closed,  the 
1 2 

system  is  undergoing  a quick-phase . 

One  of  the  inputs  to  the  AND-gate 
labeled  as  "50  milliseconds  refractory 
period"  assures  that  a period  of 

m 

50  milliseconds  must  pass  before  another 
quick-phase  can  be  made. 
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WHERE  SIGNAL 

(THOUGHT  OF  AS  CONTINUOUS  SIGNAL) 


Figure  5.3  Outline  of  the  Integrated  Model 

(s  and  s in  QP  Mode) 

1 2 
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time  of  50  milliseconds  (called  the  refractory  period)  for 
the  brain  before  making  another  quick -phase.  It  is  a well- 
known  fact  that  after  one  saccade,  a period  of  around  200 
milliseconds  must  pass  before  another  saccade  can  be  made. 
This  rule  also  seems  to  apply  to  QP's  with  50  milliseconds 
instead  of  200  milliseconds.  This  choice  of  50  milliseconds 
will  be  justified  in  Section  7.4.  Thus,  the  AND-gate  dis- 
regards the  WHEN  signal  unless  50  milliseconds  has  elapsed 
from  the  end  of  the  previous  quick-phase. 

Next,  we  consider  the  WHEN  signal  and  the  WHERE  signal. 
As  mentioned  before,  the  WHEN  signals  are  discrete  trigger 
signals  which  initiate  quick-phases,  while  the  WHERE  signal 
is  thought  of  as  a continuous  signal  which  represents  the 
target  position  of  where  the  eye  wants  to  go.  The  WHERE 
signal  is  the  output  of  a system  to  be  determined  in  detail 
in  a later  section.  For  the  moment,  we  want  to  see  how  these 
signals  shape  the  pattern  of  nystagmus. 

For  rotatory  nystagmus  the  WHERE  curve  may  be  thought 
of  conceptually  as  a smooth  curve  which  contains  the  end- 
points of  all  quick-phases  as  shown  schematically  in  Fig- 
ure 5.4A.  Since  the  WHEN  signal  terminates  a slow-phase  and 
initiates  a quick-phase,  for  each  WHEN  signal  there  is  a 
corresponding  eye  position  which  is  the  end-point  of  a slow- 
phase.  It  is  useful  to  imagine  a WHEN  curve  also  as  a 
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This  figure  shows  a typical  VOR 
response  from  our  cat  for  step 
head  input  velocity  of  about 
60  dog/sec.  The  WHEN  signals 
and  the  WHERE  signals  are  thought 
of  as  continuous  curves  that  are 
sampled  by  the  occurrence  of  a QP. 
The  WHERE  curve  C(t)  is  divided 
into  two  parts:  the  deterministic 

part  C(t),  equal  to  the  mean  value 
of  C(t),  and  noisy  part  n(t)  with 
zero  mean.  Similar  considerations 
apply  for  the  WHEN  curve. 


WHERE  CURVE 


Figure  5.4A  A Typical  VOR  Response  from  our  Cat  for  H ==  60  deg/sec 


continuous  signal.  Whenever  the  eye  reaches  it,  a quick- 
phase  is  triggered.  Consequently,  it  contains  the  end-points 
of  all  the  slow-phases  as  shown  in  Figure  5.4A.  The  WHEN 
signal  is  the  output  of  another  system  to  be  determined  in 
a later  section.  The  typical  appearances  of  the  WHERE  curve 
and  the  WHEN  curve  are  similar  to  those  shown  in  Figure  5.4A, 
which  is  a typical  VOR  response  from  our  cat  for  the  step 
head  input  velocity  of  about  60  deg/sec,  the  highest  value 
we  have  used.  Actual  responses  from  our  cat  for  this  input 
velocity  are  shown  in  Figures  5.4B,  5.4C,  5.4D,  5.4E,  and 
5.4F. 

The  WHERE  and  WHEN  curves  are  internal  signals  which 
cannot  be  observed  as  continuous  signals.  However,  each 
time  a quick-phase  occurs,  tnese  signals  are  made  explicit. 
Thus,  experimentally  wc  observe  them  as  sampled  signals 
indicated  by  the  dots  in  Figure  5.4A. 

These  signals  are  obviously  noisy.  When  head  velocity 
is  constant,  however,  nystagmus  is  more  or  less  in  a steady 
state.  In  this  situation  we  can  divide  these  signals  into 
their  mean  value  (deterministic  part)  and  a noisy  part  which 
has  zero  mean  value.  Thus,  the  curve  C(t)  can  be  divided  into 
C(t)  = C(t)  f n(t)  (5.1.1) 

where  C(t)  is  the  deterministic  WHERE  curve  and  n(t)  is  the 
noise,  (Figure  5.4A).  Similar  considerations  apply  for  the 


WHEN  curve. 


Figures:  5.4B,  5.4C,  5.4D,  5.4E  and  5.4F 
These  figures  show  more  VOR  responses  from 
our  cat  for  step  head  input  velocity  of 
about  60  deg/sec.  Figures  5.4B  and  5.4E 
show  the  responses  from  left  turn  head 
motion,  while  Figures  5.4C,  5.4D,  and  5.4F 
show  the  responses  from  right  turn  head 
motion.  Note  that  variations  in  the  QP 
end-points  (WHERE  the  eye  is  going)  are 
considerably  less  than  the  variations 
in  the  SP  end-points  (WHEN  the  QP  is 
initiated) . 
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Figure  5.4B  VOR  Response  from  Our  Cat  for  Step  Head  Input  Velocity  of  About 
60  deg/sec  (for  Left  Turn  Head  Motion) 


Figure  5.4C  VOR  Response  from  Our  Cat  for  Step  Head  Input  Velocity  of  About 
60  deg/sec  (for  Right  Turn  Head  Motion) 


NAVAL  SURFACE  NEA^ONS  CENTER  OAHLSREN  LAB  VA 
NATMEMATICAL  model  of  the  VESTXBULOOCULAR  reflex. (U) 
JUN  77  K S CHUN 


Figure  5.4F  VOR  Response  from  Our  Cat  for  Step  Head  Input  Velocity  of  About 
60  deg/sec  (for  Right  Turn  Head  Motion) 


5 . 2 The  contours  of  the  deterministic  WHEN  curve  and  VfflERE 

curve 

Referring  to  Figure  5.4  in  the  last  section,  the  most 
striking  aspect  of  the  nystagmus  pattern  is  that  the  mean 
eye  position  of  our  cat  at  the  conclusion  of  quick-phases 
in  the  steady  state  is  shifted  well  into  the  direction  of 
the  head  rotation.  Earlier  experiments  with  human  subjects 
demonstrated  that  when  a subject  in  the  dark  is  rotated  (say 
to  the  right) , quick-phases  drive  the  eye  in  the  direction 
of  turning  (to  the  right  in  this  case)  (Melvill  Jones,  1964). 
This  observation  supports  the  idea  mentioned  in  section  4.1 
that  the  purpose  of  a quick-phase  is  not  simply  to  reset  the 
eye  after  a slow-phase  but  to  allow  a subject  to  see  where 
he  is  going  by  turning  the  eye  to  look  in  that  direction. 
Thus,  a quick-phase  is  an  active  thing  seeking  new  visual 
information  which  lies,  of  course,  in  the  direction  in  which 
one  is  turning. 

Our  cat's  data  show  that  for  various  step  input  head 
velocities  the  eye  makes  quick-phases  in  the  direction  of 
the  head  motion  by  about  12  degrees.  Since  the  limit  of  the 
cat's  eye  movement  range  is  only  20  degrees,  it  is  clear 
that  even  at  moderate  head  velocities,  the  cat  is  looking 
ahead  almost  as  far  as  it  can.  From  this,  one  can  make 
several  observations  about  the  WHERE  curve  C(t).  First,  it 
obviously  must  be  derived  from  the  canal  signal  <))(t)  since 
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in  the  dark  that  is  the  only  siqnal  which  informs  the  brain 
that  rotation  is  taking  place.  Second,  it  should  promptly 
call  for  an  eye  position  in  the  direction  of  turning  cuid, 
finally,  since  the  range  of  eye  position  is  limited,  one 
might  expect  an  early  saturation  effect  by  which  C(t)  stops 
increasing  as  H continues  to  increase. 

Thus,  during  each  nystagmus  beat,  the  WHERE  curve  is 
the  "center  of  interest"  or  the  place  where  the  eye  is 
purposefully  driven  by  quick-phases.  Once  there,  the  slow- 
phase  maintains  the  visual  axis  at  that  point  in  space  in 
search  of  new  information.  This  requires  a certain  minimum 
amount  of  time  for  visual  processing  (at  least  50  milliseconds) 
but  beyond  that  the  gap  between  the  position  of  the  eye  and 
the  "center  of  interest"  C(t)  is  widening  and  eventually, 
probably  by  some  sort  of  threshold  phenomenon,  the  eye  is 
released  and  is  then  quickly  returned  to  C(t)  in  search  of 
now  visual  information. 

As  we  can  observe  in  Figure  5.4A,  the  steady  state 
deterministic  WHEN  curve  from  our  cat  seems  to  be  shifted 
by  a constant  amount  Q in  the  direction  of  the  slow-phase 
eye  motion  from  the  steady  state  deterministic  WHERE  curve 
in  such  a way  that  the  deterministic  WHEN  curve  appears  to 
approximately  coincide  with  the  zero  degree  line.  This 
suggests  that  the  difference  between  0(t)  and  C(t)  is  the 
critical  parameter  which  triggers  a quick-phase. 
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Another  observation  supports  the  idea  that  the  quick- 
phase  is  an  active  goal-seeking  process.  Even  casual 
inspection  of  Figure  5.4  shows  that  the  noise  on  the  WHERE 
curve  is  about  half  of  that  on  the  WHEN  curve  (analysis  of 
the  data  is  given  in  Chapter  6) . This  suggests  from  the 
systems  control  viewpoint  that  the  WHERE  curve  is  considered 
to  be  more  important  than  the  WHEN  curve  and  the  system 
might  be  organized  purposely  to  effect  more  control  (less 
noise)  over  the  former  than  the  latter. 

5 . 3 The  model  of  the  deterministic  '.."HERE  curve 

A set  of  the  VOR  responses  from  our  cat  for  step  input 
head  velocity  of  about  60  deg/sec  is  shown  in  Figures  5.4B, 
5.4c,  5.4D,  5.4E  eind  5.4F.  The  head  velocity  of  60  deg/sec 
is  selected  for  the  model  because  the  higher  the  input 
velocity  the  more  frequent  and  uniformly  spaced  are  the  quick- 
phases.  The  step  responses  of  the  deterministic  WHERE 
curves  are  drawn  for  best  fits  based  on  visual  inspection. 
Unfortunately,  the  sampling  rate  is  still  not  very  high  and 
one  has  difficulty  in  determining  any  details  in  the  shape 
of  any  early  transients.  However,  by  running  each  test 
several  times  one  can  observe  ensembles  of  responses  and 
meike  a better  estimate,  but  still  the  fitting  process  is  all 
based  on  visual  inspection.  Idealized  deterministic  curves 
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are  shown  in  Figure  5.5  for  a right  turn  head  motion  and  in 
Figure  5.6  for  a left  turn  head  motion.  The  step  response 
of  the  deterministic  WHERE  curve  C(t)  seems  to  consist  of 
two  parts.  The  first  part  is  the  portion  which  appears  to 
be  a response  of  a first-order  lag  system  (time  constant, 

T ) to  a step  input  with  steady  state  response  magnitude  f . 

L 0 

The  second  part  is  the  portion  represented  by  a step  function 
with  magnitude  g^. 

As  already  mentioned,  we  must  use  4>(t)  as  the  input 
signal  from  which  C(t)  is  generated.  In  fact,  since  the 
compensatory  slow-phase  model  is  a second  order  linear  system, 
the  velocity  signal  and  the  position  signal  are  the  state 
variables  of  the  system.  We  have  already  used  the  latter 
in  the  bang-bang  controller.  Consequently,  the  use  of  (})(t), 
the  velocity  signal,  as  the  other  input  to  the  mechcinism  of 
the  neural  pulse  generation  seems  appropriate  simply  from  a 
system  engineering  viewpoint. 

For  each  nominal  intended  head  input  velocity  H^,  the 
mean  values  of  actual  input  head  velocity  H,  the  slow-phase 
eye  velocity  6 , and  the  time  constcint  T of  the  first-order 

S L 

lag  system  as  well  as  f^  and  g^  are  determined  from  our  cat's 
analogue  data.  The  results  are  shown  in  Table  5.1.  H is  the 
actual  head  input  velocity  recorded  in  the  VOR  analogue  data. 

H is  not  identical  to  because  the  apparatus  used  for  our 
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Figures  5.5  and  5.6 

These  figures  show  deterministic  WHERE 
and  WHEN  curves  for  step  input  head 
velocity  for  the  right  turn  head  motion 
(Figure  5.5)  and  for  the  left  turn 
head  motion  (Figure  5.6). 
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Figure  5.5  Deterministic  WHERE  Curve,  C(t)  and  WHEN  Curve 
(For  Right  Turn  Head  Motion) 


experiment  is  a hand-operated  instrument  without  a feedback 
arrangement.  0^  is  the  mean  slow-phase  eye  velocity  in  deg/ 
sec  obtained  by  dividing  the  mean  of  all  slow-phase  magni- 
tudes of  several  runs  for  a given  H by  the  mean  of  all  slow- 
phase  durations  of  the  same  several  runs  for  the  same  H.  T 

L 

is  the  time  constant  of  the  first-order  system  determined  by 
visual  inspection  based  on  the  fact  that,  in  the  step  response 
of  the  first-order  lag  system,  the  initial  slope  line  crosses 
the  steady  state  response  position  one  time  constant  later 
(as  shown  in  Figure  5.5). 


Table  5.1 


0 

H 

0 

s 

gain 

^0 

^0 

^0  ■ 

deg/sec 

deg/sec 

deg/sec 

(6  /H) 
s 

deg 

deg 

deg 

sec 

10 

10.0 

8.9 

0.890 

4.2 

9.5 

13.7 

0.9 

20 

23.4 

21.0 

0.897 

4.2 

10.3 

14.5 

0.8 

40 

42.0 

30.8 

0.733 

3.6 

11.9 

15.5 

0.4 

60 

60.6 

46.5 

0.767 

2.3 

13.0 

15.3 

0.3 

All 

values 

are  mean 

values  of 

about 

. seven 

runs 

for  each  nominal  head  velocity  with  standard 
deviations  ranging  from  30%  to  50%  of  their 
values.  The  symbols  used  in  this  table  are 
defined  and  explained  elsewhere  in  this 
section  (5.3). 
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There  were  about  seven  runs  for  each  nominal  head 


velocity,  with  each  run  lasting  about  6 seconds,  which  is 

about  one  half  of  the  effective  time  constant  t'  for  cat. 

c 

The  values  listed  in  Table  5.1  have  standard  deviations 
ranging  from  30%  to  50%. 

The  values  from  our  cat  of  f„,  g , C„  = f„  + g„  and  T 

0 0 0 0 ^0  L 

are  plotted  in  Figure  5.7  as  a function  of  Apparently, 

f„  and  g reach  constant  saturation  values  above  0_  = 9 deg/sec. 
0 0 ^ 

• • 

The  mean  value  of  the  gain  0 /H  is  0.82.  This  is  somewhat 

s 

smaller  than  the  VOR  reflex  gain  of  0.9  which  is  adopted  in 
our  model  (see  section  3.5)  but  is  within  the  normal  range  for 
cat . 

The  values  used  in  the  model  for  f^^  and  g^  are  picked 

in  consideration  of  responses  with  step  and  sinusoidal  inputs 

in  both  deterministic  and  stochastic  modes. 

Table  5.1  gives  no  information  in  the  range  of  eye 

velocity  below  8.9  deg/sec.  Details  in  this  range  are  left 

to  a consideration  of  sinusoidal  responses  which  contain 

such  low  velocities.  To  implement  the  decrease  of  T as  6 

L s 

increases,  we  have  incorporated  a special  scheme  in  our 
model  to  be  discussed  shortly. 

Next,  we  derive  the  equation  by  which  C(t)  may  be  gen- 
erated from  (})(t)  (which  is  proportional  to  9^)  • The  deter- 
ministic WHERE  curve  C(t)  may  be  viewed  as  the  sum  of  two 
functions,  one  denoted  by  f (t)  which  appears  to  be  the  response 
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Figure  5.7 


This  figure  shows  the  plot  of  f^,  g^, 

Cq,  and  as  functions  of  the  s low- 
phase  eye  velocity  6^.  For  the  meanings 
of  f , g , C , and  T refer  to  Figures  5 

U 0 0 Li 

and  5.6.  Our  model  assumes  that  f^,  g^, 

and  C = f + g are  constants,  while  T 
U U U L 

the  time  constant  of  the  first-order  lag 
system,  decreases  as  the  input  head 
velocity  and,  hence,  the  slow-phase  eye 


velocity  increases. 


of  a first-order  system  to  a step  input,  and  the  other  denoted 
by  g(t)  which  may  be  considered  as  a step  function.  That  is. 


C(t)  = f(t)  + q(t) 


(5.3.1) 


From  Figure  5.5,  f(t)  is  in  the  form 


(5.3.2) 


with  the  time  constant  T (which  we  assume  to  be  constant 

J-i 

for  the  moment).  Then,  g(t)  has  a constant  value  g^  or 

g(t)  = g^  (5.3.3) 


We  note  that  for  t >>  T 

L 

?(t)  =-  * ■!„  • C„  (5.3 

Thus,  the  transfer  function  between  f(t)  and  4>(t)  is 
in  the  form  of 


f (s)  ^ G_ 

(Ms)  sT  + 1 

Lj 


(5.3.5) 


with  the  gain  G.  The  differential  equation  corresponding 
to  5.3.5  with  zero  initial  condition  is 


f (t) 


~ f(t)  + ^ (jXt) 
L L 


(5.3.6) 


The  solution  of  (5.3.6)  is  given  by 


93 


If  we  assume  for  a step  response  that  ())(T)  is  approximately 
a constant  denoted  by  <}),  then  (5.3.7)  reduces,  with  t^^  = 0 
and  f (t^)  = 0,  to 


f(t)  = G(J) 


(5.3.8) 


which  is  the  same  as  (5.3.2)  with  G<{>  = f^. 

For  the  purpose  of  further  analysis  and  computer  simu- 
lation we  want  to  put  this  into  a time  discrete  form.  Con- 
sider periodic  sampling  every  T seconds  or  at  t = kT  for 
k = 0,  1,  2,...,.  In  (5.3.7),  let  t = (k  +1)T  anc  t^  = kT 


Then, 


fl(k  + DTI 


.(k  + DT 
1 


[ (k  + 1)T-T] 


f(kT)  + 


— (})([)dT 


(5.3.9) 


If  we  assume  in  (5.3.9)  that 


(p(r)  = 4>(kT)  for  kT  < T < (k  + 1)T  (5.3.10) 

Then 


in  which,  for  simplification 

^k+1  = 

D = f(kT) 
k 

= <t)(kT)  (5.3.12) 
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are  used.  (5.3.11)  is  the  discrete  equation  corresponding 
to  (5.3.7). 

Now  it  is  necessary  to  deal  with  this  system's  non- 
linearities.  First,  roaches  a saturation  value  at 
4)  = 9 deg/sec.  Second,  T decreases  as  <J'  increases  from 

J_i 

0.9  seconds  at  (})  = 9 deg/sec  to  0.3  seconds  at  <{)  = 46.5  deg/sec. 
To  model  these  characteiistics,  the  arrangement  as  shown  in 
Figure  5.8  is  adopted  for  the  generation  of  the  deterministic 
WHERE  curve.  The  upper  channel  corresponds  to  5.3.11  if  the 
non-linear  limiter  is  short-circuited.  The  limiter  clips 
the  output  at  ^ f^(f^  > 0)  if  the  input  to  the  limiter 

exceeds  the  iDounds  +f^.  A possible  danger  with  the  type  of 
system  in  which  a first-order  lag  is  simply  followed  by  a 
limiter  is  that  some  internal  variable  such  as  f^^  in  (5.3.11) 
should  not  be  allowed  to  jsecome  large  compared  with  f'  for 
then,  when  input  (Ji  and,  hence,  G({)  suddenly  drops  bac)c  to  zero, 
the  system  will  momentarily  hang-up  in  saturation  until  f^^ 
has  time  to  drop  below  f'  due  to  the  effect  corresponding 
to  the  first  term  of  (5.3.7),  unless  properly  handled. 

Thus,  system  dynamics  must  be  removed  when  it  is  in 
saturation  and  this  is  guaranteed  by  the  location  of  the  non- 
linearity within  the  feedbac)c  loop  in  Figure  5.8.  Thus,  if 
the  system  is  in  saturation  (for  ())  > 0)  but  falls  bac)c 
within  the  bounds  of  the  limiter,  the  output  also 


95 


Figure  5.8 


This  figure  shows  a mathematical  scheme 

to  generate  the  deterministic  WHERE  curve 

The  upper  branch  creates  the  portion 

of  the  deterministic  WHERE  curve  which 

is  the  output  of  the  first-order  system 

to  a step  head  velocity.  The  lower 

branch  creates  the  step  portion  of  the 

curve.  With  this  arrangement,  if  the 

output  is  in  saturation  (+f , and  if 

the  input  falls  in  such  a way  as  to  make 

fall  within  the  bounds  of  the 

limiter,  the  output  f,  , also  starts 

k+1 

falling  immediately.  This  is  illustrated 
at  point  a.  Also  with  this  scheme, 
the  effective  time  constant  T of  the 

Li 

system  decreases  as  the  input  <t>  increases 
due  to  the  increasing  clipping  action  of 


the  limiter. 


starts  to  fall  at  once.  This  follows  from  the  fact  that,  if 
one  assumes  in  (5.3.11)  that  f^^  = f'  and  G4>j^  < f'f  then  it 
can  be  shown  that  f,  , < f^,  that  is,  the  output  is 
decreasing. 

However,  the  main  point  of  the  limiter  in  the  upper 
channel  is  that  it  causes  the  effective  "time  constant"  to 
decrease  as  the  input  <p  increases  due  to  the  increasingly 
earlier  clipping  action  of  the  limiter.  For  reasons  discussed 
later,  f*  - Iq  = 1^  decj  is  .diosini  tentatively. 

The  lower  channel  of  Figure  5.8  is  self-explanatory, 

consisting  solely  of  a non-linear  limiter  with  g ~ ~ 

= (b,  . For 
k 

\+i  = V - ^k-.i  = 

Finally,  C,  , , = f , . , + i is  equivalent  to  C(t),  the 
k+1  k+1  k+1 

deterministic  WHERE  curve.  Referring  to  Figure  5.8,  we  note 
that 


corresponding  to  (5.3.3),  For  |4<|^1  < g^<  *^k+l 


max 


k+1 


5 . 4 The  WHEN  curve 

Our  cat's  VOR  responses  from  various  step  input  head 
velocities  indicate  that  the  steady  state  deterministic 
WHEN  curve  coincides  approximately  with  the  zero  degree 

line  (Primary  position) . For  6^  > 30  deg/sec,  this  requires 
a threshold  value  Q of  12  deg  which  also  gives  an  acceptable 
quick-phase  frequency  from  the  model.  Our  data  also  shows 
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that  at  6^  - 10  we  need  Q = 6 to  match  the  quick-phase 

I frequency  of  our  model  with  that  of  our  cat.  Since  we 

want  Dq  for  the  step  input  of  4>  = 10  (with  ())  being  the  SCC 

output  signal)  to  coincide  with  the  zero  degree  line,  the 

1 

magnitude  of  the  steady  state  WHERE  curve  must  be  equal 

to  six,  that  is  |Cq|  = Q = 6.  Since  = G<})  + g,  we  set 

g|())|  + 1 = lOG  + 1 = 6 (5.4.1) 

in  which  G is  the  gain  of  the  transfer  function  5.3.5. 

From  5.4.1,  we  obtain  G = 0.5. 

Next,  we  want  to  find  |(}>|  at  which  the  deterministic 
WHERE  curve  saturates.  For  this  purpose,  since  <(>  would  be 
greater  than  one,  we  set 

g|4)|  + 1 = 0.5l4)|  + 1 = 12  (5.4.2) 

t ' 

from  which  we  obtain  |(t)|  =22.  Thus,  we  make  Q = 12  for 
|(j>|  ^22.  For  10  < 1(})|  £ 22,  we  set  Q = 0.5|4)|  +1.  This 
will  make  the  line  coincide  with  the  zero  degree  line  for 
any  step  input  head  velocity  which  results  in  |<))|  ^10  with 
the  understanding  that  |4>|  - |9gl  because  0^  = O.Otfi  as  a 
result  of  the  reflex  gain. 

For  0 < |(j)|  < 10,  we  do  not  have  data  to  determine  the 
I value  of  Q for  this  region.  However,  the  simulation  of  our 

model  indicates  that  a constant  value  of  Q = 6 for  this 
region  yields  acceptable  sinusoidal  response  for  the  test 
I frequencies  of  0.05,  0.25,  and  1.2  Hz,  with  peak  velocities 

I 

[ 
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of  10,  20,  and  30  deg/sec.  For  this  reason,  we  have  adopted 
Q = 6 for  0 < 4>  < 10.  The  summary  is  given  below: 

Q = 12  for  |<f|  >22 

Q = 0.5l4)|  + 1 for  22  ^ |(})|  > 10  (5.4.3) 

Q = 6 for  10  ^ 1 <()  I >0 
Equation  5.4.3  is  plotted  in  Figure  5.9. 

If  Q given  by  Equation  5.4.3  is  used  in  the  model, 
for  |4>|  > 10  would  coincide  with  the  zero  degree  line.  For 
10  ^ |(}>|  ^ 0,  Dq  would  shift  in  the  compensatory  direction 
(opposite  to  the  direction  of  the  head  motion)  by  the  amounts 
explained  below  (for  left  turn  head  motion) . 

For  1 ^ 14>1  > 0,  Dq  = -(0.5|(})|  + 1<))1)  +6  = 6-  1.5l(})l. 

For  10  >$  > 1,  = -(0.51$!  + 1)  + 6 = 5 - 0.5l({)|. 

The  plot  of  Dg  vs  )0|  is  given  in  Figure  5.10.  Note  that 
is  the  "steady  state"  deterministic  WHEN  curve  for  step 
input  head  velocity,  which  is  a constant  for  given  |(()|. 

Next,  we  consider  the  condition  for  the  quic)c-phase 
initiation.  Referring  to  Figure  5.6  (in  section  5.-3)  for 
left  turn  step  input  head  velocity,  the  equation  for  the 
deterministic  WHEN  curve  is  C(t)  + Q.  The  efference  copy 
6'(t)  of  the  eye  position  during  the  slow-phase  always 
satisfies 


Figure  5.9 

This  plot  shows  the  threshold  value  Q i 

in  degrees  as  a function  of  the  magni- 
tude |4>|  of  the  see  output  signal.  4> 
has  units  of  deg/sec.  The  dependence 

of  Q on  stimulus  intensity  yields  not  I 

only  acceptable  QP  frequencies  in  step 
responses,  but  also  reasonable 
sinusoidal  responses. 


101 


Figure  5.10 

This  plot  shows  (for  left  turn  head 
motion)  the  steady  state  value  of 
the  deterministic  WHEN  curve  (for  step 
input  head  velocity)  as  a function  of 
the  magnitude  1(})|  of  the  SCC  output 
signal.  It  shows  that  = 0 for 
I0I  >10.  However,  for  0 < 1(}>1  < 

10,  Dq  is  shifted  to  the  compensatory 
direction  (opposite  to  the  direction 
of  head  motion)  by  the  amount 
indicated  in  the  figure. 
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Figure  5.11 

This  figure  shows  the  threshold  device 
as  used  for  the  generation  of  quick- 
phases.  Whenever  A0,  which  is  the 
distance  from  the  efference  copy  6^(t) 
of  the  eye  position  to  the  WHERE  curve 
C(t),  exceeds  the  thresholds  a 
WHEN  signal  is  generated  and  intiates 
a process  (if  it  meets  the  50  milli- 
seconds refractroy  period)  which 
results  in  a generation  of  a quick- 
phase  as  explained  in  the  text. 

The  value  of  Q is  determined  according 
to  5.4.3  and  Figure  5.9. 
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TO  THE  EXTENTEO 
FINAL  COMMON 
PATHWAY 


Figure  5.11  The  Threshold  Device  for  Quick-Phases 


Thus  a quick-phase  is  initiated  whenever 

0'(t)  > C(t)  + Q (5.4.5) 

Now,  referring  to  Figure  5.5  (in  section  5.3)  for  right 
turn  step  head  input  velocity,  the  deterministic  WHEN  curve 
is  C(t)  - Q.  For  this  case,  0^(t)  during  the  slow-phase 
always  satisfies 

6'(t)  > C(t)  - Q (5.4.6) 

It  follows  that  the  quick -phase  is  initiated  whenever 

e"(t)  < C(t)  - Q (5.4.7) 

5 . 5 The  threshold  device 

The  condition  for  the  quick-phase  initiation  is  from 
(5.4.5)  and  (5,4.7). 

|c(t)  - e'(t) I > e (5.5.1) 

Equation  (5.5.1)  is  the  unified  condition  for  QP  initiation 
for  both  left  and  right  turn  head  motion  which  sets  the 
threshold  for  the  quick-phase  at  Q as  shown  in  Figure  5.11. 
Whenever  the  error  angle  A6  from  0'(t)  to  C(t)  matches  or 
exceeds  tlie  thresholds  ^ Q,  a quick-phase  is  initiated.  As 
long  as  A6  stays  within  the  boundaries  + Q and  - Q,  no  quick- 
phase  is  initiated. 

Referring  to  Figure  5.11,  a WHEN  signal  generated  by 
the  threshold  device  actuates  the  AND-gate  if  it  also  meets 
the  50  milliseconds  refractory  period  condition  (see  sec- 
tion 7.4).  Once  the  AND-gate  is  actuated,  s is  closed 

2 
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through  the  OR-gate.  Once  s is  closed,  it  remains  closed 

2 

as  long  as  the  output  from  the  amislifier  persists  because 
of  the  feedback  signal  from  the  amplifier.  The  threshold 
device  is  also  referred  to  as  the  WHEN  function  because  it 
generates  the  WHEN  signal  from  the  error  signal.  It  is  to 
be  noted  that  once  the  threshold  value  is  exceeded  (from 
either  direction  for  whatever  reasons) , the  eye  is  always 
driven  to  C(t)  by  a QP. 

5. 6 The  total  integrated  deterministic  VOR  model 

The  total  integrated  deterministic  (without  noise)  VOR 
model  is  shown  in  Figure  5.12  with  each  block  being  explained 
in  previous  sections.  The  final  values  of  parameters  of 
each  block  in  the  model  will  be  determined  after  the  overall 
considerations  of  the  responses  to  step  and  sinusoidal  inputs 
and  the  behavior  of  the  stochastic  model.  These  are  discussed 


in  a later  chapter. 


Figure  5.12 

This  figure  shows  the  total  integrated 

deterministic  VOR  model,  in  which  the 

mechanism  of  quick-phase  generation 

is  integrated  with  the  compensatory 

slow-phase  model.  During  the  slow- 

phase,  s is  closed  and  s is  open. 

1 2 

During  the  quick-phase,  s is  open 

and  s is  closed. 

2 
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Fiqure  5.12  The  Total  Integrated  Deterministic  VOR  Model 
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CHAPTER  VI 


STOCHASTIC  MODEL 


6 . 1 Introduction 

The  stochastic  model  is  necessitated  by  the  random 
behavior  of  the  nystagmus  (more  specifically  by  the  fluctu- 
ations of  the  WHEN  curves  and  the  WHERE  curves)  as  illus- 
trated in  Figure  5.4.  We  have  observed  that  as  the  eye 
velocity  increases,  fluctuations  are  more  pronounced. 

We  view  the  quick-phase  end-points  as  the  scunpled 
sequence,  C*(t),  of  the  continuous  WHERE  curve,  C(t).  Sim- 
ilarly, we  view  the  quick-phase  initiation  points  as  the 
sampled  sequence,  D* (t) , of  the  continuous  WHEN  curve,  D(t). 
Based  on  several  VOR  responses  of  our  cat  corresponding  to 
various  nominal  step  input  head  velocities,  relevant  measure- 
ments are  made  and  statistics  are  computed.  The  mean  values 
of  various  relevant  quantities  for  each  nominal  head  velocity 

are  shown  in  Table  6.1,  in  which  H„,  H,  and  9 have  the  same 

0 s 

meaning  as  in  Table  5.1  in  section  5.3  and  0^  is  the  mean 
quick-phase  velocity  in  deg/sec  obtained  by  dividing  the  mean 
of  all  quick -phase  magnitudes  of  several  runs  for  a given 
H by  the  mean  of  all  quick-phase  durations  of  several  runs 
for  the  same  ft.  Q is  the  separation  in  degrees  between  the 
steady  state  deterministic  WHEN  curve  and  WHERE  curve.  f^ 
is  the  mean  quick-phase  frequency  in  beats  per  sec.  o(D*) 
is  the  standard  deviation  of  the  s. imp)  i d secjuence  D*(t)  of  the 
WHEN  curve.  0(C*)  is  the  standard  deviation  of  the  sampled 
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Table 


sequence  C*(t)  of  the  WHERE  curve.  R is  the  linear  correla- 

! 2 

tion  coefficient  between  D* (t)  and  C*(t).  The  mathematical 
approximations  for  these  quantities  are  derived  as  follows. 

As  we  have  discussed  in  section  5.1,  the  WHERE  curve, 
C(t),  may  be  thought  of  as 

C(t)  = C{t)  + r (t)  (6.1.1) 

2 

in  which  C(t)  is  the  mean  location  of  the  deterministic 
WHERE  curve  and  r^  (t)  is  its  noise  with  zero  mean.  Similarly 
the  WHEN  curve,  D(t),  may  be  thought  of  as 

D(t)  = D(t)  + r (t)  (6.1.2) 

i 

in  which 

D = C(t)  + Q (6.1.3) 

is  the  mean  deterministic  WHEN  curve  and  r (t)  is  its  noise 

1 

with  zero  mean.  This  line  of  thought  suggests  introducing 
noise  sources  at  certain  locations  in  the  model.  The  noises 
thus  introduced  are  generated  through  filters  from  purely 
random  (white)  Gaussian  noise  sources  and  located  in  such 
places  as  to  create  noisy  WHEN  curves  and  a noisy  WHERE 
curve  with  certain  statistical  characteristics  such  as 
standard  deviations  and  a linear  correlation  coefficient, 
and  with  certain  spectral  compositions  which  would  reasonably 
match  those  of  our  cat's  VOR  step  responses.  We  have  chosen 
Gaussian  inputs  to  the  filters  because  that  distribution  is 
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mathematically  tractable  and  the  noise  histograms  constructed 
from  the  limited  data  available  from  our  cat’s  step  responses 
seem  sufficiently  bell-shaped  to  permit  this  approximation. 

6 . 2 The  stochastic  model 

Let  X (t)  and  z (t)  be  two  independent  Gaussian  noise 
2 2 

signals  both  with  zero  means  and  with  desired  variances  and 

spectral  components.  How  these  noises  are  generated  is 

discussed  in  section  6.3.  x (t)  and  z (t)  are  introduced 

2 2 

into  the  VOR  model  at  points  p and  p respectively,  as  shown 

1 2 

in  Figure  6.1. 

The  point  p is  selected  for  z (t)  because  noise  injected 
2 2 

at  this  point  would  simulate  the  resultant  effect  of  combining 
all  the  noise  sources  from  anywhere  in  the  WHERE  function 
and  the  neural  integrator  circuits.  It  is  reasonable  to 
assume  that  these  circuits  are  not  noise  free. 

Now,  referring  to  Figure  6.1,  we  note  that 

A6  = [c(t)  - e'(t)j  + z^(t) 

= [^C(t)  + - e'(t) 

= C(t)  - 0"(t)  (6.2.1) 

According  to  (6.2.1),  z (t)  may  be  treated  as  noise  on  either 

2 

C(t)  or  on  [C(t)  - 6'(t)). 
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Figure  6.1 

Noises  are  introduced  into  the  model  to 

create  a noisy  WHERE  curve  and  a noisy 

WHEN  curve.  Noise  z (t)  is  added  to  the 

2 

deterministic  WHERE  curve  C(t)  at  p to 

2 

create  the  noise  on  the  WHERE  curve. 

Another  noise,  x (t) , is  added  at  p so 
2 1 

that  the  resultant  noise  on  the  threshold, 

thus  on  the  WHEN  curve,  is  x (t)  + z (t) . 

2 2 

This  would  cause  correlation  between  the 

noises  on  the  two  curves  because  z (t)  is 

2 

common  in  noises  of  both  curves. 
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Next  we  consider  the  noise  r (t)  on  the  deterministic 

1 

WHEN  curve.  One  way  to  interpret  the  fluctuations  of  the 
WHEN  curve  is  to  attribute  them  to  the  fluctuations  of  the 
threshold  value  of  Q.  It  is  neurophysiologically  reasonable 
to  assume  that  the  decision  process  in  the  brain  is  so  com- 
plex and  unpredictable  that  it  can  be  represented  by  the 
addition  of  a random  signal.  For  this  noise,  a separate 

noise  source  x (t)  is  introduced  at  p in  addition  to  z (t) 

2 1 2 

introduced  at  p (Figure  6.1).  The  input  signal  A0  to  the 
2 T 

threshold  device  is  now 

A6  (t)  = C(t)  - 0^(t)  + X (t)  + z (t) 

.p  2 2 

= C(t)  - 0"(t)  + r (t)  (6.2.2) 

1 

so  that  X (t)  + z (t)  = r (t)  by  definition. 

2 2 1 

With  this  arrangement,  r (t)  is  correlated  with  r (t) 

1 2 

because  the  expected  value  of  r r expressed  by  E(r  r ) is 

12  12 

given  by 

E(rr)=E  (x  + z)z 
1 2 2 2 2 

= E (x  z ) + E (z  ^ ) 

2 2 2 

= E(Z  ^) 

2 

= (6.2.3) 
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since  x and  z , both  with  zero  means,  are  uncorrelated. 

2 2 

In  the  model,  the  noise  z (t)  would  cause  correlation  between 

2 

the  WHEN  curve  and  the  WHERE  curve  because  z (t)  affects 

2 

each  of  them  through  the  high-gain  amplifier  (WHERE  curve) 
and  the  threshold  device  (WHEN  curve) . This  correlation  was 
first  noted  by  casual  inspection  of  the  raw  data  and  then 
confirmed  by  measurement.  It  is  intellectually  satisfying 
that  this  model  is  not  only  compatible  with  such  a correla- 
tion, but  it  would  be  difficult  and  unrealistic  to  introduce 
noise  in  such  a way  that  there  was  no  correlation. 

Next,  we  want  to  derive  the  expression  for  the  linear 

correlation  coefficient  R between  the  eye  position  at  the 

12 

beginning  (t  ) and  the  eye  position  at  the  end  (t  ) of  each 
1 2 

quic)^ -phase.  From  (6.1.1),  (6.1.2),  and  (6.1.3)  we  have 

C*(t  ) = C(t  ) + r *(t  ) (6.2.4) 

2 2 2 2 

D*(t  ) = C(t  ) + Q + r *(t  ) (6.2.5) 

1 1 — 11 

in  which  the  superscript  * denotes  the  sampled  values  either 

at  the  beginnings  or  at  the  ends  of  quic)c -phases.  From 

(6.2.4)  and  (6.2.5),  noting  that  E(r*)  = 0 and  E(r*)  = 0, 

1 2 

we  obtain 


E(C*)  = C 


(6.2.6) 


E(D*)  = C + Q 


(6.2.7) 
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Thus,  remembering  that  t is  at  the  beginning  and  t is  at 

1 2 

the  end  of  each  quick -phase,  the  numerator  of  R is  given  by 


E I - E^D*(t^)j  C*(t^)  - E^C*(t^)j  I 

= E|r*(t  )r*(t  ) 

1 1 2 2 I 


and  the  denominator  of  R is  given  by 

1 2 


(6.2.8) 


(.1 

D*(t  ) 

- E/o*(t  )\ 

2 

E 

1 

Ie  ( 

L 

'r*(t  )y] 

\ 

E 

* 

(r*{t 

■ 

i 1 1/1 

\ 2 

2/ 

I 

1 

- 

' I 

C*(t  ) 
2 


- E(cMt^)) 


(6.2.9) 


It  follows  that 


(6.2.10) 


which  is  identical  with  the  definition  for  the  linear  cor- 
relation coefficient  between  r*  at  t (the  beginning)  and 

1 1 

r*  at  t (the  end)  of  each  quick-phase.  Equation  (6.2.10) 

2 2 

is  utilized  in  section  6.4.  Note  that  (6.2.10)  may  also  be 
approximated  by 


"k^k''=l  ( ^ 


*(k  T) 


r*(k  't) 
2 


1 2 


^ J [rJ(k'T)' 


I 

n 


^ r*(k''T) 


k =1 


k =1 


(6.2.11) 


in  which  T is  a periodic  sampling  interval  (see  section  6.3)  and 
in  which  k^T  corresponds  to  t^  in  (6.2.10)  and  k'^T  corresponds 
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to  t in  (6.2.10).  R approaches  R asymptotically,  (if 
2 12  12 

it  is  ergodic)  as  n increases. 

The  resulting  stochastic  VOR  model,  which  is  shown  in 
Figure  6.2,  is  acheived  by  simply  incorporating  the  scheme 
of  Figure  6.1  into  the  deterministic  VOR  model  as  presented 
in  Figure  5.12. 

6 . 3 Noise  generation 

We  want  the  output  from  our  stochastic  model  to  match 
the  random  behavioral  data  from  our  cat. 

Knowing  the  variance  of  the  noise  signals  is  not  enough. 

It  is  important  to  know  the  spectral  composition  of  the 
noise.  It  is  likely  that  the  actual  neural  noise  is  band- 
limited,  and  computer  simulations  in  which  very  high  or  very 
low  frequencies  were  allowed  gave  unrealistic  behavior. 

For  this  reason,  we  attempted  to  determine  the  spectral 
composition  of  the  WHERE  and  WHEN  curves  of  our  cat's  VOR 
responses  by  frequency  analysis  of  the  sampled  curves  C* (t) 
and  D*(t).  Since  the  inter-saccadic  interval  which  corre- 
sponds to  the  inter-sampling  period  is  not  uniform,  we  could 
not  use  the  discrete  Fourier  transform  which  assumes  a uni- 
form inter-scimpling  period.  We  overccime  this  difficulty 
by  using  a trigonometric  interpolation.  In  this  interpolation, 
a finite  Fourier  series  is  forced  to  pass  through  data 

values  X (with  n=l,2,---,  N)ofN  data  points,  which 
n 

may  not  be  equally  spaced,  by  proper  choices  of  the  co- 
efficients of  the  truncated  Fourier  series. 
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Figure  6.2 

The  stochastic  VOR  model  is  achieved 

by  simply  incorporating  the  scheme 

of  Figure  6.1  into  the  deterministic 

model  as  presented  in  Figure  5.12. 

This  model  will  create  the  WHERE 

curve  with  noise  r (t)  = z (t) 

2 2 

and  the  WHEN  curve  with  noise 

r (t)  = X (t)  + z (t)  . 

12  2 
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For  an  even  number  of  data  points  N with  the  data  length 
t' , the  value  of  x(t)  which  is  the  deviation  from  the  mean, 

1j 

for  0 < t < t'’,  may  be  expressed  as 


© (?) 

X(t)  = ^ A^^COS  qw^t)  + ^ B^^sin  qw^t') 


q=l 

in  which  w = 2'n/T' 


^ ^ (6.3.1) 

For  N data  points,  we  have  N equations 


with  N unknowns,  N/2  A 's  and  N/2  B 's,  to  be  solved.  The 

q q 

magnitude  F^  at  a particular  angular  frequency  (in  radians/sec) 

of  qw  = 2TTq/T^  i s 
o L 


q V q c 


(6.3.2) 


Unfortunately,  the  length  of  th®  records  from  our  cat's 
step  response  is  typically  only  about  5 to  6 sec  long  with 
each  record  consisting  of  about  a dozen  data  points.  The 
short  duration  of  the  record  is  inherently  imposed  by  the 
leaky  integration  of  the  SCC,  which  causes  the  slow-phase 
eye  velocity  to  gradually  decrease  with  an  effective  time- 
constant  of  about  12  sec. 

For  T^  = 5 and  N 12  in  (6.3.1),  the  lowest  angular 
L 

2ti 

frequency  (in  the  available  frequency  range)  is  w = — rad/sec 

1 5 

(2T\  \ 

(0.2  Hz)  for  q = 1 and  the  highest  frequency  is  w^  = \^J  6 
rad/sec  (1.2  Hz)  for  q = 6.  The  result  of  spectral  analysis 
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by  this  method  for  these  data  was  very  erratic.  However,  the 
results  showed  no  consistent  trend  in  the  power  spectrum  over 
this  bandwidth,  suggesting  that  in  the  range  of  roughly  0.2 
to  1.2  Hz,  there  is  no  reason  to  assume  the  spectrum  to  bo 
other  than  flat.  The  Fourier  analysis  described  above  would 
give  no  information  about  frequency  components  outside  this 
frequency  range. 

It  turns  out  that  the  plot  of  the  computer-simulated 
noise  with  components  of  w greater  than  fluctuates  so 
rapidly  that  it  presents  an  unphysiological  appearance,  while 
the  plot  of  the  noise  with  components  of  w less  than  w.j  is 
undesirable  because  it  has  long  periods  in  which  QP's  are  all 
too  short  or  too  long. 

Under  these  circumstances  and  constraints,  we  have 
designed  a filter  which  would  produce,  from  a purely  random 
Gaussian  sequence,  another  Gaussian  sequence  v.ith  flat  band- 
pass characteristics  (with  gain  of  1)  for  < w < w^ 
and  rejects  frequencies  above  or  below  this  range. 

The  transfer  function  of  the  filter  between  the  output 
noise  y(t)  and  the  input  noise  u(t),  which  is  the  purely 
random  Gaussian  noise  with  zero  mean,  is  given  by 


y (s) 
u (s) 


ST 

1 


(3T_  . 1)  (ST^  + 1) 


(i.3.3) 


with  T 

1 

271 
5 


w 


1 from  w T 
1 1 

1 and  w = 

2 


nd  T = 1/8  from  w T =1  using 
2 2 2 

6=8  for  computational  simplicity. 
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As  the  first  step  to  obtain  the  time-discrete  equation. 


we  define  x(s)  such  that  y(s)  = T^sx(s)  so  that 


x(s) 


u(s)  (sT^  + 1)  (sT^  + 1) 


(6.3.4) 


The  differential  equation  corresponding  to  (6.3.4)  with  zero 
initial  conditions  is  given  by 


T + T _ 

x(t)  = - — - x(t)-  :pj—  x(t)+ 

12  12  12 

(6.3.5) 

With  state  variable  notation  of  x (t)  = x(t)  and  x (t)  = x(t) 

1 2 

which  implies  that  y(t)  = T x (t) , (6.3.5)  may  be  expressed 

1 2 

in  a matrix  form  as 


0 1 \ /x^  (t) 


a a / \ X (t) 
2 1 22/  \ 2 


(■*-) 


u (t) 
2 


in  which 


21  T T 

1 2 


(6.3.6)  is  in  the  form  of 


a = - 
22 


T + T 
1 2 


T X 
1 2 


(6.3.6) 


, and  u (t)  = u(t) 
2 


(6.3.7) 


X(t)  = AX(t)  + BU(t) 


(6.3.8) 
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in  which 


X(t) 


(X  (t) 

' 

X (t) 
2 


A = 


(6.3.9) 


(6.3.10) 


B = -a 

2 1 


and 


(6.3.11) 


(6.3.13) 

in  which  as  the  transition  matrix  from  X(t  ) to  X(t). 

0 

Now,  consider  pericxJic  sampling  at  every  T sec  (compu- 
tation interval)  or  at  t = )cT  for  )c  = 0,  1,  2,  ---.  We 
see)i  to  find  the  output  of  the  filter,  y[()<  + 1)T),  given 
the  input  U(k:T)  and  the  filter  properties  A.  Substitu*-ing 

t = )cT  and  t = (It  + 1)T  in  (6.3.13),  we  obtain 
0 
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X [(k  + 1)tJ  = ^lT)X(kT)  + W^(kT)  (6.3.14) 

where  the  transition  matrix  <1>(T)  (because  <t>(k+l)T  - kT]  = <J>(T)) 
is  given  by 


'4)  (T) 

1 1 


4>  (T) 

1 2 


(J)(T)  = e 


AT  A 


2 1 


(T) 


<t>  (T) 

22 


(6.3.15) 


and 


W^(kT)  = 


r 

kT 


(k  + 1)T 


A[(k  + 1)T  - T] 


BU(T)  dT 


w (kT) 

XI 


w (kT) 
X2 


(6.3.16) 


In  (6.3.16),  both  w (kT)  and  w (kT)  are  purely  random 

XI  X2 

Gaussian  sequences  with  zero  mean  because  of  (6.3.12)  in 

which  u (t)  is  a purely  random  Gaussiein  noise  with  zero 
2 

mean  by  assumption. 

Next,  we  determine  the  elements  of  <I>(T)  = exp  (AT)  by 
the  inverse  Laplace  transform  (L  S method.  That  is, 

AT  — — iT 

4>(T)  = e = L (sI  - A)  J (6.3.17) 

From  (6.3.10),  we  get 


(si  - A) = 


s 


-a 

2 1 


-1 


s-a 

22 


(6.3.18) 
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It  follows  that 


The  fact  that  (6.3.19)  is  the  correct  expression  may  be 
confirmed  by  multiplying  it  by  (6.3.18),  which  results  in 
the  identity  matrix  as  it  should.  Using  (6.3.7), 


- a s 
22 


a = [s  + — 
21  \ T 


s + 


(6.3.20) 


Finding  the  inverse  Laplace  transforms  of  (6.3.19)  with 

(6.3.20)  using  the  standard  table,  we  obtain 


(6.3.21) 


(6.3.22) 


(6.3.23) 
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and 


(6. 3.24) 


Equation  (6.3.14)  may  now  be  expressed  in  terms  of  (6.3.21) 


to  (6.3.24)  as 


X [ ()c  + 1)T1  = 4>  (T)x  ()cT)  + ())  (T)x  ()cT)  + w ()cT) 

1 111  122  XI 

(6.3.25) 


X [ (k  + 1)T]  = <p  (T)x  (kT)  + (p  (T)x  (kT)  + w (kT) 
2 211  222  X2 

(6.3.26) 


Because  y(t)  = t x (t)  as  mentioned  before, 

1 2 

yl(k  + 1)T]  = T X ( (k  + 1)T]  (6.3.27) 

1 2 

which  is  the  output  we  are  seeking. 

The  remaining  problem  is  to  determine  w (kT)  such  that 

X2 

E(y^)  = E(x  ^)  = (x  ) , for  T =1,  matches  the  variance 
2 2 1 

obtained  from  our  cat's  VOR  responses.  E denotes  the  opera- 
tion of  taking  the  expected  value  of  its  argument  (which  is 
equal  to  taking  its  time-average  if  it  is  ergodic) . 

In  (6.3.14),  X[(k  + 1)T]  is  a Gaussian  random  sequence 
if  W^(kT)  is  a Gaussiein  purely  random  sequence  and  the  initial 
state  vector  X(0)  is  Gaussian  (Bryson  and  Ho,  1975).  In  our 
case,  E[W^(kT)]  = 0 and  E[X(kT)]  = 0.  Next,  we  want  to 
generate  a two-dimensional,  jointly  Gaussian  sequence  w'(kT) 
by  a digital  conputer  such  that  it  is  equivalent  to  W^(kT) 
of  (6.3.16),  to  the  extent  that  the  statistical  characteristics 
of  the  output  X[(k  + 1)T]  of  (6.3.14)  forced  by  w'(kT)  is 
identical  with  those  forced  by  W (kT) . 

X 
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From  (6.3.16)  using  (6.3.15)  we  obtain  (denoting  B by  b) 


W (kT) 

X 


(k  + 1)T 


KT 


b l<t>  [ (k  + 1)T  - T]  u (t)  dT 
1 2 2 


(k  + 1)T 


h I Ip  [ (k  + 1)T  - T]  u (T)  dT 
22  2 

KT 


We  want  to  determine 


:|w^(kT)W^(kT)^  = 


'w^  (kT) 
xi 


(6.3.28) 


w (kT)w  (kT) 
XI  X2 


<w  (kT)w  (kT)  (kT) 

xi  X2  X2 


(6.3.29) 

with  u (t)  representing  a continuous  white  (purely  random) 

2 

noise,  that  is 


Elu  (t)u  (T)l  = ii>6(t-T) 
2 2 


(6.3.30) 


in  which  ip  is  the  power  spectral  density  of  the  white  noise 

u (t) , and  6(t-T)  is  the  Dirac  delta  function.  Substituting 
2 


( 


130 


from  (6.3.28)  into  (6.3.29)  and  using  (6.3.30),  we  obtain 


4)^  (T-T)dT 
12 


4>^  (T  - T)  dT 
22 


(6.3.31) 


(6.3.32) 


E w 

()cT)w  ()cT) 

= ipb^ 

f 

4)  (t-t)4) 

(T-T)dT 

X2  J 

1 12 

22 

(6.3.33) 

in  w)iic)i  (t)  (T-T)  and  (p  (T-X)  are  given  by  (6.3.22)  and 

12  22 

(6.3.24)  witli  T replaced  by  the  new  argument  (T-T).  Carrying 
out  integrations  in  (6.3,31),  (6.3.32),  and  (6.3.33),  we  have 


1 

( 


i 


(6.3.34) 
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I 

1 


- e 


(6.3.36) 

where  T and  x are  originally  introduced  in  (6.3.3)/  and 
l 2 

the  computation  interval. 
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with  T = 1,  T = 0.125,  and  T = 0.01  in  the  above  three 
1 2 

equations,  we  obtain 


Q = eFw  (kT)W^(kT)l  = 
X L X X J 


w w 
11  12 


w w 
12  22 


_5  _3 

1.9918  X 10  2.9258  x 10 


2.9258  X 10  0.58556 


(6.3.37) 


Suppose  we  have  another  covariance  matrix  given  by 


q q 

Xil  ^Xi2 


q q 

Xl2  X22 


po  0 
1 2 


pa  a o 

1 2 2 


(6.3.38) 


If  0 ^ is  spiecified,  0 may  be  related  to  0 by 
2 1 2 


o ^ = c^o  ^ 


(6.3.39) 


in  which  the  value  of  is  to  be  determined  subsequently.  This 


allows  us  to  express  (6.3.38)  by 


C^O  ^ PCO  ^ 

2 2 


2 2 
pco  ^ O 
2 2 


(6.3.40) 


We  want  to  make  (6.3.40)  equal  to  (6.3.37),  that  is,  we  want 
= q\  This  requires  that 


w ij;  = 0 ^ (6.3.41) 

22  2 


w = C^O  ^ = C^w  ij;  (6.3.42) 

11  2 22 


w ijj  = pCO  ^ = pCw  (6.3.43) 

1 2 2 22 


It  follows  that,  cancelling  iji's  out. 


w 


= 

1 1 

(6.3.44) 

« L, 

w 

22 

pc  = 

W 

1 2 

(6.3.45) 

W2  2 

1 

L 1 2 _ 

w 

12 

(6.3.46) 

p = 

C w - 

22 

/ w w 
\ 1122 

For  T = 1,  T 
1 2 

= 0.125, 

and  T = 0.01,  we  get 

= 

3.4016  X 

O 

1 

C = 

5.8323  X 

-3 

10 

pc  = 

4.9966  X 

-3 

10 

(6.3.47) 

p = 0.8567 


With  the  knowledge  of  the  covariance  matrix  given  by 
(6.3.40),  w^  (kT)  and  w^  (kT)  are  generated  by  a computer 
routine  (Scheuer  eind  Stoller,  1962)  such  that  the  joint 
probability  density  of  w^  and  w^^  corresponds  to  the  jointly 
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Gaussiaii  density  function  j (Papoulis,  1965)  given  by 

(since  w and  w have  zero  mean) 

XI  X2 

f /w  , w \ = ^ r 

V XI  ' y.2  ) 2TIC7  0 (1  - p'^)  “5 

1 2 

2pw  w w ^ 

X]  X2  ^ X2 

o a 0 ^ 

12  2 _ 

(6.3.48) 

in  which  O = o (w  ) , O = O (w  ) , and  O = CO  . 

1 Xi  2 X2  12 

T 

W (kT)  = w (kT)  w (kT)  thus  generated  becomes  the 

X [ XI  X2  J 

input  of  (6.3.14).  Since  W (kT)  is  related  to  u (T)  by 

X 2 

(6.3.28)  and  u (T)  is  assumed  to  be  a purely  random  (white) 

2 

Gaussian  noise,  y((k  + 1)T]  = T x [ (k  + 1)T]  given  by  (6.3.27) 

1 2 

would  have  the  same  band-pass  characteristics  as  the  output 

of  the  filter  equation  given  by  (6.3.3). 

Next,  we  want  to  determine  the  relationship  between 

E[x^(kT)]  = O^  and  E [W^  (kT) ] = o To  do  this  (Bryson 

2 X2  X2  WX2 

and  Ho,  1975),  we  start  with  (6.3.8)  which  is  repeated 
below  as 

X(t)  = AX(t)  + BU(t)  (6.3.49) 

Since  u (t)  of  U(t)  is  a Gaussian  random  process  with  zero 
2 

mean,  X(t)  is  also  a Gaussian  random  process  with  zero  mean 
because  linear  transformations  of  a Gaussian  vector  preserve 
its  Gaussian  character. 
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We  define  P(t)  by 


P(t)  = E[X(t)x'^(t)  ] 


in  which 


P (t)  P (t) 

11  12 


P(t)  = 


P (t)  P (t) 

12  22 


We  assume  that 

P(0)  = E[X(0)x'^(0)  ] 
is  a purely  Gaussian  random  variable. 
We  define  4*  by 


E[U(t)u'^(T)  ] = 4'(S(t-T) 

in  which 


Because  of  (6.3.12)  and  (6.3.30),  we  have  from 
above  two  equations. 


0 0 


E[U(t)U  (T)l 


6(t-T) 


0 iji 


in  which  ^ = ip  is  used.  Under  these  conditions, 

22 

be  shown  that 


P(t)  = AP(t)  + P(t)A"^  + 


with  P(t  ) given. 
0 


(6.3.50) 


(6.3.51) 


(6.3.52) 


(6.3.53) 


(6.3.54) 


the 


(6.3.55) 


it  can 


(6.3.56) 
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since  the  characteristic  matrix  A has  two  eigenvalues 

both  of  which  are  real  and  negative  being  - and  - — 

1 2 

(T  >0  and  T > 0) , asymptotically  as  t approaches  infinity, 
1 2 

P(t)  in  (6.3.56)  approaches  a stat  ifjtic.illy  stationary  value 
of  constant  P making  P(t)  = 0.  In  this  case,  6.3.56  becomes 


AP  + Pa'^  + =0  (6.3.57) 
Substituting  from  (6.3.10),  (6.3.51),  (6.3.11),  and  (6.3.55) 
into  (6.3.57),  we  obtain 


_ . 2 

— Ct 

2 1 


,0  Ip 


which  gives 


2 1 


22 


2a 


22 


In  view  of  (6.3.7),  (6.3.59)  becomes 


22  2t  T ( 

1 2 1 


T (T  + T ) 


(6.3.58) 


(6.3.59) 


(6.3.60) 


Since  P = E[X^(t)]  for  large  t and  y(t)  = T x (t) , we  have, 
22  2 1 2 

denoting  E(y^(t)]  for  large  t by  O^, 
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(6.3.61) 


g2  _ .j.2p 

y 122 


It  follows  from  (6.3.60)  that 


o'  = 


2t  (X  + X ) 
2 1 2 


(6 


As  a chec)c,  by  rewritinq  (6.3.62),  since  X 

related  to  the  break-frequencies  w and  w 

1 2 

w X = 1 by  (6.  3.  3)  , 

2 2 


and  X are 
2 


by  w X = 1 
1 1 


j2  _ £ 2 

y 2 w + w 
1 2 


f6 


We  can  clearly  see  that  0 increases  as  w increases, 

y 2 

should,  because  this  increases  the  band-width,  while  o 

decreases  as  w increases,  as  it  should,  because  this 
1 

decreases  the  band-width. 

Now  dividing  (6.3.35)  by  (6.3.62)  on  both  sides, 

denoting  E[W  ' (kT) ] for  large  kT  by  0 ' , we  obtain 

X2  '^1^2 


.3.62) 


and 


.3.63) 

as  it 

2 

y 

and 


.3.64) 
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From  (6.3.64)  with  T 


= 1,  T = 0.125  and  T = 0.01,  we  obtain 
2 


O = 0.4058  a 
WX2  ■ y 


(6.3.65) 


in  which  is  related  to  the  standard  deviation  of  experi- 
mentally obtained  data.  If  (6.3.14)  is  driven  by  W (kT) 

X 

generated  by  (6.3.48)  using  o = O above,  O = CO  and 

2 WX2  1 2 

the  values  of  p and  C as  given  by  (6.3.47)  for  T =1, 

1 

= 0.125  and  T = 0.01,  then  the  output  y[(k+l)T]  as  given 

by  (6.3.27)  would  have  the  same  band-pass  characteristics 

as  the  output  y(t)  from  (6.3.3).  Note  that  since  T = 1, 

1 

y[(k+l)T]  = X [(k+l)T]  in  (6.3.27). 

2 

The  validity  of  the  computer  routine  of  generating 

w (kT)  and  w (kT)  according  to  f (w  , w ) given  by  (6.3.48) 
X X X2  X X X2 

was  confirmed  by  recomputing  the  values  of  C and  p from 

the  computer-generated  w (XT)  and  w (kT) . The  validity  of 

XI  X2 

(6.3.64)  was  confirmed  by  recomputing  the  value  of  O = T a 

y 1 X2 

based  on  x [(k+l)T]  obtained  by  simulating  (6.3.25)  and 
2 

(6.3.26) ,  which  are  driven  by  w (kt)  and  w (kT) . The  O 

XI  X2  y 

in  (6.3.65)  is  related  to  the  experimental  data  obtained 

from  our  cat  (see  the  following  section) . 

As  a check  for  (6.3.62),  we  made  T =2  and  T = 1/16, 

1 2 

thus  increasing  the  band-width  of  the  filter  given  by  (6.3.3). 
obtained  from  the  simulation  of  (6.3.25),  (6.3.26),  and 

(6.3.27)  gave  a good  match  with  the  value  predicted  by 


(6.3.62),  which  shows  the  soundness  of  (6.3.62). 
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6.4  Statistical  analysis  in  the  model 


Referring  to  Figure  6.3,  r (kT)  and  r (kT)  are  discrete 

1 2 

noise  signals  coming  out  every  computation  interval  of  T 

seconds.  However,  r*(k'Tj  is  the  sampled  noise  on  the  WHEN 

1 

curve  at  k'T  when  a is  triggered  and  r*[(k'  + L)T]  is  the 

2 

sampled  noise  on  the  WHERE  curve  at  (k'  + L)T  where  LT  is 

equal  to  the  QP  duration.  L is  an  appropriate  integer. 

For  example,  L = 11  for  a QP  duration  of  0.11  sec  because 

T = 0.01  sec.  k'  is  the  value  of  k when  each  QP  is  triggered. 

Working  backwards  from  the  outputs  r*  and  r*  to  the  inputs 

1 2 

w (kT)  and  w (kT)  in  Figure  6.3,  we  want  to  determine  the 

X2  Z2 


standard  deviations 


0 /w  \ and  o/w  \ i 
{X2j  [ Z2  J 


in  terms  of  experi- 


mentally obtained  standard  deviations  corresponding  to 


a('r^*j  and 


Referring  to  Figure  6.3,  since  r = T x + T z and 

1 12  12 

since  x and  z are  uncorrelated  with  each  other,  txath  with 
2 2 

zero  means,  we  have 


“I  ('.) 

Because 

r ■'y(h) 

(6.4 

1 

(%)  = »'(h) 

we  have 

1 

(6.4 

* » 

1 

(h) 

(6.4, 
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Figure  6.3 

This  figure  shows  an  arrangement  by 


which  the  sampled  noise  r * on  the 

1 

WHEN  curve  and  the  sampled  noise  r * 

2 

on  the  WHERE  curve  are  created  through 

the  VOR  model  such  that  respective 

variance  of  r * and  r * as  well  as 
1 2 

the  correlation  between  r * and  r * 

I 2 

would  match  the  corresponding  values 
from  our  cat's  data. 


Figure  6.3  Noise  Generation 


It  follows  from  (6.3.64)  that 


(6.4.4) 


(6.4.5) 


Once  o' 


(v  \ and  0^(vi  ^ are  thus  determined,  x [ (k  + 1)T] 

\ X2j  \ Z2  J 2 


and  z [ (k  + 1)T]  are  generated  by  filters  represented  by 
2 

equation  I and  equation  II  in  Figure  6.3  from  (kT)  and 

w (kT) , lK)th  of  which  are  purely  random  Gaussicui  sequences  with 
Z2 

zero  mecms  but  with  different  variances  of  ^and 

respectively. 

Next,  we  want  to  determine  the  linear  correlation 

coefficient  R between  r * (k^T)  and  r * [(k^  + L)T] . k^T  is 
12  1 2 

the  time  at  the  beginning  of  a quick-phase,  while  (k'  + L)T 
is  the  time  at  the  end  of  the  same  quick-phase.  Thus,  LT 


is  equal  to  the  quick-phase  duration.  R is  obtained  from 

1 2 

(6.2.10)  by  replacing  t by  k'*T  and  t^  by  (k^  + L)T  and 
designating  E|^^r^*(t^)  j by  0^(r^*)  and  E|j^r^*(t^)  j by 


(r  *) . Thus, 
2 


12 


E|r  *{k  T)  r * (k  + L)T 
2 

0(r  * ) o(r  *) 

1 2 


] 


(6.4.6) 


Now,  for  the  sampled  noise  r * of  the  WHERE  curve 

2 


0(r  *)  o(r  ) = 1 0 (z  ) 
2 212 


(6.4.7) 
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However,  for  the  sampled  noise  r * of  the  WHEN  curve  we  can  say 

1 


0(r  *)  = 0(r  ) (6.4.8) 

1 1 

only  for  a relatively  high  slow-phase  eye  velocity.  As  the 
eye  velocity  becomes  slower,  the  eye  position  tends  to  touch 
only  the  nearest  part  of  the  noisy  WHEN  curve  because  of  its 
smaller  slope.  (This  is  illustrated  in  Figure  7.6.)  This 
would  malce  0(r  *)  smaller  than  O(r^)  as  the  eye  velocity 
becomes  slower.  This  effect,  which  we  might  call  "shading," 
is  not  too  significant  although  noticeable  in  our  model 
output  at  the  lowest  available  stop  input  head  velocity  of 
about  10  deg/sec  compared  with  the  highest  input  head 
velocity  of  about  60  deg/sec  for  our  cat.  This  point  is 
further  discussed  in  Chapter  7. 

With  this  understanding,  we  have 


fr  * (k'T)r  * ( (k'  + L)T  ] 

~ Efr  (kT)r  ((k  + L)t')1 

L 1 2 V J 

L 1 2V  )} 

= T^Ej^^x^(kT)  + z^(kT))  (()^  + L)t)J 

= T^E  ()iT)z^  ^()c  + L)Tjj  (6.4.9) 

because  x and  z are  uncorrelated.  It  follows  from  (6.4.6) 


that 


R 

12 


z (kT)z  (k  + L)T 
2 2 


T.) "(%) 


(6.4.10) 
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express  (6.4.10)  in  terms  of  known  quantities,  first 

we  define  P (kT)  and  c(kT,  (k  + L)T)  such  that 
z 


P (kT)  = E Z (kT)Z  (kT) 
z 


P P 

Zii  Z12 


P P 

Zi2  Z22 


(6.4.11) 


C^kT,  (k  + L)T^  = E jz  (kT)z'^((k  + L)T)j 


c c 
11  12 


c c 
2 1 22 


(6.4.12) 


Now,  (Bryson  and  Ho,  1975) , 


c(kT,  (k  + L)t)  = P^(kT)<I>'^  ((k  + L)T,  kT) 


= P (kT)4>  (LT) 
z 


(6.4.13) 


We  are  interested  in  the  c element  of  (6.4.12)  which  is, 

22 

using  (6.4.11),  (6.4.12),  and  (6.4.13), 


c (kT,  (k  + L)T]  = e|z  (kT)Z  (k  + L)t] 
22  [.2  2 J 


= (j)  (LT)p.,  + (})  (LT)p 

21  Zl2  22  Z22 


Equation  (6.4.14)  becomes,  using  (6  3.23)  and  (6.3.24)  with 


T replaced  by  LT, 
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LT 


(6.4.15) 


For  L = 0,  (6.4.15)  reduces  to  p , as  it  should,  because 

22  2 

c (kT,  kT)  = Elz  (kT)  z, (kT)  = P . For  a typical  quick- 

22  L 2 ^ J 222 

phase  duration  of  0.05  sec  to  0.2  sec,  L is  5 ~ 20  with 

T = 0.01  sec.  For  these  values  of  L with  T =1  and  T = 1/8, 

1 2 

the  first  term  inside  the  square  bracket  of  the  right  hand 
side  of  (6.4.15)  is  negligible  compared  with  the  second  term 
because  it  turns  out  that  Ip  | <<  jp  | for  T = 0.01. 


c [kT,  (k  +L)T] 
22 


LT 

T 


T 

1 


T e 
2 


(6.4.16) 


It  follows  from  (6.4.10)  with  (6.4.14)  and  (6.4.16)  that 


(6.4.17) 
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because  p 


1 


Equation  (6.4.17)  allows  us  to  take  the  physiologically 

measured  variables  a(r  ),  a(r  ) and  LT  and  the  estimated  frequency 

1 2 

limits  of  the  noise,  ^ and  — , and  allows  us  to  predict 

1 2 

the  correlation  coefficient  R . The  success  of  this  pre- 


1 2 

diction  is  discussed  in  the  next  chapter. 


6. 5 Summary  (Refer  to  Figure  6.3) 

The  noise  on  the  WHEN  curve  corresponds  to  r * and  the 

1 

noise  on  the  WHERE  curve  corresponds  to  r * . From  our  cat's 

2 

step  responses,  we  determine  variances  0^ (r  *)  and  0^  (r* ) as 

1 2 

well  as  the  linear  correlation  coefficients  R between  the 


1 2 

eye  position  at  the  beginning  and  that  at  the  end  of  each 

quick-phase.  From  the  experimentally  determined  values  of 

(r  *)  and  (r  *),  we  estimate  (r  ) and  0^ (r  ) by  (6.4.8) 
12  12 

and  (6.4.7).  Next,  we  determine  0^ (z  ) and  (x  ) by  (6.4.2) 

2 2 

and  (6.4.1).  After  that,  we  determine  (w  ) and  (w  ) by 

X2  22 


(6.4.4)  and  (6.4.5) 


w (kT)  and  w (kT)  are  generated  by  a computer  routine 

Xl  X2 

in  such  a way  as  to  meet  the  probability  density  function 


f (w  , w ) given  by  (6.3.48).  w (kT)  and  w (kT)  are 

Xi  X2  Zi  Z2 

generated  in  a similar  way  according  to  f(w^^,  (6.3.48) 

with  w and  w replaced  by  w and  w and  with  the 

Xi  X2  Zi  Z2 
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understanding  that  O = a (w  ) and  0 (w  ).  For  given 

1 zi  2 22 

0(r  *)  and  O (r  *)  corresponding  to  experimentally  obtained 
1 2 

standard  deviations,  the  computer  program  takes  care  of  the 

rest  of  the  procedure  which  yields  W^(kT)  and  W^(kT). 

Referring  to  Figure  6.3,  x [ (k  + 1)T]  and  z [ (k  + 1)T] 

2 2 

are  generated  respectively  through  equations  (I)  and  (II) 
from  W (kT)  and  W (kT) . x and  z have  the  Scune  frequency 

X 2 2 2 

band'pass  characteristics  as  the  filter  equation  (6.3.4). 

a(r  *)  and  o(r  *)  as  well  as  R (see  section  6.2)  computed 
1 2 12 

from  our  VOR  model  would  be  compared  with  those  obtained 

from  our  cat's  data.  R may  be  estimated  by  (6.4.17),  in 

1 2 

which  LT  may  be  interpreted  as  the  mean  quick-phase  dura- 
tion. In  the  computer  program,  R is  determined  by  (6.4.6) 

12 

using  running  averages.  The  R from  the  VOR  model  has  shown 

1 2 

a reasonable  match  with  our  cat's  data. 


148 


CHAPTER  VII 


RESULTS 


7.1  General  approach 


The  parameters  of  the  model  developed  thus  far  are 
adjusted  so  that  the  output  of  the  model  obtained  by  simu- 
lation would  approximately  match  the  responses  from  our  cat. 
This  is  done  in  two  stages.  Initially,  the  parameters  are 
adjusted  for  the  deterministic  model  without  noise  inputs. 
After  this,  the  parameters  of  the  noise  inputs  for  the 
stochastic  model  are  adjusted  to  meet  the  statistical 
characteristics  from  our  cat's  data. 

First,  the  quick-phase  frequency  from  the  determin- 
istic model  for  various  step  input  head  velocities  is  made 
to  match  approximately  with  the  mean  f^  of  our  cat's  VOR 
responses  by  adjusting  the  threshold  value  Q.  Next,  the 
deterministic  model  thus  obtained  is  further  modified  to 
obtain  reasonable  sinusoidal  responses  by  slightly  readjust- 
ing Q and  adjusting  the  time  constant  T of  the  lag  network 

L 

in  the  WHERE  function,  so  that  the  shapes  of  deteirministic 
sinusoidal  responses  at  test  frequencies  of  1.2,  0.25,  and 
0.05  Hz  are  reasonably  comparable  with  typical  sinusoidal 
responses  from  our  cat  at  these  frequencies. 

The  final  results  and  comparison  with  our  cat's  data 


are  shown  towards  the  end  of  this  chapter. 
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Many  of  the  system's  parameters  are  fixed  by  physiological 
measurements  and  are  not  adjustable.  Thus,  we  note  in  Fig- 
ure 3.1  that  t'  = 12  sec,  T = 25  sec,  T ~ 0.15  sec  fall  into 
c N e 

this  category.  In  the  high-gain  amplifier  of  the  neural  pulse 

generator,  the  slope  a = 16.2  and  the  intercept  point  b = 40 

(Doth  in  Figure  4.7)  are  also  based  on  neuroptiysiological 

observations  as  explained  previously  and,  in  any  event,  would 

not  influence  the  main  nystagmus  pattern.  This  leaves  the 

threshold  value  Q (Figure  5.11),  and  G and  T (Figure  5.8) 

L 

for  malting  relatively  small  amounts  of  adjustments.  However, 

G is  constrained  by  (5.4.1)  to  the  value  of  0.5.  Thus  we 
are  left  with  only  Q and  T as  available  for  minor  adjustments. 
Adjustment  of  these  two  parameters  is  discussed  in  the 
following  two  sections. 


7 . 2 The  adjustment  of  the  threshold  value  Q 

The  initial  estimate  of  the  threshold  value  Q is  adjusted 
so  that  the  quick-phase  frequency  (in  beats  per  sec)  f^  for 
the  model  step  response  and  the  quick-phase  frequency  per 
cycle  f^  for  the  model  sinusoidal  response  would  match 
reasonably  those  from  our  cat. 

After  some  trials  and  errors,  we  have  obtained  the  final 
profile  of  Q as  a function  of  the  SCC  output  signal  ({)  greater 
than  10  deg/sec  as  shown  in  Figure  5.9  and  explained  in 
section  5.4.  For  0 < |4>1  <10,  although  we  do  not  have  data 
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from  step  responses  to  determine  the  value  of  Q for  this 
region,  we  found  a constant  value  of  Q = 6 yields  acceptable 
sinusoidal  responses  as  explained  in  section  5.4. 

According  to  Table  5.1  in  Chapter  5,  the  mean  value 
of  the  steady  state  deterministic  WHERE  curve  is  about  15  deg. 
Yet,  in  our  model  we  saturated  the  deterministic  WHERE  curve 
at  12  deg  for  a canal  signal  greater  than  22  deg/sec.  This 
modification  is  concerned  with  the  location  of  the  WHEN 
curve  which  is  4^  Q. 

Although  for  the  step  inputs  the  deterministic  WHEN 
curve  appears  to  be  shifted  about  a few  degrees  beyond  0 deg 
(e.g.,  3 deg)  in  the  direction  of  head  motion,  for  the 
sinusoidal  responses  at  0.05  Hz  with  input  peak  head  velocity 
of  about  30  deg/sec  the  mean  initiation  point  of  quic)t -phases 
appears  to  be  on  the  0 degree  line.  Because  of  the  large 
amount  of  uncertainty  in  the  determination  of  and  the 
location  of  the  WHEN  curve  for  step  inputs  and  sine  waves, 
and  also  because  of  the  uncertainty  in  the  determination  of 
the  0 degree  line  and  the  smallness  of  the  discrepancy,  we 
arbitrarily  chose  as  12  deg. 

The  comparison  of  f^  (for  step  input)  and  f^  (for 
sinusoidal  input)  from  the  model  and  from  our  cat  is  shown 
in  Table  7.1  and  7.2,  respectively.  The  results  from  the 
model  are  obtained  with  the  parameter  values  as  specified 
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in  section  7,1  and  with  the  noises  (as  described  in  section  7.4) 
added.  Since  the  cat's  data  is  noisy,  the  noises  must  bo  added 
to  the  model  for  consistent  comparison  on  the  same  basis.  In 
Table  7.1  the  fact  that  f^  matches  fairly  well  between  the  cat 
and  the  model  output  indicates  that  our  original  estimation  of 
the  threshold  value  Q is  approximately  correct.  This  is  because 
f^  is  approximately  inversely  proportional  to  Q at  a constant 
slow-phase  velocity  since 


-2 


s 

Q 


(7.2.1) 


in  which  and  V are  respectively  QP  and  SP  mean  velocity. 


and  in  which  V >>  V 
' Q ' ' s 


Equation  7.2.1  shows  that  f^  is 


proportional  to  slow-phase  eye  velocity.  Because  of  (7.2.1) 
we  note  that  f^,  for  velocity  steps,  is  not  a very  demanding 
criteria.  All  one  needs  to  do  is  to  get  Q correct. 

Table  7.1 
Step  Responses 
Cat  Model 


6 

s 

deg/ sec 

-1 

sec 

sec 

8.9 

1.2 

1,3 

21.0 

1.7 

1.8 

30.8 

2.2 

2.6 

46,5 

2,5 

2.9 

slow-phase  eye  velocity  in  deg/sec 
quick-phase  frequency  per  sec 
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Table  7.2 


Sinusoidal  Response 

Cat  Model 


f 

Hz 

H 

M 

deg/sec 

f' 

Q 

cycle  ^ 

cycle 

1.2 

30 

1.1 

2.0 

0.25 

30 

6.6 

7.9 

0.05 

30 

23.0 

24.0 

f frequency  of  the  input  sinusoidal  head  velocity  in  Hz 

amplitude  of  the  input  sinusoidal  head  velocity  in  deg/sec 
f^  quick-phase  frequency  (in  beats)  per  cycle  of  the  input 
frequency 


In  Table  7.2,  f^  also  matches  fairly  well  except  at 
f = 1.2  Hz.  Although  the  model  makes  a quick-phase  at  each 
half-cycle  without  fail  at  f = 1.2  Hz,  the  cat  makes  on  the 
average  only  one  quick-phase  each  cycle.  The  reason  that 
might  account  for  this  will  be  discussed  in  a later  section 
on  hysteresis. 


7.3  The  adjustment  of  the  time  constant  T in  the  WHERE  function 

L 


The 

function 


transfer  function  of  the  lagged  portion  of  t!ie  WHERE 
is  given  by  (5.3.5)  as 


f (s)  _ G 
4>(s)  ~ ST^  + 1 


(7.3.1) 
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For  G = 0.5  and  T near  0.6  as  a rough  estimate,  the  break- 

Li 

frequency  of  (7.3.1)  is  about  0.27  Hz.  Thus,  for  an  input 
frequency  of  f =1.2  Hz,  (7.3.1)  reduces  to 


f (s)  ^ G_  1 
4)(s)  T s 

Li 


(7.3.2) 


which  acts  as  a pure  integrator.  For  this  reason,  at  f = 1.2  Hz 
the  phase  shift  is  close  to  90  degrees  and  its  change  due  to 
variation  of  T is  small  compared  with  the  amplitude  change  due 

Li 

to  the  same  cause. 

Figure  7.1  describes  how  the  component  of  the  deterministic 
WHERE  curve  C(t)  is  shaped  for  an  input  (})(t)  given  by 


i(t)  = A sin  27ift 


(7.3.3) 


with  f = 1.2.  In  Figure  7.1,  g(t)  is  the  step  portion  of  C(t) 
as  described  in  the  lower  portion  of  Figure  5.8  with  g^  = 1. 

Thus,  C(t)  = f(t)  + g(t). 

In  Figure  7.2,  the  interaction  of  the  eye  position  with 
C(t)  and  the  deterministic  WHEN  curve  D(t)  is  shown  [it  is 
drawn  parallel  to  C(t)  for  simplicity  although  in  the  model  it 
is  not  quite  (see  section  5.4)].  As  T is  increased,  the 

Li 

quick-phase  end-points  q , q , and  q would  move  closer  to 

1 2 3 

the  line  C(t)  = 0,  and  the  slow-phase  end-points  S , S , and  S 

1 2 3 

would  move  further  away  from  it.  If  we  denote  by  | 0 | the 

^ M 

maximum  magnitude  of  the  slow-phase  eye  position  change  in  the 
transition  region  (Figure  7.2),  this  would  make  |6  | bigger, 
thus  altering  the  shape  of  the  VOR  response.  Figure  7.3  shows  the 


effect  of  increasing  the  value  of  T by  a factor  of  about  two. 

L 
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Figure  7.1 

This  figure  describes  a sinusoidal 
see  signal  around  1.2  Hz,  and 

the  components  f(t)  and  g(t)  of  the 
corresponding  deterministic  WHERE 
curve  which  is  generated  by  <j)(t). 

Two  curves  are  shown  for  f{t).  The 
dashed  curve  shows  what  happens  if 
T is  doubled.  At  this  frequency 

Xj 

the  amplitude  of  f(t)  is  so  small, 
even  for  A equal  to  30  deg/sec,  that 
it  does  not  reach  the  clipping  level. 
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$(t) 


(The  time  constant  in  f(t)  for  the  dotted  curve  is 
about  twice  that  for  the  solid  curve) 


g(t) 
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Figure  7.2 


This  figure  shows  the  deterministic  WHERE 
curve  and  WHEN  curves  corresponding  to 
Figure  7.1  (with  f(t)  represented  by  the 
solid  line) , and  also  shows  the  interaction 
between  eye  position  and  the  deterministic 
WHERE  and  WHEN  curves. 


D(t) 


Figure  7,2  Interaction  Between  Eye  Position  and  WHERE  and 
and  WHEN  Curves  (corresponds  to  the  solid  curve 
of  f(t)  in  Figure  7.1) 
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Figure  7.3 

This  figure  shows  the  deterministic 
WHERE  curve  and  WHEN  curves  corresponding 
to  Figure  7.1  (with  f(t)  represented  by 
the  dotted  line) , and  also  shows  the 
interaction  between  eye  position  and  the 
deterministic  WHERE  and  WHEN  curves. 
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D(t) 


Figure  7.3  Interaction  Between  Eye  Position  and  WHERE  and 
WHEN  curves  (corresponds  to  dotted  curve  of 
f(t)  in  Figure  7.1) 
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0.6 


After  some  trials  and  errors,  we  have  found  that  T = 

gives  a good  approximation  to  the  shape  (ratio  — ^ ) of  the 

VOR  response  for  f = 1.2  Hz  with  A = 30  in  (7.3.3).  The 
sinusoidal  VOR  response  is  more  sensitive  to  T at  f = 1.2  Hz 
because  of  (7.3.2).  At  f = 0.25  Hz  and  f = 0.05  Hz,  | sT  | is 

i-i 

not  much  greater  than  1 in  (7.3.1)  and  the  nystagmus  is 

relatively  independent  of  T . In  the  step  responses,  the 

L 

determination  of  T was  based  on  rough  visual  inspection 
L 

of  curves  such  as  in  Figure  5.4  B-F.  This  was  an  uncertain 
way  to  estimate  T . For  this  reason,  the  value  of  T is 

Li  L 

adjusted  at  1.2  Hz  and  T =0.6  was  adopted  for  our  model. 

L 

7.4  The  variance  of  the  noise  input  to  the  model 

Referring  to  Figure  6.3,  we  can  see  that  the  standard 
deviations  0(r  *)  and  0(r  *)  depend  on  O (W  ) , a (W  ) , a (W  ), 

12  Xi  X2  zi 

and  O (W  ).  Because  of  (6.3.39),  we  note  that  O (W  ) = Ca (W  ) 
Z2  X2 

and  0(W  ) = CO (W  ) in  which  C is  given  by  (6.3.47).  For 

Zi  Z2 

given  (experimentally  obtained)  0(r  *)  and  0(r  *) , the  values 

1 2 

of  0(W  ) and  0(W  ) are  determined  by  means  of  (6.4.1)  to 

X2  Z2 

(6.4.5)  along  with  (6.3.64). 

According  to  the  data  from  our  cat's  responses  for  step 

input  head  velocities,  0(r  *)=  5 deg  and  0(r  *)  = 2.5  deg 

1 2 

for  |())1  22  deg/sec.  That  is,  the  noise  on  the  WHEN  curve 

is  about  twice  as  big  as  the  noise  on  the  WHERE  curve.  This 


161 


would  naturally  make  the  linear  correlation  coefficient 
between  each  slow-phase  magnitude  and  the  subsequent  quick- 
phase  magnitude  bigger  than  that  between  each  quick-phase 
magnitude  and  subsequent  slow-phase  magnitude.  Indeed,  the 
value  of  the  former  is  0.846,  while  that  of  the  latter  is 
0.326. 

However,  when  we  run  our  model  with  the  computed  values 

of  0(W  ) and  a (W  ) , we  have  observed  some  "backward"  quick- 

X2  Z2 

phases  in  the  same  direction  as  the  ongoing  slow-phases  in 
the  step  responses  from  our  model.  Around  (f)  = 22  deg/sec, 
we  have  observed  more  than  a few  of  them  for  each  computer 
run . 

The  backward  quick-phase  occurs  when  the  slow-phase  eye 
movement  crosses  the  WHEN  curve  in  the  region  in  which  the 
WHERE  curve  and  the  VvwEN  curve  overlap  as  shown  in  Figure  7.4. 
Note  that,  in  our  model,  whenever  the  eye  crosses  the  WHEN 
curve  it  seeks  the  WHERE  curve  regardless  of  where  that  curve 
may  be  located.  The  probability  of  the  two  curves  being 
overlapped  increases  as  (<}>|  falls  below  22  deg/sec  because 
of  the  decreasing  value  of  Q (see  Figure  5.9). 
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Figure  7.4 

This  figure  shows  the  situation  in 
which  the  backward  quick-phase  occurs 
when  the  slow-phase  eye  movement 
crosses  the  WHEN  curve  in  the  region 
in  which  the  WHERE  curve  and  the  WHEN 
curve  overlap. 
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Because  we  have  not  observed  any  backward  quick-phases 

in  our  cat's  step  responses,  with  several  trials  for  head 

velocities  equal  to  or  over  10  deg/sec,  we  have  decided  to 

eliminate  them  by  controlling  the  noise  generation  at  its 

source  level,  rather  than  at  its  output  level.  That  is,  we 

have  decided  to  eliminate  them  by  limiting  the  magnitude  of 

the  four  one-dimensional  independent  Gaussian  sequences 

(with  zero  mean  and  variance  1)  from  which  W and  W are 

X z 

generated,  until  no  backward  quick-phases  are  observed  in 
the  model  output  for  step  input  head  velocities  above 
22  deg/sec.  We  have  observed  that  by  limiting  the  magnitudes 
of  the  Gaussian  sequence  sources  to  1.65  times  their 
standard  deviations  (which  are  all  equal  to  one) , we  could 
eliminate  the  backward  quick-phase  for  | (])  | ^22.  For  a Gaussian 
distribution,  with  this  limitation,  90%  of  all  possible 
numbers  are  still  included.  Every  time  this  limit  is 
exceeded,  the  computer  brings  out  another  number  from  the 
Gaussian  distribution  until  the  requirement  is  met.  This 
results  in  another  distribution  which  is  approximately 
Gaussian  with  slightly  smaller  variance  than  the  original ^ 
one.  This  would  somewhat  decrease  the  predicted  output 
variances.  These  variances  are  subject  to  a further  decrease 
due  to  another  problem  which  we  call  the  "shading  effect"  to 
be  discussed  later  in  this  section  (see  Table  7.3).  Due  to 


1/65 


sa- 


the  SCC's  12-socond  time  constant,  the  number  of  data  points 
in  each  step  response  from  our  cat  is  limited  to  not  more 
than  about  a dozen  points.  In  addition,  the  platform  on 
which  our  cat  is  rotated  was  not  servo-controlled  and,  there- 
fore, could  not  maintain  a constant  velocity  within  an  indi- 
vidual run  or  several  runs  intended  for  the  same  nominal  head 
velocity.  Under  these  situations,  statistics  have  large 
variabilities  and  we  have  not  insisted  too  much  on  an  exact 
numerical  match  between  model  and  experimental  statistics. 

The  probability  of  the  two  curves  being  overlapped 
increases  as  | (J)  | falls  below  22  deg/sec  because  of  the 
decreasing  threshold  value  Q,  as  mentioned  before.  This 

necessitates  decreasing  the  values  of  a(W  ) and  a(W  ) to 

X2  Z2 

prevent  backward  quick-phases  at  lower  velocities.  In  addi- 
tion, for  values  of  ()>  near  zero,  it  was  necessary  to  decrease 

the  values  of  o (W  ) and  O (W  ) for  the  following  additional 
X2  Z2 

reason.  For  a sinusoidal  input  head  velocity  of  0.05  Hz 

(with  an  amplitude  of  about  30  deg/sec) , the  transition 

region  within  which  the  slow-phase  of  the  VOR  response  goes 

through  zero  velocity  (t  in  Figure  7.21)  is  about  2 sec- 

TR 

onds.  Too  much  noise  near  zero  velocity  interrupts  this 

transition  region  with  quick-phases,  making  the  resulting 

duration  t„,  much  too  short. 

TR 


166 


One  problem  in  this  regard  is  that  cats  make  saccades 


in  the  dark.  Perhaps  the  rapid  eye  movements  that  occur 

for  very  small  head  velocities,  whether  in  the  step  responses 

or  in  the  transition  region  of  low  frequency  sinusoidal 

responses,  are  really  saccades  and  not  quick-phases.  However, 

for  modeling  purposes  only,  we  assume  this  model,  which  has 

no  saccadic  system,  should  not  make  quick-phases  for  very 

small  head  velocities.  In  practice,  when  the  head  is  still, 

the  model  should  produce  no  rapid  eye  movements.  To  do  this, 

it  is  necessary  to  greatly  reduce  the  noise  when  |(})|  =0. 

One  way  to  meet  these  requirements  is  to  reduce  0(x  ) 

2 

and  a (z  ) of  X and  z (see  Figure  6.3)  for  low  velocities. 

2 2 2 

It  is  often  the  case  in  neurophysiology  that  a reduced 

intensity  of  a neural  signal  is  accompanied  by  a reduction 

in  its  noise.  In  the  model,  0(x  ) and  0(z  ) for  I <})  < 22 

2 2 

are  made  to  be  proportional  to  | 4> | ^ as  explained  below.  We 

have  assumed  that  a(x  ) and  O (z  ) should  be  monotonically 

2 2 

decreasing  as  |(J)|  approaches  zero,  and  should  become  zero 

at  l4i|  =0,  although  it  does  not  have  to  become  zero  at 

||})1=0.  At  1(1)1  =22,  we  require  that  0(z  ) = — o(r  *)  = 2.5 

2 X 1 2 

for  X =1  [sec  (6.4.2)]  and  0(x  ) = — (r  *)  - (r  *)  = 

1 2 T 1 2 

f ^ 

V 5^  - 2.5^  for  X =1  (see  (6.4.3)]  to  meet  the  boundary 
1 

condition  at  | <1)  | = 22.  Note  that  a(x  ) and  0(z  ),  for 

2 2 

K|  > 22,  are  constant  and  have  the  same  values  as  at  |({)1  = 22. 
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The  profiles  of  o(x  ) and  o(z  ) are  shown  in  Figure  7.5,  in 
2 2 

which  the  values  of  K and  K arc  determined  by  the  boundary 

1 2 

conditions  at  | (})  | -22.  We  have  used  a cubic  function  of 

I (J)  I because  a linear  or  square  function  does  not  reduce  the 

noise  level  fast  enough  to  meet  our  requirement  for  the 

duration  of  the  transition  region  t„  at  f = 0.05  Hz.  This 

cubic  noise  condition  also  eliminates  all  backward  quick- 

phases  for  I (|)  I < 22.  The  profiles  of  the  variances  of  0(x  ) 

2 

and  0(z  ) in  Figure  7.5  are  created  in  the  model  by  similarly 
2 

shaping  the  variances  of  the  noise  sources  a and  o 

wx  wz 

That  is,  in  the  model  O = 0. 40587 , O (x. ) and  O 

' WZ2  1 2 WZ2 

= 0. 4058T  jO  (Zj, ) based  on  (6.4.4)  and  (6.4.5).  According  to 
(6.3.64)  and  (6.3.65),  [Fd^,  1^,  T)  ] is  0.4058  for  ij  = 1, 

= 0.125  and  T = 0.01. 

We  recognize  that  these  limitations  on  noise  would  cause 

decreases  in  0(r  *)  and  0(r  *)  computed  from  the  model  output 
1 2 

at  low  velocities  in  addition  to  decreases  in  0(r  *)  and 

1 

o(r  *)  caused  by  clipping  off  the  peaks  of  the  one-dimensional 
2 

Gaussian  sources  for  all  velocities.  There  is  still  another 

reason  that  0(r  *)  would  be  further  decreased.  We  call  it 
1 

the  "shading  effect"  (Figure  7.6) . As  slow-phase  velocity 
becomes  smaller,  the  eye  position  tends  to  touch  the  WHEN 
curve  more  often  near  its  lower  peaks  as  shown  in  Figure  7.6. 
This  would  naturally  make  the  sampled  WHEN  noise  less 


Figure  7.5 

This  figure  shows  the  profiles  of  the 

standard  deviations  o(x  ) and  0(z  ) of 

2 2 

noises  x (t)  and  z (t)  used  in  our 
2 2 

stochastic  VOR  model. 


I 


I 
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Figure  7.6 

This  figure  shows  what  we  call  "shading  effect" 
in  which  as  slow-phase  velocity  becomes  smaller, 
the  eye  position  tends  to  touch  the  WHEN  curve 
more  often  near  its  lower  peaks. 
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igure 


fluctuating,  thus  making  oCr^*)  smaller.  However,  as  the  eye  ve- 
locity becomes  faster,  a{r^*)  would  gradually  increase  approaching 
the  standard  deviation  0{r^)  of  the  noise  itself.  Thus,  calcu- 
lated values  of  oCr^)  overestimate  the  actual  values  of  0(rj*) 
produced  by  the  model.  The  "shading  effect"  is  only  a problem  for 
the  sampled  WHEN  curve,  and  not  the  sampled  WHERE  curve. 

Before  closing  this  section,  we  want  to  make  a comment  on  the 
refractory  period  of  50  msec  for  quick-phases  instead  of  200  msec. 
Originally  we  started  with  a 200  msec  refractory  period  in  our 
model  because  we  assumed  at  first  that  the  refractory  period  for 
quick-phases  must  be  similar  to  that  for  saccades  which  is  recog- 
nized as  about  200  msec.  However,  with  a 200  msec  refractory 

period  the  linear  correlation  coefficient  R (6.4.6)  from  the 

1 2 

model  was  much  too  low  compared  with  those  obtained  from  the  cat's 
data.  The  reason  was  that  when  the  noise  demanded  a short  duration 
SP  (e.g.,  6 deg  in  size  at  a velocity  of  40  deg/sec  which  takes 
only  0.15  sec)  the  refractory  period  would  not  allow  it.  This 
causes  A9  to  cross  the  WHEN  curve  and  the  quick-phase  occurred 
later  at  an  eye  position  that  was  not  on  the  WHEN  curve,  and  which 
was  consequently  uncorrelated  with  the  WHERE  curve.  With  this 
in  mind,  we  have  carefully  re-examined  the  cat's  analogue  data, 
and  found  that  the  refractory  period  for  our  cat  is  actually  only 
50  msec  instead  of  200  msec,  contrary  to  our  previous  assumption. 
After  we  have  changed  the  refractory  period  to  the  new  value  of 
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50  msec,  the  values  of  R increased,  giving  acceptable 

1 2 

results  compared  with  that  from  our  cat  as  shown  in  Table  7.3. 

7 . 5 Results 

The  final  results  are  obtained  with  parameter  values  and 
noise  models  determined  as  previously  described.  In  what 
follows,  the  results  obtained  from  the  VOR  model  output  for 
various  step  input  head  velocities  are  shown  in  subsection  (A) , 
and  the  results  for  sinusoidal  input  head  velocities  are  shown 
in  subsection  (B) . The  discussion  of  the  results  and  other 
comments  are  given  in  section  7.6  and  section  7.7. 

(A)  Results  - step  input  head  velocity 

The  typical  step  responses  from  our  cat  are  shown  in 
Figures  7.7,  7.8,  7.9,  and  7.10.  The  corresponding  step 
responses  from  the  stochastic  (with  noise)  VOR  model  for 
similar  input  head  velocities  are  shown  in  Figures  7.11,  7.12, 

7.13,  and  7.14. 

The  numerical  results  from  our  cat  and  from  our  VOR  model 
for  various  step  input  head  velocities  are  shown  in  Table  7.3, 
in  which  various  symbols  have  the  following  meanings: 

0^:  slow-phase  eye  velocity  in  deg/sec 

0^:  standard  deviation  of  the  scunpled  noisy  WHEN  signal  in  deg 

a : Standard  deviation  of  the  sampled  noisy  WHERE  signal  in  deg 
WR 

R : linear  correlation  coefficient  between  the  WHEN  and 

1 2 

WHERE  signal  [same  as  R of  (6.4.6)] 

12 
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Table  7.3.  Step  Responses 


e 

s 

(deg/sec) 

a 

WN 

cat 

(deg) 

a 

WN 

model 

(deg) 

o 

WR 

cat 

(deg) 

a 

WR 

model 

(deg) 

cat 

model 

8.9 

3.8 

0.5 

2.5 

0.4 

0.46 

0.67 

21.0 

5.0 

3.0 

2.7 

1.6 

0.34 

0.32 

30.8 

4.7 

3.4 

2.1 

1.8 

0.  30 

0.30 

46.5 

4.9 

3.5 

2.8 

1.7 

0.27 

0.16 

I 
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Figure  7.7,  7.8,  7.9,  and  7.10 
These  figures  show  the  responses  from  our 
cat  for  various  step  input  head  velocities. 


176 


A042  342 


JNCLASSIFIEO 

3<r3 

1^42342 


NAVAL  SURFACE  NEaPONS  CENTER  DAHLSREN  LAB  VA 
mathematical  model  of  the  VESTXBULOOCULAR  reflex. (U) 
JUN  77  K S CHUN 

NSmC/OL-TR-3669 


Figure  7.9  Step  Response  from  Our  Cat  (with  mean  slow-phase  velocity  of  30.8  deg/sec) 


Figure  7.10  Step 


Figures  7.11,  7.12,  17.13,  cmd  7.14 
These  figures  show  responses  from  our  stochastic 
VOR  model  for  various  step  input  head  velocities 
similar  to  the  input  velocities  applied  to  our 
cat. 
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slow-phase  velocity  of  8.9  deg/sec) 


Figure  7.12  Step  Response  from  Our  Stochastic  VOR  Model 

(with  mean  slow-phase  velocity  of  21.0  deg/sec) 


Figure  7.13  Step  Response  from  Our  Stochastic  VOR  Model 

(with  mean  slow-phase  velocity  of  30.8  deg/sec) 


Figiire  7.14  Step  Response  from  Our  Stochastic  VOR  Model 

(v,ith  mean  slow-phase  velocity  of  46.5  deg/sec) 


O and  0 „ from  the  model  correspond  to  o(r  *)  and  0(r  *)  of 
WN  WR  12 

r * and  r * in  Figure  6.3.  More  specif ically , O and  0 
1 2 WN  WR 

from  the  model  are  computed  by  the  comimter  (using  running 
means)  based  on  the  efference  copy  of  the  eye  position  9^ 
(Figure  6.2)  at  the  beginning  and  at  the  end  of  quiclc-phases, 
respectively. 

(B)  Results  - sinusoidal  input  head  velocity 

The  typical  sinusoidal  responses  from  our  cat  are  given 
in  Figures  7.15  (for  1.2  Hz),  7.16  (for  0.25  Hz),  and  7.17 
(for  0.05  Hz)  all  with  the  pea)<  input  head  velocity  of 
30  deg/sec.  The  corresponding  sinusoidal  responses  from 
our  stochastic  VOR  model  are  given  in  Figures  7.18  (for 
1.2  Hz),  7.19  (for  0.25  Hz),  and  7.20  (for  0.05  Hz)  all  with 
the  same  pea)t  input  velocity  of  30  deg/sec.  The  numerical 
results  of  the  sinusoidal  responses  from  both  the  cat  and 
the  model  are  tabulated  in  Table  7.4.  Some  of  the  symbols 
used  in  Table  7.4  are  depicted  in  Figure  7.21.  Referring 
to  Figure  7.21,  is  not  generally  equal  to 

may  happen  that  > |0^1  for  the  cat  while 

|9„l  1 0^1  for  the  model  and  vice  versa.  We  have  defined 

1 0^1  by  1 0^1  = — ^~2 ^ simply  to  get  an  estimate  of 

Q when  head  velocity  is  near  zero. 
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figures  7.15,  7.16,  and  1.11 
These  figures  show  the  sinusoidal  responses 
from  our  cat  with  input  frequencies  of 
1.2  Hz,  0.25  Hz,  and  Hz,  and  with 

amplitude  of  30  deg/sec. 
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Lgure  7.16  Sinusoidal  Response  from  Our  Cat  (with 
of  0.25  Hz  and  amplitude  of  30  deg/sec 


Figure  7.17  Sinusoidal  Response  from  Our  Cat  (with  input  frequency 
of  0.05  Hz  and  amplitude  of  30  deg/sec) 


Figures  7.18,  7.19,  and  7,20 
These  figures  show  the  sinusoidal  responses 
from  our  stochastic  VOR  model  with  input 
frequencies  of  1.2  Hz,  0.25  Hz,  and  0.05  Hz, 
and  with  amplitude  of  30  deg/sec. 
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Figure  7.18  Sinusoidal  Response  from  Our  Stochastic  VOR  Model  (with  input 
frequency  of  1.2  Hz  and  amplitude  of  30  deg/sec) 


Figure  7.19  Sinusoidal  Response  from  Our  Stochastic  VOR  Model  (with  input 
frequency  of  0.25  Hz  and  amplitude  of  30  deg/sec) 


Figure  7.20  Sinusoidal  Response  from  Our  Stochastic  VOR  Model  (with  input 
frequency  of  0.05  Hz  and  amplitude  of  30  deg/sec) 


In  Table  7.4,  the  various  symbols  have  the  following 


meanings: 


f : 


t 


TR' 


frequency  of  Lhe  input  sinusoidal  head  velocity 
in  Hz 

amplitude  of  the  input  sinusoidal  head  velocity 
in  deg/sec 

quick-phase  frequency  (in  beats)  per  cycle  of  the 
input  frequency  (not  per  sec) 

time  duration  in  sec  of  the  transition  region  in 
which  the  slow-phase  eye  velocity  goes  through 
zero  without  an  interruption  (Figure  7.21) 
the  maximum  magnitude  of  the  eye  position  dis- 
placement (trough  to  peak)  in  the  transition 
region  caused  by  the  slow-phase  in  deg  (Fig- 
ure 7.21) 

average  magnitude  of  the  two  quick-phase  dis- 
placements immediately  before  and  immediately 
after  the  transition  region  in  deg  (Figure  7.21) 


7.6  Discussion  of  results 

First  we  consider  the  step  responses  (refer  to  Table  7.3). 

For  0 = 8.9,  the  differences  in  a„,  and  those  in  0,  between 

s WN  WR 

the  cat  and  the  model  are  significant.  The  low  values  of 

O, and  from  the  model  at  6 =8.9  are  caused  by  the 

WN  WR  s 


196 


Figure  7.21 

This  figure  depicts  schematically  t » |6  ] 
and  1 0^1  used  in  Table  7.4.  These  are 
defined  in  the  text  in  connection  with 
Table  7.4.  Two  frequencies  are  shown  one 
low,  one  high,  but  the  symbols  are  the  same 
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particular  noise  models  we  have  adopted  (see  Figure  7.5). 

As  discussed  in  section  7.4,  if  the  noise  is  not  decreased  in 

this  region,  many  backward  quick-phases  would  occur.  How  the 

cat  manages  to  have  a of  3.8  deg  and  a 0^^  of  2.5  deg 

without  making  backward  quick-phases  is  a problem  unresolved 

in  this  study.  For  6^  ^22,  both  0^  and  for  the  model 

are  consistently  lower  than  those  for  the  cat.  This  is 

partially  due  to  the  limitations  of  the  magnitudes  of  the 

independent  Gaussian  numbers  (with  zero  mean  and  variance 

one)  to  the  value  of  1.65  (see  section  7.4).  For  example, 

in  one  simulation  run  with  6 = 46.5  with  this  noise  limita- 

s 

tion,  a.,,  was  3.4  deg,  while  without  the  noise  limitation 

WN 

it  was  3.9  deg.  With  the  noise  limitation  <7  was  1.7, 

WR 

while  without  the  noise  limitation  it  was  2.6. 

The  discrepancies  in  and  0,„,  for  6 > 22  deg/sec 

WR  WN 

could  probably  be  eliminated  by;  successive  readjustment 

of  the  noise  bandwidth  T and  T ; the  noise  intensities 

1 2 

O ,0  ; the  shapes  of  0(l4)|)  (Figure  7.5)  and  Q(|<t>|) 

WX  2 VfZ  2 

(Figure  5.9);  and  using  skewed  rather  than  Gaussian  noise 

distributions.  This  was  not  done  because  the  limited  amount 

of  physiological  data  keeps  the  experimental  values  of 

and  a from  being  much  more  than  rough  estimates  of  the 
WN 

true  values. 

As  for  R , we  consider  the  match  to  be  reasonable 
12 

considering  limited  samples  from  both  the  computer  and  the 
cat. 
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Next,  we  consider  the  sinusoidal  responses  (Table  7.4). 

Overall,  in  view  of  limited  samples  from  our  cat  (only  one  or  two 

runs  at  a particular  frequency  and  amplitude) , we  consider  that 

the  responses  from  the  model  compare  reasonably  with  those  from 

the  cat.  For  f = 1.2  Hz,  the  match  between  the  cat  and  the  model 

is  reasonable  except  for  f^.  A possible  explanation  for  this 

discrepancy  is  hysteresis  as  discussed  in  section  7.7. 

For  f = 0.25  Hz,  the  mismatch  both  in  t and  l9  | is  siq- 

TR  s'. 

nificant.  This  is  partially  the  result  of  our  modeling  approach. 

We  have  adjusted  the  model  parameters  (Q  and  T ) and  noises  to 

L 

meet  the  requirements  for  f = 1.2  Hz  and  f = 0.05  Hz,  the  two 

frequency  extremes,  and  hoped  for  good  luc)<.  at  f = 0.25  Hz.  We 

did  not  get  it.  However,  the  hysteresis  effect  may  again  be  to 

blame.  This  will  be  described  in  the  next  section,  but  briefly, 

the  values  of  t and  |0  I are  sensitive  to  H . We  have  observed 

TR  ' s ' M 

M 

in  the  model  that,  for  example,  when  H is  decreased  from  30  to 

M 

20  deg/sec,  the  value  of  t and  [O  | jumps  discontinuously 

TR  s ^ 

becuse  one  quiclc-phase  may  drop  out  of  each  half  cycle.  This  can 
increase  these  parameters  by  as  much  as  a factor  of  2. 

Referring  to  Table  7.4,  we  regard  the  comparison  in  |9q1 
between  the  cat  and  model  as  rough  justification  for  the  dependence 
of  Q on  j)  shown  in  Figure  5.9,  especially  at  low  head  velocities 
which  occur  near  the  zero  crossings  at  0.05  Hz. 


200 


7.7  Comments 


First,  we  will  discuss  the  so-called  "hysteresis  effect"  (Fig- 
ures 7.22  and  7.23).  For  f = 1.2  Hz,  for  example,  as  the  ampli- 
tude H of  the  peak  head  velocity  gradually  increases,  suddenly 
M 

a quick-phase  is  initiated  at  a certain  amplitude  h'  and  then 

M 

occurs  in  all  subsequent  cycles.  However,  as  the  amplitude 

gradually  decreases,  the  quick-phases  suddenly  disappear  at 

another  amplitude,  h"  . It  is  always  the  case  that  h'  > h"  , 

M MM 

similar  to  a magnetic  hysteresis  loop.  In  our  model,  (at 

• • 

1.2  Hz)  h'’  occurs  at  14.75  deg/sec  while  n'/  occurs  at  7.77  deg/sec 
M M 

for  the  noise-free  model  as  shown  in  Figure  7.22.  For  the 

stochastic  model,  the  values  of  h'  and  ii"  fluctuate  from  trial 

M M 

to  trial,  h'  fluctuates  about  + 3 deg/sec  around  a mean  value 
M ~ 

of  14.75  deg/sec,  while  H^'  fluctuates  about  + 1 deg/sec  around 

a mean  value  of  7.96  deg/sec.  For  our  cat,  h'  appears  to  be 

M 

about  10  deg/sec,  while  appears  to  be  about  7 deg/sec. 

When  the  experimental  cat  data  was  obtained,  hysteresis  was  not 

explicitly  looked  at  so  the  data  is  limited  and  H and  H are 

M M 

very  rough.  The  match  between  model  and  cat  is  thus  maybe 

better  thcin  one  might  expect. 

Because  of  this  phenomenon,  if  or  f'  were 

plotted  against  H , they  would  change  by  a jump  each  time 
M 

a new  quick-phase  was  added  to  the  output.  This  means  that 

these  measures  of  performance  are  sensitive  to  H and  could 

M 
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Figure  7.22 

This  figure  depicts  what  we  call  the 
"hvsteresis  effect".  As  the  magnitude  of  the 
input  head  velocity  gradually  increases, 
suddenly  a quick -phase  is  initiated  and  then 
occurs  in  all  subsequent  cycles.  As 
the  magnitude  gradually  decreases,  the  quick- 
phases  suddenly  disappear  at  another 
magnitude.  It  is  always  the  case  that  the 
magnitude  of  the  former  is  greater  than  the 
magnitude  of  the  latter. 


Figure  7.23 

This  figure  depicts  the  effects  of  the 

initial  eye  position  on  t and  | 9 | 

TR  ^ H 

(defined  in  Figure  7.21).  Sinusoidal 
head  velocity  is  at  0.25  Hz.  Initial 
position  can  modify  the  quick-phase 
frequency  per  cycle. 
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I I I I 


on 


explain  in  (Table  7.4)  the  differences  between  cat  and  model 

in  f'  at  1.2  Hz,  and  in  t and  |0  I at  0.25  Hz. 

Q TR  ' s' 

M 

As  usual  with  a hysteresis  phenomenon,  the  behavior  depends 

on  ititial  position.  If  a sinusoidal  stimulus  is  applied  with 

H such  that  there  are  two  possible  values  for  the  number  of  QPs, 
M 

because  of  hysteresis,  which  solution  will  result  depends  on 

i.nitial  eye  position.  Figure  7.23  shows  an  example  at  0.25  Hz. 

The  two  solutions  have  quite  different  values  for  t and  |9  | . 

® M 

If  one  unfortunately  chooses  a value  of  H.,,  for  cat  or  model, 

M 

that  is  near  such  a jump  discontinuity  point,  the  measures  we 

have  chosen  will  fluctuate  considerably  and  create  more  apparent 

discrepancy  between  cat  and  model  that  may  actually  exist. 

I Peak  head  velocity  H was  not  systematically  studied 

M 

in  this  investigation.  One  conclusion  from  our  model  is 

that,  for  future  investigations,  the  parameters  in  Table  7.4 

are  not  so  important  by  themselves  at  any  one  of  f and 

as  the  way  in'  which  they  vary  with  H^. 

Next,  we  will  make  a comment  on  the  sensitivity  of  R , 

1 2 

approximated  by  (6.4.17),  to  changes  in  quick-phase  duration 

and  changes  in  the  noise  band  edge  frequencies  — and  — . 

T j 

' These  values  of  R are  given  in  Table  7.5.  Judging  from 

1 2 

p,  1 our  cat's  data,  most  realistic  values  of  L lie  between  5 and 

20  with  10  the  nominal  value.  Referring  to  the  trapezoid- 
shaped  filter  in  Figure  7.24,  if  the  lower  frequency,  Wj , is 
decreased  by  a factor  of  two,  thus  making  Tj  =2,  there  is  very 
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little  change  in  R . This  is  not  surprising.  The  correlation 
1 2 

should  depend  more  on  high  rather  than  low  frequencies.  In  fact, 

if  the  higher  break-frequency,  w , is  increased  by  a factor  of 

2 

two,  thus  making  T = the  value  of  R for  L = 10  drops 
2 16  12 

sharply.  This  means  that  the  choice  of  T is  fairly  critical. 

2 

Note  that  the  values  of  R in  Table  7.5  (e.g.,  the  value 

1 2 

0.19  in  the  nominal  case  for  L = 10)  are  determined  from  (6.4.17) 

with  0(r  ) = 5.0  and  a(r.)  = 2.5,  while  the  values  of  R (for 
* 12 

the  model)  in  Table  7.3  (typically  0.3)  are  computed  based  on 
the  simulation  output  of  our  stochastic  model  with  its  associated 
effects  and  its  noise  limitations.  Consequently,  the  main  value 
of  6.4.17  is  to  provide  a rough  quantitative  estimate  of  experi- 
mental data  and,  more  important,  to  show  how  R should  depend 
. 12 

on  parameters  like  T , T , T,  and  L. 

1 2 
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Table  7.5 


Nominal  Bandwidth  High  Frequencies 


(for  c 

1 

L 

= 1,  7 = 1/8) 

2 

R 

1 2 

(for  T = 1, 
1 

L 

T = 1/16) 
2 

R 

1 2 

0 

0.50 

0 

0.50 

1 

0.46 

1 

0.42 

5 

0.31 

5 

0.21 

10 

0.19 

10 

0.08 

15 

0.11 

15 

0.02 

20 

0.06 

20 

~0 

LOW  Frequencies 

Both 

(for  X = 2,  X = 
1 2 

1/8) 

(for  X = 2,  X 
1 2 

= 1/16) 

L 

R 

L 

R 

12 

1 2 

0 

0.50 

0 

0.6 

1 

0.46 

1 

0.42 

5 

0.32 

5 

0.22 

10 

0.21 

10 

0.09 

15 

0.13 

15 

0.03 

20 

0.08 

20 

-0 
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Figure  7.24 

This  figure  depicts  the  effect  of  the 

band-pass  characteristic  when  T =1 

is  increased  by  a factor  of  two  and 

T = ^ is  decreased  by  a factor  of 
2 o 

two. 
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CHAPTER  VIII 


DISCUSSIONS  AND  COMMENTS 

8. 1 Modification  of  the  VOR  model  for  saccades 

In  evolution,  when  the  fovea  developed  along  with  the 
need  for  saccades,  the  saccadic  system  appears  to  have  utilized 
the  quick-phase  circuits  because  the  two  movements  are  identical 
and  share  the  same  neural  circuits.  For  this  reason, -a.  good 
VOR  model  should  be  able  to  accommodate  the  signals  for 
saccadic  eye  movements  without  too  much  modification.  The 
modifications  in  our  VOR  model  to  do  this  are  shown  in 
Figure  8.1. 

We  want  our  model  to  be  able  to  make  a saccadic  eye 
movement  when  perceptual  systems  have  selected  a visual  target. 
Referring  to  Figure  8.1,  the  selection  of  a visual  target  eind 
the  decision  to  make  a saccade  requires  moving  switch  S from 

3 

s to  s , and  switch  S from  s to  s . These  arrangements 

31  32  4 41  42 

are  made  for  the  following  reasons.  For  a saccade,  the  vesti- 
bular  WHERE  signal  must  be  replaced  by  a visual  WHERE  signal 
which  corresponds  to  the  position  of  the  target  in  space. 

Since  the  decision  to  acquire  a visual  target  dictates  an 
immediate  rapid  eye  movement  to  the  target  position,  provided 
it  meets  the  refractory  period  of  200  milliseconds  (see 
section  5.1),  the  WHEN  function  must  also  be  bypassed  by 
setting  switch  S to  s from  s . The  difference  in  the 
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Figure  8.1 

This  figure  shows  the  modified  VOR 
model  with  arrangement  for  saccades. 
The  switch-positions  for  the  VOR  mode 
are  S at  s and  S at  s , while 

3 3 1 “♦  **  1 

those  for  the  saccadic  mode  are  S 

3 

at  s and  S at  s . 

3 2 “*  2 


•W. 
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Figure  8.1  Our  Modified  VOR  Model  with  Arramgement  for  Saccades 


refractory  periods  between  the  saccades  and  quick-phases  may  be 
attributable  to  the  signal  processing  time  required  by  the  visual 
system  which  separates  saccades  by  200  msec.  This  might  be  bypassed 
by  the  vestibular  system  for  large  fast  head  movements  when  the  de- 
mand for  a second  quick-phase  overrides  the  value  of  visual  analysis 
at  the  end  of  the  first. 

Essentially,  in  Figure  8.1,  the  basic  motor  machinery  is  shared 
by  saccades  and  quick -phases  but  the  control  signals,  WHEN  and 
WHERE,  are  just  taken  over  by  one  system  or  the  other.  In  foveate 
animals,  in  the  light,  the  saccadic  system  is  usually  in  control. 
During  rotation  in  the  dark,  the  VOR  system  is  usually  in  control. 
But  there  are  many  intermediate  situations  in  which  they  might  com- 
pete. This  means  that  the  switches  are  an  oversimplification  and 
they  really  represent  some  sort  of  nonlinear,  multi-input,  decision 
process.  However,  there  have  been  no  studies  of  how  saccades  and 
quick-phases  interact  with  each  other  so  it  is  not  possible  to 
improve  on  the  schevie  in  Figure  8.1. 

A new  AND-gate  II  assures  a 200-millisecond  refractory  period 
for  saccades,  instead  of  the  50-millisecond  refractory  period  for 
quick-phases  used  in  AND-gate  I.  Again,  this  artificial  separation 
into  two  refractory  periods  points  out  lack  of  knowledge  about 
saccade  and  quick-phase  interaction.  If  a quick-phase  occurs,  is 
it  possible  to  make  a saccade  within  50  milliseconds?  Probably 
not.  If  a saccade  occurs,  is  it  possible  to  make  a quick-phase 
within  200  milliseconds?  This  is  not  known  because  the  situation 
has  not  been  experimentally  studied. 


As  for  the  locations  of  noises,  that  of  x (t)  is  not 

2 

altered,  while  that  of  z (t)  is  shifted  to  a new  location 

2 

(Figure  8.1).  The  noise  x (t)  is  not  involved  in  the  saccadic 

2 

mode  because  the  WHEN  function  is  bypassed.  Since  saccades 

are  accurate,  not  noisy,  the  noise  z (t)  is  shifted  to  a new 

2 

location  as  shown  in  Figure  8.1  so  that  it  is  associated  only 
with  the  VOR  mode.  Conceptually,  we  may  regard  the  portion 
enclosed  by  uhe  dotted  lines  in  Figure  8.1  as  the  "New  WHERE 
Function." 

A basic  problem  that  is  made  clear  by  Figure  8.1  is: 

What  is  the  difference  between  a saccade  and  a quick-phase? 

Is  it  possible  to  say  with  certainty,  this  is  a saccade, 
that  is  a quick-phase?  Four  possible  situations  for  rapid 
eye  movements  seem  to  exist  in  this  regard.  They  are: 
active  head  rotation  in  the  light;  active  head  rotation  in 
the  dark;  passive  head  rotation  in  the  light;  and  passive 
head  rotation  in  the  dark. 

It  is  our  opinion  that,  regardless  of  the  light-dark 
condition  under  which  a rapid  eye  movement  is  made,  the 
criteria  (for  whether  a rapid  eye  movement  is  a quick-phase 
or  a saccade)  should  indicate  whether  it  is  caused  by  a 
passive  head  rotation  or  by  the  active  acquisition  of  the 
target.  For  the  former  case,  it  is  a quick-phase.  For  the 
latter  case,  it  is  a saccade.  Thus,  in  a coordinated  head-eye 
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movement  to  a seen  target  in  the  light  or  an  auditory, 
felt,  or  even  remembered  target  in  the  dark,  the  movement  is 
a saccade.  It  is  a saccade  whether  the  eye  moves  before 
or  during  the  active  head  movement.  Conversely,  if  the  head 
is  passively  rotated,  the  rapid  eye  movement  is  a quick-phase. 

In  the  dark  that  is  by  definition,  but  even  in  the  light  we 
suggest  that  if  there  were  no  conscious  attempt  to  acquire  a 
target,  then  the  movement  should  be  called  a quick-phase. 

This  idet.  underlies  Figure  8.1,  because  only  when  higher 
perceptual  levels  are  trying  to  acquire  a target  are  the 
switches  thrown. 

8 . 2 Discussion  on  head  velocity  around  zero  value. 

For  head  velocities  near  zero  (e.g.,  less  than  10  deg/ 
sec),  we  do  not  know  much  about  the  behavior  of  our  cat's 
eye  movements.  First  of  all,  we  do  not  know  what  the  threshold 
value  Q should  be.  In  fact,  when  the  head  velocity  is  zero, 
the  question  becomes  almost  meaningless.  Also,  we  do  not 
know  how  to  handle  the  noise  level  around  zero  velocity  since 
it  can  cause  backward  "quick-phases"  when  the  noisy  WHEN  and 
WHERE  curves  intersect  each  other.  Here  again,  the  terms 
quick-phase  and  backward  become  meaningless  when  the  velocity 
is  zero.  In  the  absence  of  supporting  data,  we  have  arbitrarily 
picked  a value  for  Q (6  deg)  and  allowed  the  level  of  the  noise 
at  low  velocities  to  approach  zero  to  try  and  fit  the  sinusoidal 
behavior  particularly  in  the  region,  at  f = 0.05  Hz,  where  head 
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velocity  passes  through  zero.  The  main  problem  there  was 
backward  quick-phases.  For  this  reason,  if  we  are  to  extend 
the  research,  we  would  investigate  the  cat's  VOR  behavior 
at  low  velocities  more  thoroughly. 

Our  model  is  clearly  not  designed  for  operation  around 
zero  head  velocity.  One  alternative  is  to  suggest  that  the 
VOR  system  is  simply  disconnected  when  the  head  is  stationary; 
that  is,  in  some  neighborhood  of  zero  head  velocity.  One  way 
to  inclement  this  is  to  always  throw  the  switches  into  the 
saccadic  mode  below  some  level  of  4i(t)  in  Figure  8.1. 

Although  it  is  poorly  designed  for  operation  in  the 
neighborhood  of  zero  head  velocity  because  of  the  lack  of 
physiological  data,  our  VOR  model,  as  it  stands,  would  behave 
for  a zero  canal  signal  in  the  following  way.  The  noise  on 
the  WHEN  and  WHERE  curves  is  zero  (refer  to  Figure  7.5). 

The  WHERE  curve  is  located  at  zero  (section  5.3  and  Figure  5.8), 
and  the  WHEN  curves  are  located  at  +6  deg  (section  5.4  and 
Figure  5.9).  If  the  initial  position  of  the  eye  is  within 
deg,  the  eye  will  remain  there  as  long  as  the  canal  signal 
remains  zero.  If  the  cat  makes  a saccade  to  a point  beyond 
the  +6  deg  boundary  and  then  loses  "interest"  in  this  real  or 
imagined  target  so  that  the  switches  in  Figure  8.1  return  to 
the  \0R  mode,  the  eye  will  return  to  the  WHERE  curve  (zero), 
in  this  case  by  a quick-phase,  and  remain  there.  This  is 


not  unreasonable  behavior  but,  again,  we  come  to  the  problem  of 
what  distinguishes  a saccade  from  a quick-phase.  The  former 
is  "voluntary"  or  "active"  and  assumes  the  cat  is  "interested," 
all  of  which  are  subjective  terms  and  hard  to  deal  with  in  a 
model.  Thus,  even  if  more  data  were  available,  eye  movement 
records  at  very  low  head  velocities  would  bo  filled  with  rapid 
eye  movements  and  it  would  be  difficult  or  imijossible  to 
decide  which  ones  were  saccades  and  which  ones  were  quick-phases. 

8 . 3 Comments  on  other  VOR  models. 

To  our  knowledge,  there  exist  three  other  VOR  models 
concerned  with  quick-phases.  They  are  by  Sugie  and  Melvill  Jones 
(1971),  Schmid  and  Lardini  (1976),  and  Barnes  (1977).  The 
latter  two  models  were  published  while  this  research  was 
in  progress. 

The  model  by  Sugie  and  Melvill  Jones  was  the  first  published 
model  that  drew  attention  to  the  anti-compensatory  nature  of 
quick-phases.  To  compare  it  with  our  model,  its  essential 
features  are  shown  in  Figure  8.2.  In  this  figure,  our  present 
nomenclature  has  been  used;  the  canal  output  <))(t),  the  pulse 
P(t),  and  the  eye  position  9(t).  It  uses  periodic  sa.mpling 
every  200  milliseconds  to  initiate  quick-phases  which,  of 
course,  is  a bad  oversimplification.  The  transfer 
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Figure  8.2 

This  figure  shows  the  essential  features 
of  the  VOR  model  by  Sugie  and  Melvill  Jones 
(1971)  with  associated  waveforms.  For  a 
step  input  it  would  produce  the 

output  0(t)  as  shown. 
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Figure  8.2  The 


T 

function  ^ corresponds  to  our  neural  integrator.  They 

3 

used  a time  constant  T = 1 sec.  This  time  constant  was 

3 

determined  from  cats  under  light  ether  anesthesia.  In  Fig- 
ure 5.12,  the  element  in  our  model  most  susceptible  to  cuies- 
thesia  is  the  neural  integrator.  Anesthesia,  or  even  loss 
of  alertness,  causes  its  time  constant  to  decrease.  It  is 
generally  agreed  now  that  1 sec  is  a correct  value  for  the 
anesthetized  situation,  but  does  not  apply  to  alert  animals. 
This  short  time  constant  of  1 sec  would  cause  the  eye  position 
to  decay  exponentially  after  the  completion  of  a quick-phase 
so  that  a quick-phase  itself  creates  a slow-phase  component 
unrelated  to  <j>(t).  This,  in  fact,  became  the  whole  point  of 
their  model.  It  now  appears  that  it  was  based  on  a situation 

that  applies  only  under  anesthesia.  In  our  model,  the  time 

T 

constant  T^  of  the  neural  integrator,  ^ , is  25  sec. 

n 

Thus  in  essence,  it  acts  like  a pure  integrator.  In  addition, 
this  model  did  not  deal  with  the  stochastic  behavior  of  VOR. 

A sinplified  version  of  the  VOR  model  by  Schmid  euid 
Lardini  (1976)  is  shown  in  Figure  8.3.  Variables  analogous 
to  thost;  in  our  model  such  as  ())(t),  C(t),  6(t),  A0(t),  and 
P(t)  have  buen  so  labelled.  The  output  C(t)  of  the  transfer 
function  0 corresponds  to  our  deterministic  WHERE  curve  except 
that  a is  a scalar.  This  means  that  C(t)  is  linearly  propor- 
tional to  4>(t)  and  it  does  not  saturate  like  ours.  This,  of 
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course,  contradicts  with  experimental  observations  because 
it  does  saturate  in  reality.  If  A6  exceeds  +Q,  a quick-phase 

e~Ts 

IS  triggered.  is  a separate  neural  integrator  with  a 

delay  of  T = 0.2  sec  for  a refractory  period  although  it  can  bo 
as  low  as  50  milliseconds,  which  is  not  a trivial  matter. 

Figure  8.3  is  misleading.  What  actually  happens  in  the  quick- 
phase  mechanism,  when  |AG|  is  bigger  than  Q,  is  that  a sample 
is  taken.  Then,  the  sampler  is  inhibited  for  200  milliseconds. 
Then,  it  is  enabled,  so  that  another  sample  can  be  taken  if 
|A6|  is  bigger  than  Q.  That  is  what  Schmid  and  Lardini  meant 
to  depict.  This  integrator  is  in  addition  to  another  neural 
integrator,  , for  the  slow-phase  with  a time  constant 

of  150  sec.  Although  this  model  has  two  neural  integrators, 
there  are  several  reasons  both  theoretical  and  neurophysiologi- 
cal why  both  quick-phase  and  slow-phase  eye  movements  must 
shaire  the  same  neural  integrator  as  summarized  below  (Robin- 
son, 1975)  : 

A.  This  is  a theoretical  argument.  There  is  no  reason 
to  suppose  a saccade  differs  (from  a quick-phase  (Ron  et  al . , 
1972).  The  best  evidence  for  this  is  that  all  neurons  seen 
so  far  throughout  the  brain  whicJi  discharge  in  a saccade- 
correlated  way  do  the  same  thing  for  quick-phases  in  the  dark. 
Now,  during  rotatory  nystagmus  in  the  dark,  if  there  were 
separate  integrators,  then  (see  0^.  and  0^  in  Figure  8.3) 
one  would  accumulate  all  the  slow-phases  in  one  direction. 
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SAMPLER 

CONTROLLER 


Figure  8.3  Simplified  Version  of  the  VOR  Model  by  Schmid  and  Lardini  (1976) 


the  other,  all  the  quick-phases  in  the  other.  After  a short 
period  of  time,  each  integrator  would  have  stored  within  it 
the  equivalent  of  many  complete  revolutions  of  the  eye  in  the 
head.  This  unlikely  idea  can  be  dismissed  simply  on  con- 
ceptual grounds.  It  is  quite  clear  that  each  quick -phase 
pulse  must  reset  the  neural  integrator  to  a discharge  rate 
proportional  to  the  new  eye  position.  But  that,  of  course, 
is  equivalent  to  putting  the  pulse  into  the  integrator  as 
shown  in  Figure  4.1  and  Figure  5.12. 

B.  Neurons  which  behave  in  the  unlikely  fashion  just 
described  (the  two  integrator  hypothesis)  have  never  been 
observed  in  the  brains  of  alert  monkeys. 

C.  Cerebellar  lesions  in  cat  (as  well  as  in  man  and 
monkey)  make  the  eye  position  at  the  end  of  a saccade  slip 
back  with  a time  constant  of  1.3  sec.  It  turns  out  that  for 
the  same  cat  the  slow-phase  of  the  nystagmus  is  also  leaky  (does 
not  have  a constant  velocity)  and  velocity  falls  off  with  the 
same  time  constant  of  1.3  sec.  It  is  unlikely  that  this 

is  a coincidence  and  the  simplest  explanation  is  that  there 
is  only  one  integrator  now  made  leaky  by  cerebellectomy . 

D.  Finally,  if  a single  integrator  was  responsible  for 
creating  all  version  eye  movements,  a lesion  that  destroyed 
it  would  abolish  all  such  movements.  It  is  well  known  in 
neuro-ophthalmology  that  the  pprf  is  just  such  a location. 
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In  addition,  this  model  does  not  deal  with  the  stochastic 
behavior  of  the  VOR.  A comment  on  the  sampling  system  is 
given  later  in  this  section. 

A VOR  model  by  Barnes  (1977)  is  shown  in  Figure  8.4. 

Again,  we  have  substituted  our  own  terminology  to  help  in 
comparing  this  model  to  our  own.  The  output  C(t)  of  the 
transfer  function  A corresponds  to  our  deterministic  WHERE 
curve.  But,  unli)ce  ours,  it  does  not  saturate,  which  is 
contrary  to  observed  behavior.  Since  lc(t)  - 0^^(t)|  is 
the  input  to  the  threshold  device  T,  the  quantity  d"  more 
or  less  corresponds  to  our  efference  copy  of  the  eye  position 
(9^  in  Figure  5.12),  while  the  transfer  function  F seems  to 
convert  velocity  to  position  (more  or  less) , and  so  corresponds 
to  our  neural  integrator.  The  transfer  function  B provides 
a phase-lead  to  compensate  for  the  phase-lag  caused  by  the 
eye  dynamics.  This  could  easily  be  traded  off  to  a direct 
path  around  the  neural  integrator  as  in  Figure  3.1  and  Fig- 
ure 5.12.  The  discharge  pattern  similar  to  the  output  signals 
P(t)  and  Q"  (t)  from  B and  F are  not  observable  among  a variety 
of  brainstem  neurons  (see,  for  example.  Figure  5 of  Robinson, 
1975) , while  in  our  model  all  internal  signals  are  consistent 
with  the  neurophysiological  observations.  This  model  has  two 
neural  integrators  if  we  view  F as  an  integrator,  in  this  case 
both  are  reset  by  quic)c-phases.  Thus,  one  can  only  object  to 
two  integrators  on  the  grounds  of  parsimony.  In  addition,  this 
model  does  not  deal  with  the  stochastic  behavior  of  the  VOR. 
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Figure  8.4 

This  figure  shows  the  VOR  model  by 
Barnes  (1977).  Q"{t),  the  output 

of  F,  corresponds  to  the  efference 
copy  of  the  eye  position  Q' {t)  in  our 
model  (Figure  5.12),  and  T corresponds 
to  our  high-gain  amplifier  in  the 
bang-bang  control  system  of  our  model. 


( 
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In  summary,  if  the  goal  of  the  VOR  model  is  to  represent 
real  neurophysiology,  we  believe  our  model  is  more  realistic 
than  others  for  the  following  reasons. 

First,  our  model  deals  with  the  stochastic  behavior  of  the 
VOR,  while  the  other  three  models  do  not.  This  is  important 
simply  because  the  cat's  real  data  is  noisy,  and  the  noise 
statistics  are  an  important  part  of  the  data  base. 

However,  a really  significant  criteria  of  the  merit  of 
the  model  is  whether  the  internal  signals  of  the  model  are 
similar  to  the  neural  discharge  patterns  observed  in  the 
brain.  The  internal  signals  in  our  model  do  correspond  to 
discharge  patterns  of  neurons  in  the  brainstem.  In  other 
models,  this  is  not  the  case  as  explained  before  in  connection, 
for  excunple,  with  P(t)  and  d"  (t)  in  Figure  8.4.  The  model 
by  Schmid  and  Lardini  (Figure  8.3)  and  the  model  by  Barnes 
(Figure  8.4)  both  have  two  neural  integrators,  which  is  not 
correct  for  reasons  given  before.  The  model  by  Sugie  and 
Melvill  Jones  (Figure  8.2)  and  the  model  by  Schmid  and 
Lardini  (Figure  8.3)  both  use  a sampler  which  is  an  open-loop, 
ballistic  pulse  generator,  instead  of  a closed-loop  bang-bang 
controller  (pulse  generator)  as  in  our  model.  We  have  explained 
in  Chapter  4 why  we  chose  the  latter  over  the  former.  Further- 
more, the  sampling  system  has  inherent  problems  associated 
• with  it  because,  being  open-loop  and  ballistic,  it  has  to 
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specify  the  amplitude  of  the  sampler  output  so  that  the  end 

. point  of  quick-phase  reaches  the  predetermined  goal.  If  one 

{ 

now  adds  noise  in  the  system,  in  particular  to  the  threshold 
device,  it  becomes  complicated  and  a little  unnatural  as  to 
how  to  realize  a sampler  with  neurons  that  will  respond  to 
the  noise  in  an  approximate  manner.  The  closed-loop  continuous 
controller  avoids  these  problems. 

Of  course,  there  are  similarities  between  our  model  and 
the  other  models.  It  is  generally  accepted  that  quick-phases 
are  made  by  a pulse,  which  therefore  necessitates  an  inte- 
grator. The  pulse  must  be  triggered,  so  some  form  of  threshold 
i j device  is  needed  (save  for  the  Sugie  and  Melvill  Jones  model) 

into  which  some  time  varying  signal  must  go.  Almost  any  model 
would  of  necessity  have  such  elements.  All  modelers  since 
1971  have  recognized  the  fact  that  quick-phases  carry  the  eye 
in  the  anti-compensatory  half  of  the  motor  field  and  all  have 
used  a signal  from  the  canals  [C(t)]  to  effect  this.  Such  j 

elements  are  fairly  basic  and  it  is  not  surprising  that  they 
are  common  to  all  models  in  one  form  or  another. 

Hov^ever,  the  model  presented  here  does  differ  in  other  ways 
than  just  circuit  details.  The  approach  has  been  teleological. 

The  question  has  always  been  what  was  nature  trying  to  accom- 
•«  plish  in  designing  the  system,  and  then,  what  is  the  simplest 

I 

neural  circuit  to  achieve  that  goal?  This  approach  is  sub- 
jective and  disliked  by  some  but  to  us  it  has  more  appeal  than 
simply  wiring  together  black  boxes  to  fit  a data  base. 
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Probably  the  largest  difference  between  this  study  and 
others  is  the  quantitative  testing.  Other  modelers  are  often 
content  to  show  one  or  two  pictures  of  nystagmus  created  by 
the  model  and  comment  on  how  similar  they  looked  to  physio- 
logical records.  While  our  data  base  was  not  as  large  as  we 
realized  later  it  should  have  been,  the  analysis  of  both  cat 
and  model  behavior  is  much  more  extensive  than  anything 
attempted  before.  Thus,  this  model,  despite  its  obvious 
imperfections,  is  much  closer  quantitatively  to  physiological 
behavior  over  a wide  range  of  vestibular  stimuli  than  is  the 
case  for  other  models  so  far. 


231 


BIBLIOGRAPHY 


Barnes,  G.  R.  (1977)  The  Generation  of  Saaoadio  Eye  Movements 
in  VestihuXar  Nystagmus , RAF  Institute  of  Aviation  Medicine, 
Great  Britain  (in  press) . 

Bec)cer,  w.  and  Klein,  H.  M.  (1973)  Aaouracy  of  Saooadia  Eye 
Movements  and  Maintenance  of  Eooentria  Eye  Positions  in  the 
Dark,  vision  Res.  13,  1021-1034. 

Bryson,  A.  E. , Jr.  and  Ho,  Y.  C.  (1975)  Applied  Optimal 
Control,  John  Wiley  & Sons,  New  Yorlc,  NY. 

Cheng,  M.  (1972)  Point  Process  Analysis  in  the  Physiological 
Study  of  Human  Nystagmus , Ph.D.  Thesis,  McGill  University, 
Montreal,  Canada. 

Cohen,  B.  and  Henn,  V.  (1972)  The  Origin  of  Quick-Phases  of 
Nystagmus  in  the  Horizontal  Plane,  Bibl.  Ophthal.  82,  36-55. 

Fernandez,  c.  and  Goldberg,  J.  M.  (1971)  Physiology  of  Periph- 
eral Neurons  Innervating  Semicircular  Canals  of  the 
Squirrel  Monkey:  II  Response  to  Sinusoidal  Stimulation  c — ^ 

Dynamics  of  Peripheral  Vestibular  System,  J.  Neurophysic  .,gy 
34,  661-684. 

Fuchs,  A.  F.  and  Kimm,  J.  (1975)  Horizontal  Angular  Acceler- 
ation and  Eye  Movement,  J.  Neurophysiology  38,  1140-1161. 

Fuchs,  A.  P.  and  Luschei,  E.  S.  (1970)  Firing  Patterns  of 
Abducens  Neurons  of  Alert  Monkeys  in  Relationship  to 
Horizontal  Eye  Movement,  J.  Neurophysiology  33,  382-392. 

Fuchs,  A.  F.  and  Robinson,  D.  A.  (1966)  A Method  for  Measuring 
Horizontal  and  Vertical  Eye  Movement  Chronically  in  the 
Monkey,  J.  Appl.  Physiology  21,  1068-1070. 

Hallet,  P.  E.  and  Lightstone,  A.  D.  (1976)  Saccade  to  Flashed 
Targets,  Vision  Res.  16,  107-114. 

Jurgens,  R.  and  Bec)cer,  w.  (1975)  Is  There  a Linear  Addition 
of  Saccades  on  Pursuit  Movement?  Basic  Mechanisms  of  Ocular 
Motility  and  Their  Clinical  Implications,  Pergamon  Press, 
Oxford  and  New  Yor)c. 


PFECEDIlfi  PAGE  PLANK-NOT  PILHCD 


i- 


K 


KflK't,  t.  L.  <l'i74)  :'u-ticipation  of  the  Medial  Pontine 
,r»- \n  Fife  Movement  Generation  in  Monkey, 

J.  Nfuro^ihys  ioioi^y  37,  316-332. 

Keller,  L.  I,,  and  Robinson,  D.  A.  (1971)  Absence  of  a Stretch 
refl^jr  'n  Estraocular  Muscles  of  Monkey,  J.  Neurophysiology 
34,  908-919. 

.Melvill  Jones,  G.  (1964)  Predominance  of  Anti-Compensatory 
Oculomotor  Response  During  Rapid  Head  Rotation,  Aerospace 
Medicine  35,  965-968. 

Melvill  Jones,  G.  and  Milsum,  J.  H.  (1965)  Spatial  and 
Dynmic  Aspects  of  Visual  Fixation,  IEEE  Trans.  Bio-med. 

Eng.  BME  12(2),  54-62. 

Melvill  Jones,  G.  and  Milsum,  J.  H.  (1971)  Frequency  - 
Response  Analysis  of  Central  Vestibular  Unit  Activity 
Resulting  from  Rotational  Stimulation  of  the  Semicircular 
Canals,  J.  Physiology  218,  1-25. 

Morasso,  P.  (1973)  Adjustment  of  Saccade  Characteristics  During 
Head  Movements,  Experimental  Brain  Research  16,  492-500. 

Outerbridge,  J.  s.  (1969)  Experimental  and  Theoretical 
Investigation  of  Vestibular-Driven  Head  and  Eye  Movement, 

Ph.D.  Thesis,  McGill  University,  Montreal,  Canada. 

Papoulis,  A.  (1965)  Probability,  Random  Variables,  and  Stochastic 
Processes , McGraw-Hill  Boo)c  Company,  New  Yorlt,  NY. 

Robinson,  D.  A.  (1963)  A Method  of  Measuring  Eye  Movement 
Using  a Scleral  Search  Coil  in  a Magnetic  Field,  IEEE 
Transaction  on  Bio-medical  Electronics  10,  137-145. 

Robinson,  D.  A.  (1970)  Oculomotor  Unit  Behavior  in  the 
Monkey,  J.  Neurophysiology  33,  393-404. 

Robinson,  D.  A.  (1973)  Models  of  the  Saccadic  Eye  Movement 
Control  System,  Kybernetic  14,  71-83. 

Robinson,  D.  A.  (1975)  Oculomotor  Control  Signals,  Basic 
Mechanism  of  Ocular  Motility  and  Their  Clinical  Implications, 
Pergeimon  Press,  Oxford  and  New  Yorlc. 

Robinson,  D.  A.  (1976)  Adaptive  Gain  Control  of  Vestibuloocular 
Reflex  by  the  Cerebellum,  J.  Neurophysiology  39,954-969. 


Robinson,  D.  A.  and  Keller,  K.  L.  (l'J72'  The  Behavior  of  Eye 
Movement  Motoneuron  in  * hi  Ah'rt  M mh  _i , ‘'ihl.  Ophthal.  82, 

7-16. 

Ron,  S. , Robinson,  D.  A.  , jnd  Skavenski,  A.  A.  (1972)  Saoaodes 
and  the  Quick  Phase  nf  Hystaamus , Vision  Res.  12,  2018-2022. 

Scheuer,  E.  M.  and  stoller,  D.  s.  (1962)  On  the  Generation  of 
Normal  Random  Vectors,  Toclmometrics  4,  278-281. 

Schmid,  R.  and  Lardini,  F.  (1976)  On  the  lOredominance  of 
Anti-Corrq'ensatory  Eye  Movements  in  Vestibular  Nystagmus , 

Biol.  Cybernetics  2b,  135-148. 

Shimazu,  n.  and  Precht,  w.  (1966)  Inhibition  of  'cntral 
Vestibular  Neurons  from  the  Contralateral  Labyrinth  and  Its 
Mediating  Pathway , J.  Neurophysiology  29,  467-492. 

Skavenski,  A.  A.,  Haddad,  G. , and  Steinman,  R.  M.  (1972)  The 
Extraretinal  Signal  for  the  Visual  Perception  of  Direction, 
Perception  and  Psychophysics  11,  287-290. 

Skavenski,  A.  A.  and  Robinson,  D.  A.  (1973)  Hole  of  Abducens 
Neurons  in  the  VestibuloocuZar  Reflex,  J.  Neurophysiology  36, 
724-738. 

Sparks,  L.  L.  and  Travis,  R.  P. , Jr.  (1971)  Firing  Patterns  of 
Reticular  Formation  Neurons  During  Horizontal  Eye  Movements, 

Brain  Res.  33,  477-481. 

Sugie,  N.  and  Melvill  Jones,  G.  (1971)  A Method  of  Eye  Movements 
Induced  by  Head  Rotation,  IEEE  Transactions  on  Systems,  Man,  and 
Cybernetics,  Vol.  SMC  1(3),  251-260. 

Young,  L.  R.  (1969)  The  Current  Status  of  Vestibular  System 
Models,  Automatica  5,  369-383. 

Young,  L.  R.  and  Stark,  L.  (1963)  Variable  Feedback  Experiments 
Testing  a Sample  Data  Model  for  Eye  Tracking  Movements , IEEE 
Trans  of  the  Prof.  Tech.  Grp.  on  Human  Factors  in  Elect.  HPE  4, 
38-51. 

Zee,  D.  S.  Optican,  L.  M. , Cook,  J.  D. , Robinson,  D.  A.,  and 
Engel,  w.  K.  (1976)  Slow  Saccades  in  Spinocerebellar  Degeneration, 
Arcliives  of  Neurology  33,  243-251. 


r 


DISTRIBUTION 


Library  of  Congress 
Washington,  U.  C.  20540 

Attention:  (lift  and  Exchange  Division  (4) 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  (22 

Strategic  Systems  Project  Office 
Department  of  tlie  Navy 
Washington,  D.  C.  20360 
Attention:  D.  Gold  (SP-230) 

LCDR  Ferrlter  (SP-23110) 

General  Electric  Ordnance  Systems 
100  Plastics  Avenue 
Pittsfield,  Massachusetts  01201 
Attention:  B.  Flood 

P.  Schubert 

Charles  Stark  Draper  Laboratory 
555  Technology  Square 
Cambridge,  Massachusetts  02139 
Attention:  P.  Howard 

C.  Hrbek 
K.  Vincent 

The  Analytic  Sciences  Corporation 
Six  Jacob  Way 

Reading,  Massachusetts  01867 
Attention:  J.  Bortz 

J.  Kasper,  Jr. 

Chief  of  Naval  Researcii 
Department  of  the  Navy 
Arlington,  Virginia  22217 

Applied  Pliysics  Laboratory 
Johns  Hopkins  University 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Department  of  Ophthalmology 
School  of  Medicine 
Johns  Hopkins  University 
Baltimore,  Maryland  21205 

Attention:  Professor  D.  A,  Robinson  (9) 


' •'  *.  . . , V . 


fiimed 


Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technolog/ 
Cambridge,  Massachusetts  02139 
Attention:  Professor  Yao  Tzu  LI 

Professor  Laurence  R.  Young 

National  Aeronautics  and  Space  Administration 
Washington,  D.  C.  20547 

Attention:  Technical  Library 

National  Research  Institute 
National  Naval  Medical  Center 
Bethesda,  Maryland  20014 

De^'ense  Printing  Service 
Washington  Navy  Yard 
Washington,  D.  C.  20374 

Depart \ent  of  Electrical  Engineering 
Vlrgir . i Polytechnic  Institute  and 
Stai  Univei jlty 
Blacksburg,  Virginia  24061 

Attention:  Professor  Van  Landingham 


Local : 

CC 

DK-56 

CD 

DK-60 

DD 

DK-70 

DK-71 

DF 

DN 

DF-30 

DN-23 

DC 

DT 

DC- 40 

WD 

DK 

WA 

DK-01 

I7E 

DK-02 

WR 

DK-05 

WU 

DK-10 

(2) 

DK-20 

DX-21 

DK-30 

(2) 

DX-222 

DK-40 

DX-40 

DK-50 

(3) 

DK-51 

DK-53 

DK-54 

DK-55 

DK-55  (Chun) 

(35) 

(2) 

(6) 


